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ABSTRACT  fC*nr(nu*  «n  rmvttit  i id*  it  nictiim  an4  idrnitty by  Woe*  nvmbtt) 

Representations  of  space/time  signals  as  stationary,  homogeneous  random  processes 
are  used  to  formulate  the  theory  of  optimum  array  processing.  Generalized  Founer  expan- 
sions for  one,  two,  and  three  dimensions  are  discussed.  These  representations  include  direc- 
tional, isotropic  signal  fields,  as  well  as  all  previously  published  ambient  noise  models  Re- 
ceiving arrays  are  modeled  in  terms  of  how  they  observe  signal  fields  with  the  roles  of  their 
geometry  and  the  observation  noise  analyzed.  The  design  of  optimum  beamformers  is 


1473  COITION  Of  1 NOV  SS  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  Of  THIS  PAGI 


UNCLASSIFIED 


-UUlTY  CLASSIFICATION  OF  THIS  PAGEftOlui  Ef“*r»d) 


20  (Continued) 

formulated  and  coupled  to  the  representation  of  the  signal  field  and  the  receiving  array 
The  structure  of  optimum  beamforms  are  analyzed  and  related  to  beam  displacement,  null 
placement,  and  superdirective  methods  Processing  performance  obtained  by  exact  methods 
and  approximate  methods  are  compared  for  various  sign'd  fields  and  array  geometries  The 
closeness  of  the  rerults  indicate  the  usefulness  ot  the  representations  and  array  analysis  in 
predicting  the  performance  of  optimum  systems  Finally,  the  structure  of  preformed  beams 
and  clustered  receiving  arrays  are  analyzed  with  these  methods. 


iii 


CONTENTS 


i 


i i. 

2. 

I 

3. 

4 

5. 

6. 
7. 
8 


Part  1 Characterizing  Space/Time  Random  Processes 
Introduction  1 

Space/Time  Random  Processes  2 

Part  2 Response  of  Arrays  to  Various  Noise  Fields 

Receiving  Apertures  40 

Optimum  Array  Processing  for  Plane  Wave  Signals  64 
Part  3 Processing  Methods 
Linear  Arrays  78 

Optimum  Array  Processing  for  Directional  Noise  Fields  130 
Array  Geometries  in  Higher  Dimensions  133 
Array  Processing  Systems  1 98 


ill 


c3 


ti 


PART  1 - CHARACTERIZING  SPACE/TIME  RANDOM  PROCESSES 


1.  INTRODUCTION 

The  analysis  of  arrays  of  transducers,  or  antennas,  has  had  such  an  extensive  develop- 
ment by  now  thai  one  anticipates  little  which  is  novel  Yet,  m spite  of  the  widespread  use 
of  statistical  methods  in  communication  systems,  it  is  compantively  recent  that  array,  or 
antenna,  processing  methods,  based  on  the  statistical  properties  of  both  the  signal  and  the 
mt  Tn  mg  noise,  were  introduced  However,  sufficient  time  has  elapsed  such  that  a rather 
formidable  literature  is  available 

Unfortunately,  most  contributions  to  this  literature  are  quite  mathematically  oriented— 
usually  couched  in  a maze  of  covariance  matrix  manipulations  While  a certain  amount  of  this 
may  be  necessary  to  formally  establish  the  logic  of  the  approach,  the  underlying  principles 
of  operation  are  quite  intuitive  and  can  be  theoretically  explained  in  terms  of  an  appropriate 
combination  of  classical  antenna  theory  and  some  fairly  straightforward  spatial  random  pro- 
cessing Development  of  this  intuitive  understanding  is  particularly  important  since  the  pro- 
cessing gam,  or  signal-to-noise  ratio,  often  carries  a high  premium  for  normal  array  processing 
applications  Often,  direct  implementation  of  “optimum”  processing  methods  is  not  possible 
either  from  computational  or  physical  considerations  Consequently,  the  typical  concern  is 
with  suboptimum  methods  which  exploit  the  essential  features  of  the  “optimum”  processing 
to  obtain  all  but  "epsilon"  of  the  ultimate  performance  In  addition  to  understanding  how 
spatial  processors  obtain  their  optimum  performance,  it  is  important  to  determine  just  what 
this  performance  is  since  it  sets  the  standard,  or  benchmark,  against  which  to  compare  the 
performance  of  any  suboptimal,  but  implementable, stem. 

The  use  of  an  array  of  transducers,  or  an  antenna,  is  typically  motivated  by  two 
reasons,  both  of  which  are  generally  coupled  together.  First,  the  particular  spatial  coherency 
of  a propagating  signal  is  exploited  to  improve  the  signal  level  versus  the  background  noise 
The  dominant  features  governing  the  degree  of  improvement  are  the  array  geometry  and  the 
spatial  structure  of  the  background  noise.  Sometimes,  the  signal  can  be  improved  with  simple 
statistical  smoothing,  while  at  other  times  the  coherence  of  the  background  can  be  used  to 
combat  background  noise.  In  the  second  application,  the  directional  characteristics  of  the 
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ambient  field  of  signals  must  be  determined  Again,  the  degree  to  which  a source  position  can 
be  located  is  influenced  by  the  same  factors  which  determine  its  detectability 

This  report  is  divided  into  three  parts  Part  I considers  the  problem  of  characterizing 
space/time  random  processes  The  study  of  this  problem  is  fundamental  to  the  analysis  of 
the  operation  of  arrays,  or  antennas,  with  random  inputs 

In  Part  2 we  describe  the  response  ol  airays  to  vauous  noise  fields  We  consider  the 
optimization  and  resulting  optimum  performance  of  processing  methods  ior  signals  and  noise- 
fields  of  interest 

Part  3 considers  a sequence  of  topics  based  upon  the  results  developed  in  Parts  I and 
2 These  focus  on  the  intuitive  concepts  of  the  processing  methods  in  the  context  of  some- 
selected  examples  of  noise  fields  and  array  geometries  and  on  the  introduction  of  some 
selected  suboptimum  processing  methods 

Two  comments  are  appropriate  here  First,  we  assume  tnat  the  reader  has  a background 
in  temporal  random  processes  and  a knowledge  of  Fourier  methods  Appropriate  references 
are  given  by  Papouhs  [1  ].  Davenport  and  Root  [2],  and  Van  Tree,  (3]  Second,  we  formulate 
our  arrays  in  terms  of  a continuum,  or  an  aperture,  and  a representation  of  spatial  signals 
observed  across  an  aperture  This  leads  to  a notation  which  is  somewhat  more  abstract  than 
that  which  corresponds  to  discrete  element  arrays  However,  it  does  lead  to  a more  funda- 
mental consideration  of  the  issues  in  spatial  processing  and  to  results  which  are  much  more 
intuitive  and  uncluttered  by  the  spatial  sampling  ,ssue  introduced  by  a discrete  formulation 
Once  these  results  are  established  with  a continuous  formulation,  we  consider  the  issues  needed 
to  relate  them  to  discrete  arrays  In  essence,  we  are  asserting  that  lor  many  situations  a contin- 
uous analysis  is  more  fundamental  and  easier,  although  we  recognize  that  ultimately  most 
arrays  of  interest  here  are  discrete  and  our  results  must  be  related  to  them 


2.  SPACE/TIME  RANDOM  PROCESSES 

In  the  theory  of  temporal  processes  there  are  several  ways  to  characterize  or  describe 
the  random  processes  of  interest  Ultimately,  one  would  want  to  have  a complete  specifica- 
tion of  the  processes  in  terms  of  their  probability  distribution  functions  for  an  arbitrary 
number  of  points  in  the  space/time  field  With  the  prominent  exception  of  Gaussian  and 
Poisson  processes,  this  information  is  generally  not  available  and  one  must  settle  for  a less 
complete  or  partial  characterization  Most  analyses  involve  a second  moment  characterization 
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in  which  one  specifies  the  mean  and  covariance  function  for  the  processes  of  interest  As  is 
welt  known,  this  information  is  alro  sufficient  to  provide  a complete  characterization  for 
Gaussian  processes  since  its  probability  distribution  of  arbitrary  order  can  be  related  to  these 
two  moments  We  pursue  an  approach  that  uses  the  second  moment  characterization  which 
suffices  for  most  of  our  applications  Further,  we  assume  that  the  processes  are  also  Gaussian 
when  we  ask  probability  questions,  e g , the  performance  ot  optimal  detectors 

A second  method  of  analysis  involves  a desenption  of  how  the  random  processes  are 
generated  While  this  has  found  extensive  use  for  purely  temporal  processes.  i>  euucnced  by 
the  wealth  of  material  in  the  literature  regarding  state  variable  estimation  procedures,  it  has 
not,  as  yet,  found  extensive  use  for  space/time  processes  The  role  of  these  processes  and  their 
relation  to  districuted  systems  is  clearly  an  emerging  field  of  interest,  however,  the  analysis  of 
the  partial  differential  equations  introduced  to  date  is  quite  complicated  and  is  not  yet 
particularly  fruitful 

2.1  SPACE/TIME  COVARIANCES  AND  FREQUENCY  WAVE  NUMBER  FUNCTIONS 

Let  us  assume  that  we  are  concerned  with  a signal  y(t  z)  that  is  defined  over  both  time 
and  space  1 We  model  y(t  z)  as  a random  process  defined  over  a field  with  a temporal  domain 
or  index  set  t e {TQ,  Tf(  and  a spatial  domain  or  index  set  ze  RN  with  a mean 

my(t  z)  = E[y(t.z)l  (2  1) 

an.'  a space/time  covariance 


Ky(t,T.z,f)  = E[(y(t  z)-my  (t  z))  (y(r  £)-my(r:£))  *J 


(we  often  assume  zero  mean  processes  for  simplicity) 

In  this  most  general  case  the  processes  are  both  temporally  non-^tationary  and  spatially 
non-stationary,  or  non-homogeneous  Analyzing  problems  that  involve  this  class  of  * ovananct- 
functions  using  a second  moment  characterization  introduces  many  issues  that  arise  with 

^ We  use  scalar  processes  and  introduce  vector  on "s  only  when  necessary 

lt~*Complex  Conjugate  We  use  a complex  formulation  throughout  This  u more  appropriate  for  a sonar  analysu  than  it 
is  for  seismic  data,  however,  the  correspondence  for  strati}  real  processes  is  direct  Sec  Ref  3 foi  a complete  development 
o( complex  process  representation 
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problems  involving  temporally  non-stationary  processes  In  general,  the  analyses  complicate 
the  mathematical  detail,  but  the  conceptual  issues  are  unchanged 

In  many  practical  situations  the  processes  encountered  are  temporally  stationary,  at 
least  m a wide  sense,  such  that  we  have 

my(t.z)  = itiy(z)  (2  3) 


Ky(t,T.z.£)  = K (t-r  ip 


This  assumption  may  not  be  strictly  satisfied  over  long  intervals  of  time  However,  when  the 
processing  interval  is  short  compared  to  the  time  period  over  which  the  covariance  changes,  it 
is  a realistic  practical  assumption 

We  also  frequently  encounter  situations  in  which  the  space/time  processes  are  homo- 
geneous or  spatially  stationary  The  coy '.nance  is  only  a function  of  the  vector  difference  of 
the  spatial  arguments 


m.,(t:z)  = mv(t) 


Ky(t,r  z.f)  = Ky(t-r  z-f)  (2  6) 

Again,  this  may  not  be  strictly  satisfied  over  large  distances,  although  it  is  a realistic  approxima- 
tion when  the  receiving  aperture  is  small  compared  to  the  distances  over  which  the  covariance 
is  non-homogeneous.  We  point  out,  however,  that  several  problems  which  involve  signals  in  the 
near  field  and  array  shading  effects  need  to  be  discussed  in  the  context  ot  a non-homogeneous 
covariance,  so  the  general  model  is  certainly  not  void  of  practical  interest. 

Finally,  we  have  the  situation  in  which  a stationary,  homogeneous  process  exists  Com- 
bining the  two  above  assumptions  we  have 


my(t.z)  = my 
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and 


Ky(t,r.z,£)  = Ky(t-rz-£)  (2  8) 

We  dnect  our  attention  principally  to  problems  involving  processes  which  fit  within  this 
general  class  of  assumption.  We  also  drop  the  mean  from  further  consideration 

In  many  of  our  analyses  it  is  useful  to  consider  a spectral  decomposifion  of  the  random 
processes  of  interest.  By  pursuing  an  analysis  in  a frequency  domain,  either  temporal  or 
spatial,  we  can  often  obtain  a significant  amount  of  insight  into  the  operations  of  our  pro- 
cessors. To  do  this  we  define  a series  of  Fourier  transform  operations  The  temporal  frequency 
spatial  correlation  function  associated  with  a stationary  process  is  given  by 


>:z,£)  - / Ky(r  z,£] 

too 


r)e-j“rd  T 


When  z = J,  this  yields  the  temporal  spectrum  at  a point  in  space.  If  we  consider  sampling 
at  a number  of  points  in  space,  so  as  to  form  a vector  process,  we  can  order  the  points 
and  create  a matrix  format  the  spectral  covariance  of  the  vector  process  by  sampling 
Sy(co:z,f),  appropriately. 

Similarly,  for  a homogeneous  process,  we  can  define  the  temporal  correlation,  spatial 
wave  number  function  as 


Fy(t,rk)  = 


noo 

dzKyft.rizJe1-'- 

O 


Note  that  the  integration  is  performed  over  the  entire  spatial  domain  of  the  process 
Consequently,  there  is  a fundamental  diffeience  among  the  transforms  over  one,  two,  or 
three  dimensions  We  do  not  use  this  function  very  frequently  because  there  are  better 
alternate  representations  for  the  signal  in  the  temporal  frequency  domain 
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For  <1  stationary  homogeneous  process,  the  frequency  wave  number  function  is 


defined  as 


Py(ork)  = dr  dz  Kytr'zle')^7^'^ 


(Note  the  sign  difference  with  respect  to  the  spatial  variable) 

This  function  has  a very  appealing  interpretation  When  we  pursue  an  analysis  similar  to 
that  used  for  the  spectral  representation  theorem  for  temporal  processes,  we  find  that  the 
process  y(t.z)  can  be  represented  in  a Sticltjes  integral  form  [4]  as 


.».rrr 

Joo  'loo  "-QO 


ej(wt-k-z)dY(w 


where  k has  dimension  N,  the  dimension  of  the  space  over  which  the  process  is  defined.  For 
the  decomposition  of  Y(w  k)  we  have  the  fol'owing  properties' 


E[|Y(w2i2)-Y(co1'i1)|2] 


Pis  f2Ji_ 

l 2"  i <WN 


P (co  k)  (2  13a) 


or  formally  E[ldY(w  k)l2l  = — — Pv(w:k). 

2jt  (2w)N  » 


E[(Y(ca2:k2)  - Y(coj  k, ))(Y(oo4  j^)  - Y(u3.k3))*)  = 0, 


v/hen  ((cj],w2)  x (k],k2))  n ((w3,co4)  x (k3,k4))  = 0, 


i.e  , disjoint  frequency  wave  number  bands  are  uncorrelatcd  Consequently,  wc  may  consider 
this  class  of  space/time  random  processes  to  be  composed  of  a superposition  of  plane  waves  of 
radian  frequency  u and  wave  number  J<  This  is  the  spatial  generalization  of  a spectral  decom- 
position of  a random  process 
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If,  in  addition,  we  impose  the  condition  that  the  process  y(t  z)  satisfies  the  homogeneous 
wave  equation,  then  we  have  the  constraint 


r« 


(2.15a) 


2jrf 


— = fX  = c (velocity  of  propagation) 

|k| 


(2  15b) 


This  defines  a structure  which  signals  propagating  in  a medium  must  satisfy  Further  this 
implies  that  at  ai.v  particular  cj  the  frequency  wave  number  function  is  nonzero  only  on  a 
sphere  of  radius  2ir/X  = |k|  = cj/c  This  constraint  is  not  necessary  for  y(t  t)  to  correspond 
to  a propagating  process  For  example,  one  may  have  propagation  in  three  dimensions  and  yet 
consider  a representation  over  a two-dimensional  surface  since  the  two-dimensional  wave 
number  value  corresponds  unambiguously  to  a three-dimensional  surface  except  for  its  sign 

The  inherent  advantage  of  this  representation  for  homogeneous  fields  is  that  it  allows  us 
to  make  many  statements  regardmg  the  processing  of  our  signals  which  are  completely  dual  to 
analogous  frequency  domain  operations  when  processing  temporal  waveforms  Naturally,  there 
is  a close  relationship  among  the  various  transforms  discussed,  which  are  summarized  in  Figure 
2 1 The  paper  and  text  by  Yaglom  [5,61  are  particularly  appropriate  references  for 
the  harmonic  analysis  of  space/time  processes.  With  these  references  the  rigor  in  using 
the  harmonic  analysis  can  be  pursued  at  length.  For  our  purposes  the  major  problem 
relates  to  the  notation  of  representing  y(t  z)  in  either  temporal  or  spatial  frequency  domains 
This  can  be  quite  awkward  if  pursued  too  far  We  shall  use  the  Stieltjes  notation  dY(co  k)  when 
we  wish  to  do  this  with  the  important  properties  of  Y(co  k)  summarized  in  Eqs  2 1 3 and 
2 1 4 above. 

Often,  it  is  convenient  to  use  different  frtquency  measures  The  temporal  frequency,  f,  in 
cycles  per  sec,  or  Hertz,  is  given  by 


f = oo/  2rr 


(2  16) 


Similarly,  the  spatial  wave  vector  v,  in  cycles  per  unit  d istance,  is  given  by 


v = k/2ir 


(2  17) 


& 

■lii 


v.v 


V*V 

AS 


i 


tv  'ji  'y+Fjr'y'jpy  *>  \v  *y*i  ■;*  W W.V*  V 

1$ 


. > v_VV  J-V--7;%>7.“7.-  -V_.  V. 


Figure  2-1  Relationship  among  the  various  second  moment  representations  for  stationary  - 
homogeneous  space  time  random  processes 


For  plane  waves 


S 

i 


|k/2jr|  =|n|  = f/c  = l/\ 

(2.18) 

It  is  also  convenient  to  define  the  etfective  wavelength  along  each  axis  in  lerms  of  the  com- 
ponents of  v.  In  Cartesian  coordinates  we  have 

"x  = ‘Ax 

(2  19a) 

•V-lAy 

(2.19b) 

V/~-  1 /X2 

(2  19c) 

(Note  that  one  does  not  obtain  these  effective  wavelengths  by  simply 

projecting  the  total  wave- 

length  oriented  along  the  propagation  direction  upon  the  respective  axes ) 

To  illustrate  these  concepts  and  to  develop  some  results  that  we  need  later,  we  present 
some  examples  First,  we  shall  consider  the  analysis  of  some  of  the  various  models  popularly 
used  in  the  literature.  Next  we  consider  a genera!  discussion  of  representing  space/time  processes 
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Here  we  also  draw  upon  some  results  from  electromagnetics  which  apparently  have  not  found 
their  way  into  the  sonar/seismic  literature 

Generally,  n is  easiest  to  discuss  these  space/timc  processes  in  terms  of  either  the 
temporal  frequency  spatial  coirelation  or  the  frequency  wave  number  function.  Consequently, 
many  results  are  specified  at  a single  temporj  frequency,  or  for  a narrow  freuuency  band  For 
a broadband  analysis,  the  integration  over  a specified  frequency  range 's  implied. 

Example  1 Directional  Signal 

The  simplest  signal  c f mteiest  is  a plane  wave  propagating  in  a direction  with  speed  c 
The  space/time  process  has  the  form 

y(t.z)  = y0[t-(a-i/c)],.  (2  20) 

as  illustrated  in  Figure  2 2 
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Figure  2-2  Model  of  plane  wave  propagation 


If  vve  assume  that  y0(t)  is  a stationary  random  process,  then  the  space  time  correlation 
function  is  given  by  (assuming  zero  means) 


Ky(Ar  Az)  = E[y(t  z)  y*(t-Ar  z-Az))  = E[y0(t-(a-z/c))  y0*(t-AT-(a-(z-Az)/c))] 

= Kyo(Ar-(a-Az/c))  (2.21) 

The  temporal  frequency  spatial  correlation  function  follows  the  Fourier  transform  relation 
of  Eq  2 9.  We  find 

S„(co.Az)  = Sv  (<o)e-J<"/c)te‘AJ)  (2  22) 

y - Vq 

Consequently,  we  have  cross-spectrum  of  a plane  curve  at  two  different  locations  which  is 
related  to  the  spectrum  of  the  process  at  a single  point  by  a simple  linear  phase  shift  that 
reflects  the  propagation  phase  between  the  two  points.  Therefore,  the  elements  of  spectral 
matrices  that  involve  pure  plane  waves  consist  of  a common  amplitude  and  a set  of  phase  shifts 
between  the  various  sensors 

If  we  now  transform  this  function  with  respect  to  the  spatial  variable,  we  obtain  the  fre- 
quency wave  number  function 

Py(w:k)=  JJ dzSy(co:z)e&'z  = Syo(w)u0^-^  aj  (2  23) 

RN  ’ 

where  u0(k-  j aj  = uQ(kx  - j ax)  u0(ky  - " ayj  u^kz  - j (2  24) 


is  a three-dimensional  impulse  function  (or  two-dimensional  for  a planar  analysis),  which  m 
Cartesian  coordinates  can  be  written  as  a product  of  impulse  functions  of  a single  variable.  (In 
discussing  spatial  impulse  functions,  one  needs  to  be  very  careful  in  defining  them,  especially 
with  regard  to  the  particular  coordinate  system  used  The  common  limiting  sequences  often 
yield  paradoxical  results,  and  an  operational  definition  that  uses  generalized  functions  is 
appropriate  [7] .) 


l*£2 
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We  find  that  a plane  wave  propagating  with  speed  c and  direction  a^has  an  impulsive  frequency 
wave  number  spectrum  located  at 


ro 

c 


a 


(2  25) 


in  wave  number  space.  (Unless  the  wave  is  monochromatic,  1 e , at  a single  temporal  frequency, 
the  magnitude  of  the  wave  number  vector  changes  as  a function  of  the  temporal  frequency  ) We 
again  point  out  that  the  space  over  which  this  transform  is  taken  must  be  specified  Theoret- 
ically, m a homogeneous,  isotropic  medium  the  magnitude  of  the  wave  number  k at  a 
particular  frequency  should  be  a constant,  indicating  that  for  three  dimensions  the  function 
Py(co  k)  is  nonzero  only  on  a sphere,  nonzero  cn  a circle  in  two  dimensions,  or  nonzero  at  two 
points  opposite  in  sign  in  one  dimension 

However,  in  many  cases  one  considers  representing  the  waves  in  a three-dimensional  space 
as  projected  on  a two-dimensional  geometry.  This  leads  to  an  analysis  in  which  the  magnitude 
of  the  wave  number  is  not  a constant  at  a particular  temporal  frequency.  It  should  be  obvious 
that  any  set  of  statistically  independent  directional  signals  propagating  as  plane  waves  can  be 
represented  with  an  impulsive  frequency  wave  number  function.  However,  such  a representation 
is  not  possible  if  the  components  are  correlated  as  could  possibly  be  envisioned  in  some  coherent 
multipath  situations  The  fundamental  difficulty  stems  from  the  requirement  that  components 
from  disjoint  regions  of  wave  number  space  must  be  uncorrelated  as  specified  by  Eq.  2 14, 
therefore  such  a process  would  not  be  homogeneous. 

Another  commonly  used  mode!  is  isotropic  noise.  This  noise  process  is  commonly 
advocated  for  ambient  sea  noise.  It  can  be  viewed  as  the  superposition  of  plane  waves  propagating 
from  all  directions  with  a uniform  statistical  level.  This  has  a spectral  covariance  structure  of 
the  form 

Sy( w z - f)  = S0(<o)  sine  |z  - £|\  (2  26)1 


Since  this  particular  process  fits  within  a more  general  context  of  a plane  wave  process,  we 
discuss  this  more  general  representation 


l 


a 

sinc(x)  = 

x 


^ j.  -j*  v r". --1  *-r -."  '- ' 


2.2  REPRESENTATIONS  FOR  PROPAGATING  SIGNAL  PROCESSES 

2.2.1  Processes  in  Three  Dimensions 

2.2.1. 1 Plane  Wave.  A commonly  employed  model  for  the  ambient  noise  background  is 
isotropic  or  omnidirectional  noise  We  first  describe  this  process  and  many  others  in  terms  of 
a model  for  generating  its  temporal  frequency  spatial  correlation  function,  after  which  we 
find  the  associated  frequency  wave  number  function  Later,  we  specify  the  frequency  wave 
number  function  directly.  We  will  find  that  this  is  an  easier  and  much  more  intuitive  approach 

At  a temporal  frequency  <o,  isotropic  noise  in  three  dimensions  is  modeled  as  a super- 
position of  infinitesimal  plane  wavt  processes,  all  radiating  towards  a common  point  [8]  These 
waves  may  be  considered  to  be  geneiated  on  the  surface  of  a sphere  whose  radius  is  large  com- 
pared to  any  geometries  or  wavelengths  of  interest  Using  the  integrated  transform  representa- 
tion we  have 


dY(w0,z) 


sin(d)  dd> 
4ir 


dY(co0  0,0)e 


-jkoar(0,Wz 


dY(co  z)ejwt,- 


and  kQ  = 2tr/X  with  ar(0,0)  as  a unit  vector  in  the  radial  direction.  Consequently. 

forms  a propagation  vector  k(d,</>)  at  a temporal  frequency  co0  radiating  towards  the  center 

of  the  sphere,  or 


iLi.: 


2rr  “o 

k (0,0)  = -kpa r(0,0)  = - y Sr(fl.«)  = - — 2r(0,0) 
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We  assume  that  disjoint  regions  of  the  sphere  radiate  uncorrelated  components  so  that 


E[dY(w0  0K0,)dY*(coo  Q242)]  =So(wo.0i>01) 


o(O)-O2)uo(01-^2)\  /dco\ 
(sin(9j)/4jr)  /Y1’) 


(2  29) 


(The  impulse  terms  should  be  interpreted  formally  and  should  operate  simultaneously.  The 
factor  of  (4ir/sin  8 j ) needs  to  be  introduced  because  of  the  use  of  spherical  coordinates ) The 
temporal  frequency  spatial  correlation  function  is  given  with  some  abuse  of  notation  by 


Sy(w0:Az)  = E(dY(w0:z)dY*(w0  z-Az)) 


(2.30) 
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When  we  use  Eq  2 27a,  we  have 


.r_%v. 
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S,,(cj0  Az) 


ft  sin(tf  j)  P1  fix  sm(0->) 

il,ii  —“•IH  ^ 


. e‘jko2r(0 1 .^l)‘£+  jkoar(02.02)’<z-Az) 


E[dY<£00.(? , ,0, )dY*(co,  02,02)j 


(231) 


Equation  2 29  implies  that  disjoint  (0,0)  is  uncorrelated  and  yields  the  desired  result 


Sy(w0  Az) 


•'o  J0 


s|n(0)  -jk  a (0,0)-Az 

d0So(«o  0,0)e  ° r 


(2  32) 


The  frequency  wave  number  function  follows  from  the  Fourier  transform  relationship  It  is 
useful  to  define  the  wave-number  k in  both  spherical  and  Cartesian  coordinates 


k = kxax  + kyay  + kzaz  = kfaf  (0k,0k) 
where  ar  is  a unit  radial  vector  in  the  same  direction  as  k We  have 

Py(cu0  k)  = J J J dAz e>h.'AlJ  c!0  * d0  SQ(wo.0,0)e 

=f«r 

Ja  Jo 


(2  33) 


sin(O) 

4ir 


d0  S0(c oo:0,0) 


j(  krar(0  k ,0k)-kQar(0 ,0))  • Az 
(2  34) 
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The  evaluation  of  the  last  integral  leads  to  an  impulse,  or  delta  function,  in  wave  number  space, 
i.e , 


„ „oo 

f/J 


dAzeJ(k*(<,k^k)-ko2r(e.«)-^  = (2ir)3 


uo(kr-ko)uo(0k^)uo(^k^) 


k0"  sin  (0) 


(2  35) 


When  we  substitute  Eq  2 35  into  Eq.  2 34,  we  obtain 


Py«ook)  = (2i03  So«oo:0k,«k) 


u0(kr-ko> 

4ir  k 2 
o 


(2  36) 


We  have  the  intuitive  interpretation  that  the  resulting  frequency  wave  number  function  has  the 
same  distribution  as  the  spectra  of  the  plane  waves  at  the  various  locations  on  the  sphere.  The 
delta  function  arises  because  these  waves  are  modeled  as  pure  plane  waves  analogous  to  pure  tones 
temporally.  By  starting  with  this  specification  of  the  noise  field  we  can  model  a large  number  of 
ambient  fields  that  may  be  encountered. 


2.2.1. 2 Series  Expansions  in  Terms  of  Spherical  Harmonics  for  Space/Time  Covariances  for 
Plane  Wave  Processes.  At  this  point  we  observe  that  once  S0(co0  0,0)  is  specified,  Py(co0k) 
follows  directly.  One  can  describe  the  distribution  of  plane  wave  power  in  the  signal  field  quite 
intuitively  in  terms  of  either  function.  We  now  demonstrate  a series  representation  for 
S0(w0  0-6).  which  enables  us  to  find  the  temporal  frequency  spatial  correlation  function 
Sy(w  Az'  quite  conveniently  Since  many  of  the  signal  processing  techniques  involve  this  func- 
tion, these  results,  coupled  with  a previous  analysis  of  purely  directional  signals,  form  a 
reasonably  complete  method  of  analyzing  ambient  signal  fields  with  a propagating  stmcture. 
The  basic  techniques  that  we  employ  draw  upon  some  results  in  spherical  expansions  which 
were  first  used  in  analyzing  similar  problems  in  electromagnetics  Our  principle  reference  is 
Stratton[9]. 
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The  expression  for  the  temporal  spectrum-spatial  covariance  function  when  represented 
m spherical  coordinates  is 

Sy(co0:z)  = f dd  d<p  S0(w0  0,0) 

J0  J0 

.e-jkorzfcos<0)cos(0z)+sin(0)sin(0z)  (sin(0)sin(0^)fcos(0)cosHz))|  (2.37a) 


where 


rz=\l\  (2  37b) 

Using  results  for  spherical  expansions,  we  can  expand  So(too.0,0)  in  a series  of  the  form  [9.  Eqs 
(399H420)] . 


°°^  n 

sq(coo.0,0)  = ? Y Snm(wo)Pl"ll(cos0)e)m^,  (2.38a) 

n*0  m— n 


where  the  coefficient  5nm(w0)  is  given  by 

2n+l  (n-imi)-'  [*  [*"  |ml 

bnm(wo)=  "jp  jJJJJ ~ T / d£?  / sm0d0  So(<o.0,0)  P n (cos0)e'jmi^  (2  38b) 

Jo  Jo 

and  the  function  * (cos 0)  is  a Legendre  function  of  the  argument  cos  0 These  functions  have 
been  tabulated  extensively,  alternatively  one  can  generate  them  with  the  Rodriguez’s  formula  (10) 


P™(x)  = (x2-l)m/2 


dmPn°(x)  (x2  |)m/2  (jm+n(x2.]jn 
dxm  2nn!  dxm+n 


(2  39) 
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Substituting  the  expansion  of  Eq  2 38a  into  the  expression  for  Sy(co0  z)  we  have,  after 
exchanging  the  integration  and  summations. 


°°  m=n 

Sy(U>0  Z)  = L E 

n=0  m=-n 


plml 

n 


(cos0)e 


-jkorzcos(0)cos(0z) 


— ^ d0  /m^Vzsin^«^)c°s(Wz)) 


(2  40) 


We  now  perform  the  integration  with  respect  to  0.  We  change  integration  variables  first  We  define 


— 0-^z  + 


n_ 

2 


(2.41) 


Substituting  this  change  of  variables  yields 


eJ(m0-korzsin(0)  sin(0z)cos  (0-0z)) 


_1_ 

2tt 


ejm(0z-ir/2  '-Vz  sm<®)  Sln<  'z> sin^ ’>1 


= (.j)mejm0z  Jm[korzs,n(0)sln(02)j 


(2  42) 


where  we  have  used  the  tabulated  integral!  10) 


(2  43) 


1 [*  if" 

Jnl(x)  - — I cos(x  sin(0)-m$)d0  = — I cos (x  sni(0)-m0)d9> 

” -V,  •/„ 


cos(x  sm(0)-m0)d0 
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.j|ni0-x  sm(0)] 


We  substitute  this  result  into  Eq  2 40,  the  expression  for  Sy(co0,z),  and  obtain 


w II 

sy<wo-£)  = X]  X]  S»»'(w«)(i)"ll'J' 


• ■jl  dfl  sin(O)  F (cosO)  Jm(kQrz)  sin(0)  sm(0z)e 


To  evaluate  this  integral  we  use  the  following  addition  formulae  for  spherical  harmonics  from 
Stratton  (9,  p 411  (Eq  2 69a)]  - 


,(k0R)Pm  (cosa)  = -j-  I dO  P^(cosO)  Jm[kR  sin(a)  sintOHel^^005^1-05^^, 
*n 


where  jn(x)  is  a spherical  Bessel  function,  which  has  also  been  tabulated  extensively.  We  have 
Rayleigh’s  formula!  1 1 ] 
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Making  the  identification  of  the  respective  variables  and  changing  the  sign  of  kQ,  we  have 
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f sin(fl)  P™  (cosO)  Jn[k0rz  s.n(0) 


Jml 


= (-J)n  Jn(korz>  P (COS0Z) 


In  summary,  we  finally  obtain 


°°  m=n 


n=0  m— n 


(2  47) 


Sy(wc  z)  = Z E W"o)(-j)"+m  Jn(korZ)  P,^11fcos<02)|eim0z  , (2  48) 
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where  we  have  from  Eqs  2 28a  and  2 28b 
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S0<“o°  ■*>  = £]  £ Snm(coo)  P^'(cos0)eJ^ 


n=0  m=-n 


and 


W")' 


2n+l 


/ (n-m)! 


47r  yn+m)! 


)f-r 

'4(>  40 


(2  49a) 


d0  sin  0 Sn(w  0,0)  p'™'  (cosOfe*11^ 


(2.49b) 
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The  expressions  for  finding  the  Legendre  functions  P (cos0)  and  the  spherical  Bessel  functions 

n 

jn(  |k0llzl ) are  given  by,  respectively,  Eqs.  2.29  and  2 46  We  point  out  that  both  of  these 
expressions  have  been  tabulated  We  have  confined  our  attention  to  plane  wave  processes 
that  originated  on  a sphere  in  wave  number  space,  i e , |k|  = 2ir/X.  Using  the  references  cited,  it 
is  straightforward  to  develop  an  appropriate  theory  when  we  relax  this  restriction. 
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Example  2a  - Isotropic  Noise 


The  simplest  noise  field  in  terms  of  this  representation  arises  when 


®nm^w^  ~ l 

/ 0,  otherwise 


We  have  from  Eq.  2.49 

So(<o:0,0)  = SQ(<o) , 0 < 0 < it,  O<0<2ir  (251) 

This  is  a commonly  used  model  for  the  ambient  noise  present  in  the  ocean  The  temporal  fre- 
quency spatial  correlation  function  follows  directly  from  Eqs.  2.46  and  2 48 


sy(re0  z)  =S0(w0)j0(k0rz)P  [cos(0z)[  = S0(<o)  sinc(kQrz)  (2.52) 


This  result  has  been  derived  in  many  places,  but  in  a much  less  general  context  [8.1 2] 

Example  2b  - Surface  ( or  Bottom  Noise) 

With  this  noise  field  we  assume  that  there  is  a much  stronger  intensity  for  the  noise  field 
in  the  vicinity  of  6 = 0 and  that  there  is  an  azimuta!  symmetry  To  do  these  we  choose 

S0( <o.0,0)  = S0(<o)  [ 1 + ocos(0)I  (2.53) 

such  that  the  noise  level  as  a function  of  0 appears  as  in  Figures  2-4  and  2-5.  This  type  of  noise 
can  be  used  to  describe  a high  intensity  noise  field  which  is  present  from  the  surface  (or 
bottom).  We  have 

!0 , m 3*=  0,  or  n > 2 

1,  n = a m=0  (2  54) 

a,,  n = 1 ,•  m=0 
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Consequently,  the  temporal  frequency  spatial  correlation  function  becomes 


Sy(to:z)  = S0M  J0(korz)  Pq0)  [cos(0z)J  ,'aj!  (k0rz)P(i°)  (cos(0z)] 


- S0(oj0)  ) slnc(k0rz)  + jo  - — (smc(k0rz)-cos(k0iz)]  cos(0z)[ 


Example  2c  - Layer  Noise 

With  this  model  of  a noise  field  we  assume  that  there  is  a much  stronger  intensity  in  the 
Vicinity  of  0 = 90°,  and  that  there  is  azimuthal  symmetry  Specifically,  we  choose 


So(<o.0,«)  = So(<uj  (l-<*P2(cos0))  =S0(co)  [1  -^a-^acos(20)]  (2.56 


The  intensity  as  a function  of  0 is  illustrated  in  Figures  2-6  and  2-7.  We  have 


10,  m st  0 
1,  n = 0,  m = 0 
0,  n = 1,  m = 0 
n = 2,  m = 0 
0 n > 3,  m = 0 


The  temporal  frequency  spatial  correlation  function  becomes 


Sy(u0:z)  = S0(co0)  j0(k0rz)  P ^°)(cos(0i.)l  + a j2(k0rz)P  ^(cosfO,)) 


= S0(oi0) 


| ["/  3 \ 3cos(k  rz)l  j"3cos(20z)+ll ) 

Jsinc(k0r z)+«  ,-lJ  sinc(k0rz)-  1 
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7 Relative  parts  of  real  components  for  the  covariance  function  for  noise  distribution 
in  Figure  2-6. 


It  should  be  clear  that  we  can  now  model  a large  class  of  processes  with  this  method  More 
complex  elevation,  or  0,  dependence  can  be  introduced  and  sectors  of  noise  in  azimuth,  or  0, 
can  also  be  incorporated 

The  model  is  particularly  convenient  when  there  are  relatively  few  significant  spherical 
harmonics.  Fortunately,  this  is  a common  situation  for  a typical  ambient  noise  field  in  many 
sonar  applications. 

It  is  worthwhile  to  point  out  that  the  analysis  in  terms  of  sphencal  harmonics  demonstrates 
that  many  ambient  field  models  can  lead  to  the  similar  oscillatory  behavior  lor  the  temporal 
frequency  spatial  covariance  function  This  stems  from  the  ideal  bandhmited  behavior  that  the 
finite  velocity  of  propagation  imposes  upon  the  wave  number  spectrum  The  oscillatory  behavior 
certainly  cannot  be  interpreted  as  being  particularly  unique  to  an  isotropic  noise  model 

2.2.2  Representations  for  Signal  Processes  in  Two  Dimensions 

In  the  discussion  of  section  2 2 we  considered  the  representation  of  signals  propagating  in 
three  dimensions  whose  wave  number  was  limited  in  magnitude  to  2ir/X.  In  this  section  we 
analyze  the  related  problem  of  the  representation  of  signals  on  a surface.  In  many  problems,  this 
representation  may  be  more  convenient  than  the  one  in  the  previous  section  Problems  in  this 
context  may  arise  in  several  applications.  In  many,  the  medium  of  interest  supports  surface 
waves;  for  example,  seismic  Rayleigh  and  Love  waves  or  internal  ocean  waves.  In  other  applica- 
tions, one  observes  signals  which  propagate  in  three  dimensions  but  are  observed  on  a two- 
dimensional  surface 

If  we  confine  our  attention  to  plane  waves,  there  is  little  ambiguity  since  only  two  points 
on  the  sphere  model  map  to  an  identical  point  in  the  plane.  If  we  assume  that  this  ambiguity 
can  be  resolved,  possibly  via  the  characteristics  of  the  receivers,  the  geometry  of  the  model,  or 
its  inrelevancy,  one  can  proceed  with  a two-dimensional  analysis  which  is  often  considerably 
simpler  than  the  corresponding  three-dimensional  one  [12,  1 3, 1 4]1 . We  point  out  here  that  in 
contrast  to  the  previous  section,  we  will  not  constrain  the  magnitude  of  the  wave  number,  although 
in  most  applications  it  will  have  a finite  upper  limit 


1 This  model  is  particularly  rele I ant  to  processing  seismic  data  for  here  the  signals  are  all  incident  from  beneath  the  earth 's 
surface,  as  partially  discussed  by  Berg.  Gaarder,  or  Capon 
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The  fundamental  relations  between  the  temporal  frequency-spatial  correlation  function 
and  the  wave  number  spectrum  remain  virtually  identical  with  the  exception  of  the  integration 
region  which  becomes  a plane  rather  than  spherical  surface.  We  have  for  the  frequency  wave 
number  function 


•»-/7 


Sv(ej-2)ei-'-  dz  (Cartesian  coordinates)  (2  59a) 


Sy(to.z)e' 


z)/krrz){co^k)  cos^z>  + sin^k>  s,n^z)|rzdrzd0z 


(Polar  coordinates) 


Similarly,  for  the  inverse  transform  we  have 


„(<oz)  = / / Pv(cj  k)ed-’-  — d~-  (Cartesian  coordinates)  (2  60a) 

y .L  J air)2 


-fT 

l")2J 0 ./o 


Py«o:ye-j(k*)  [COS(0k)  ”«»**<*$  Sin»z)l  kfdkr#k 
(Polar  coordinates)  (2.60b) 


We  can  proceed  with  an  analysis  parallel  to  that  which  we  did  for  the  spherical  representation. 
First,  we  demonstrate  that  the  analysis  of  a directional  signal  remains  unchanged,  and  then  we 
discuss  a series  representation  which  leads  to  a convenient  method  of  analysis. 

Example  3 - Directional  Signals 

One  ot  the  principal  differences  which  often  arises  in  a two-dimensional  analysis  is  that 
the  propagation  velocity  as  projected  upon  the  surface  need  not  be  a constant  unless  one  con- 
fines his  attention  to  surface  waves.  For  a directional  signal  we  have 


y(t.z)  = y0  It-  (<rz/c(ks)l 
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which  describes  the  signal  propagation  where  c(ks)  is  the  propagation  velocity  as  a function  of 
the  signal  wave  number,  k$.  This  leads  to  a temporal  frequency  spatial  correlation  function  of 


c / , c / x -jlco/c(k  >]  £*z  -jk  *z 

Sy(wz)  = S0(w)e  s =S0(co)e  s , 


where  we  have  the  relation 


, u 

k.  = a 

~s  - 


The  frequency  wave  number  function  is  still  impulsive,  or 


Py(w  k)  = S0(<o)u0(k-ks)  = S0(w)u0(kx-ks  )u0(ky-ks  ) 

x y 


To  consider  an  analysis  of  some  more  general  processes  using  a series  representation,  we  specify 
Py(oo.k)  in  terms  of  the  polar  coordinates  kf  and  We  expand  the  frequency  wave  number 
function  Py(c«j  k)  in  the  series-integral  form 


Py(ook)=  JT 


Pm(“  ki 


•'O 


Fy(Cki  k)e  Jlr"^d0|s,-  kr  = lk|  = const 


This  is  simply  a Fourier  series  decomposition  at  a specific  kr 


j 

i 


» 

i 
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To  proceed  further,  we  have  the  following  sequence  of  Founer-Bessel  or  Hankel  transforms 
For  any  function  f(p),  p > 0 with  a bounded  first  moment,  the  nth  order  (n  > - I /2)  Hankel  trans- 
form is  given  by 


Fn(X)  =/  f(p)Jn(pX)pdp, 

Jo 


while  the  inverse  transform  relationship  is 


f(p)  = / Fn(X)Jn(pX)XdX 

Jo 


(2.67a) 


(2.67b) 


(We  note  that  one  can  generally  use  Hankel  functions  rather  than  Bessel  functions,  and,  in 
particular,  when  n = 1 /2  one  has  the  common  Fourier  transform  pair  [15,  Vol  2,  p 73) ) 

We  define  the  mth  order  Founer-Bessel  transform  pair  for  pm(ca  kr)  to  be  fm(w:X)  such 
that  we  have 


fm<«:X)H 


pm^w  M'W^'r^r^r 


(2.68a) 


•-< 


and 


Pm(“'kr)  = / fm(“  X)Jm(Xkr)XdX 

J 0 


(2  68b) 


(One  could  choose  any  order  for  the  Bessel  function  or  transform  order,  however,  it  will  be 
convenient  to  choose  it  to  be  m as  indicated.)  The  net  result  of  our  decomposition  is  that  we 
can  express  Py(co:k)  in  the  series  integral  form 

£ f fm(^-X)Jm(Xkr)XdXeJm^ 
m=-<»  J o 


Py(w:k)  = 


(2.69) 


T3 


(• 


This  is  simply  a Founer-Bessel  representation  for  a cylindrical  coordinate  system,  and  it  has 
been  used  extensively  in  the  study  of  electromagnetic  fields[9] 

The  temporal  frequency  spatial  correlation  function  follows,  using  this  representation  We 
have  in  polar  coordinates  from  Eq  2 60b 


= __L  f dkr  C7'  d,k  kr  p (w  k)e'i(  Vz)cos(  W (2.70) 


Sy(W  Z) 


Substituting  Eq  2.65  yields 


Sy(u.z)=— l—  r dkrf  d(5kkr  ^ Pm(»  b.r)e'!-^k*1 
(2tt)  Jq  Jq  m=-°° 


cos(0k-07)] 


(j 7)/"  dkrkr  X]  pnito  kr)eJmW'!r/2)jn,<krrz) 
k ' J0  m— “ 


jm(0  -tr/2) 


V*  r / s jm^z' 

= Z-<  wwz)(  - 


In  summary,  we  have  the  relationships 


Py(<o:k)  = pm(w:kr)eJI 


fSl 


Sy(to  z)  = £ f 


pm(w  M4rv-'^m(w  V 


form  an  mth  order,  Fourier— Bessel,  or  Hankel  transform  pair 

One  can  use  the  numerous  transform  pair  relations  that  have  been  tabulated  by  Erdelyi 
The  essential  point  here  is  that  a large  amount  of  literature  and  tabulated  integrals  are 
available  for  these  transforms  [10] 

For  convenience,  we  point  out  a possible,  but  certainly  not  exclusive,  approach  to  the 
analysis  of  this  class  of  fields  Let  us  assume  that  the  magnitude  of  the  wave  number  is  bounded 
such  that  kr  < kQ.  Just  as  we  did  earlier,  wc  can  expand  either  pm(w  kr)  or  kfpm(w  kr)  in  an 
orthogonal  senes  of  polynominals  in  the  variable  kr/kQ  In  the  following,  we  consider  a 
Tchebychev  expansion  of  krpm(a>  kr).  We  have[  10] 


where  Tr  (x)  is  an  nth  order  Tchebychev  polynomial,  and 


Cnm  = ^ / krPm«-kr>  Tn  {j) 

* /q  Ko  Ko 

2 f°  /^r\ 

com  = 7 / krPm<w  M n=0 

" Jn  0 


We  now  have  the  tabulated  integral  [15,  Vol  2,  p 42(1)) 


f'  Tn(p)  _ ir  (\\  /A\ 

J TTTlTI  Jm<Xp)dP  = T J(m+2)/2\7;J(m-2)/2b  / 
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Substituting  into  the  transform  integral  and  using  Eq  2 74,  we  obtain 


“Ko 

fm(^:z)  = — 2_, 


cnm  -*(m+n)/2' 


/korz\  /korz\ 

\ 2 / («)/2\  2 / 


such  that 


■HsXrfe  § c- 


z\  jm(0  -tr/2) 


Since  the  Tchebychev  polynominals  are  orthogonal  over  (-1,1)  we  can  obtain  some  degree  of 
freedom  since  we  are  interested  only  in  the  region  (0,1)  Consequently,  by  choosing  pm(cj  kf) 
appropriately  for  kr  < 0 we  can  minimize  the  number  of  terms  needed  in  a finite  term 
approximation 

Alternatively,  one  could  approach  the  transform  using  Eq  2 74  by  expai  ding  Pm(co  kf) 
and  using  the  following  recurrence  relationship  in  the  transform  of  Eq  2.72a. 


xTn(x)  = 


Tn+l(x) + Tn-((x) 


Again,  the  important  point  is  not  so  much  this  particular  method  of  expansion,  but  that  there  is 
a wealth  of  results  available  using  tabulated  special  functions  which  can  be  tailored  for  an 
individual  application 

Example  4 - Circle  Noise 

A particular  noise  model  often  used  follows  when  we  choose 


S0(oj)u0(kr-k0) 

Py(<ui)  = (2rr)2 — 

y 2rrk0 
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i.e.,  this  is  a ring  of  plane  waves  radiating  towards  the  center.  We  have 


Pm(eo;kr)  = 


! So(<0)uo(kr-ko) 
i 2i r m=0 

V If 

Ko 


When  we  use  the  zero  order  Bessel  transform,  as  indicated  by  Eq  2 68a,  we  have 


S0(co)J0(k0rz),  m=0 


so  that  the  temporal  frequency  spatial  correlation  function  becomes 


Sy(oj-z)  = S0(oj)J0(k0rz) 


Example  5a  - Two-dimensional  representation  for  isotropic  noise 
Let  us  assume  a frequency  wave  number  function  of  the  form 


Pv(w:Js)  = r 1 ■ 


Due  to  the  angular  symmetry,  the  function  pm(w.kr)  is  nonzero  only  for  in=0  We  have 


Consequently,  we  have 


Sy(co  z) 


_ So(ciJ>  2,.  , . (k°Tl\ 

= 2nk0  rkoJ'/2\T/ 

7 r / 4 \ /k orz\  /korz\ 

= So(c°)  2 l^kT/Sin  \T~lCOS  \~/  So(w)sinc(korz> 


(2  86) 


We  observe  that  we  are  led  to  the  same  temporal  frequency  spatial  correlation  function  as  with 
-.pherically  symmetric  noise  (see  Eq  2 52)  In  a subsequent  discussion,  we  explore  the  reason  tor  this 
Example  5b  - Two-dimensional  representation  for  noise  with  a high  concentration  of  low 
wave  number  components 

Let  us  assume  that  the  frequency  wave  number  function  Py(co  k)  has  th»’  form 


Py(W 


1 /kr\  r /krV 

*?\ y [-y 


1/2  -1 


(2  87) 


This  wave  number  function,  shown  in  Figure  2-9,  corresponds  to  a high  intensity  source  near  kr  = 0 
which  may  be  due  to  a strong  component  normal  to  the  surface  Pursuing  the  same  analysis  we  find 


krPm(w:kr>  = S0(w)  — — 

ir-k„ 


I" 


ft)'  "'2 


(2.88a) 


SoC«) 


» kn 


n=m=0 


(2  88b) 


ns^O  or  ms*® 


Therefore,  we  have 


Sy(<o 


:£)  = So  («)Jof?) 


(2  89) 
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Figure  2-9.  Two-dimensional  noise  with  a high  concentration  of  low 
wave  number  components 


We  could  pursue  examples  at  length  using  this  procedure,  the  generality  of  the  approach, 
however,  should  now  be  apparent 


2.2.3  Representations  for  Three-Dimensional  Plane  Wave  Noise  Projected  on  a Two-Dimensional 


At  this  point  we  investigate  the  relationship  between  the  three  dimensional  representation 
discussed  earlier  and  the  two-dimensional  one  just  established,  and  which  is  appropriate  when 
three-dimensional  noise  is  observed  on  a two-dimensional  surface.  We  do  this  by  considering 
what  happens  when  we  confine  our  attention  to  a plane  m space  or,  in  particular,  whan  zz=0 
We  define 


zs  = 2 = zx-x  + zySy 


ks  = k = kxax  + kygy 
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In  terms  of  a three-dimensional  wave  number  function 


Sy(W  Zs)  = Sy(W 


1 -*!D 


Py(w  k)e"J<kxZx  + k>'Zy)  dkxdkydkz 


_ OO  r OC 

-Iff 


Py(cj  k)  e J 2 ' dkxdky/(27r)2  (2  91) 


Therefore,  the  two-dimensional  wave  number  function  is  given  by 


dk 

P2y(<-  Ks)  = / P3y(“  M) 


We  have  the  two-dimensional  Fourier  transform  pair 


Sy(W.2g)- 


j p3y("*)fr 


P3y(<0  k)  = P2y(w  ks) 


Let  us  now  examine  what  happens  when  we  project  our  plane  wave  model  on  a two- 
dimensional  surface  An  easy  way  to  do  this,  which  eliminates  many  of  the  issues  regarding  the 
delta  functions  in  several  dimensions,  is  to  establish  the  Fourier  transform  of  Sy(w  zs)  We 
have  from  Eq  2 37a  with  0z=ir/2 


mi 


ir  2ir 


r r sln(Wk)  „ „ , -jk  sin(Ok)r  cos(0z-0k) 

Syfw  ZS1  -j  dOk  J d«k-^-  So(co-Ok,0k)e  ^ -- 


0 0 
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We  first  perform  the  integration  with  respect  to  0k  over  the  regions  [O.rr/:]  and  [irl2.it] 
(This  is  necessary  since  two  points  on  the  surface  of  the  sphere  for  the  plane  wave  model 
project  to  the  same  point  on  the  two-dimensional  surface  ) We  now  change  variables  in 
each  region  Setting 


sm(0k)=  — 0 <0k  < tr/2 


(2  95a) 


sm(7r-0k)=  — , ir/2  <0k  <7r 


(2.95b) 


we  obtain 


fe)  = , 

4rrk,  2 


si 


#kkr  'Jkrrz  c°s(0z-^k) 


Kr  \2 

"■fc 


1/2 


{ S°[“  0J  +so[w-,r-sin'1(^)  . h]  } 


(2  96) 


Comparing  Eq.  2 96  with  our  two-dimensional  polar  coordinate  transform  (see  Eq  2 70), 
we  have 


P(w  ks)  = (2tr)2 


S°[“'  Sin~'fc)  + S°  ["  ~ Sm''  » J 


4irk02 


l-(^2 


1/2 


(2  97) 


The  numerator  terms  represent  those  points  on  the  sphere  which  project  to  the  point  on  a 
surface  with  a two-dimensional  wave  number  ks.  the  denominator  term  is  the  Jacobian  of  the 
transformation  Intuitively,  the  Jacobian  plot  implies  that  noise  spread  over  a unit  surface  area 
on  the  sphere  leads  to  a more  intense  value  of  the  wave  number  function  when  it  is  in  the 
horizontal  direction  than  when  in  the  vertical  direction 


It  should  be  apparent  that  we  can  also  go  from  our  surface  wave  model  to  a three- 
dimensional  model,  with  possibly  some  ambiguity  as  to  which  hemisphere  We  also  see  that  the 
factor  1/Vl  — (kr/k0)~  introduced  in  our  Tchebychev  expansion  arises  quite  naturally  as  a 
Jacobian  in  our  transformation  Consequently,  both  the  expansions  in  spherical  harmonics  and 
in  our  two-dimensional  analysis  have  many  common  results 

We  again  point  out  that  while  the  spherical  harmonic  expansion  is  quite  natural  lor  our 
plane  wave  model,  the  Tchebychev  expansion  was  only  a suggested  possible  approach  For  some 
types  of  noise  fields,  another  expansion  may  be  much  more  concise,  and  one  should  bring  the 
special  function  literature,  which  we  have  not  discussed  at  all  extensively,  to  bear. 

We  have  discussed  the  representation  of  plane  wave  signals  in  detail,  however,  we 
emphasize  the  approaches  taken,  not  the  specifics  of  a particular  example  We  now  turn  our 
attention  to  receiver  apertures  for  observing  the  signal  field 


■ r t VV  rji  tj. 


PART  2-RESPONSE  OF  ARRAYS  TO  VARIOUS 
NOISE  FIELDS 
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3.  RECEIVING  APERTURES 


In  this  chapter  we  discuss  the  general  properties  of  receiving  apertures,  or  antennas, 


with  the  emphasis  on  characterizing  their  response  to  random  excitations  Our  intent  is  to 
understand  how  knowledge  of  the  statistical  characteristics  of  the  ambient  noise  field  can  be 
used  to  achieve  enhanced  performance.  Toward  this  goal,  classical  beamformmg  theory  is  quite 
useful  so  we  introduce  results  from  this  literature  that  are  relevant  to  our  analysis. 

We  first  consider  describing  the  aperture  response.  The  development  parallels  that  of 
Section  2 in  that  we  introduce  a temporal-spatial  domain  description,  followed  by  temporal 
frequency-spatial  domain  one,  and  finally  a frequency-wave  number  one  of  the  response.  This 
is  done  by  discussing  several  commonly  used  arrays  and  analyzing  their  responses  Next  we 
study  the  role  of  sensor  noise,  which  is  extremely  important  in  the  analysis  of  the  statistical 
characteristics  of  array  processing  The  noise  structure  is  closely  coupled  to  the  array  geometry 
and  usually  does  not  lend  itself  to  a separate  analysis  Finally,  we  consider  the  issues  separating 
continuous  arrays,  i.e.,  an  aperture,  and  an  array  of  sensors,  1 e , discrete  array  We  simply  state 
that  an  array  of  sensors  is  the  spatial  dual  to  a discrete  time,  or  sampling  problem,  with  some 
added  complexities  Since  sampling  in  the  time  domain  tends  to  obscure  the  more  fundamental 
issues,  we  choose  not  to  introduce  the  corresponding  difficulties  in  the  analysis  of  spatial 
processing,  and,  we  devote  a separate  section  to  the  study  of  the  sampling  questions  introduced 
by  a discrete  array. 


3.1  ARRAY  RESPONSE  CHARACTERIZATION 


Intuitively,  it  appears  that  the  description  of  a receiver  aperture  should  be  quite  simple. 
For  example,  if  we  have  a linear  array,  as  shown  in  Figure  3-1 , which  observes  a signal  field 
y(t,z)  over  a specified  time  duration  we  would  describe  our  received  signal  as 

(T0  < t < Tf 

r(t,C)  = y(t,8aa)  for  < (3.1) 

(n  |s  | < L/2 

We  see  that,  just  as  for  a purely  temporal  process,  we  need  to  specify  the  length  of  the  obser- 
vation, however,  we  note  that  we  also  need  to  specify  the  orientation  in  space  by  means  of  a unit 
fector  aa  which  is  tangent  to  the  array. 

As  a second  example,  we  have  a circular  array  as  illustrated  in  Figure  3-2.  The  received 
signal  could  be  described  as 


1 

m 


r(t  ,d>)  = y(t 


(R(cos0o  cos 0O  cos$-sm0o  sin^)ax 
+(R(cos0o  sin0o  costf>+sin0o  sin0)ay 
+(R(-sin0o  cos<p)az 


, 0 < 0 < 2ir  (32) 


The  formal  description  of  the  observation  process  can  be  quite  tedious  As  a result,  it  is  usually 
necessary  to  keep  the  geometry  simple  in  order  to  obtain  an  intuitive  understanding  of  optimum 
array  processing  Similar  statements  can  be  made  when  describing  the  operation  of  a discrete 
array,  where  one  must  specify  the  location  of  all  the  individual  sensor  elements  For  our 
purposes  we  denote  the  array  location  by  J2,  l e . we  consider  observation  points  for  zefi  In  the 
applications  of  interest  to  us,  the  observed  signal  is  weighted,  or  shaded,  and  filtered  at  each 
point,  and  then  collected,  or  summed,  together  to  form  an  output  signal,  usually  called  a beam 
Figure  3 3 displays  this  operation  graphically  while  Eq  3 3 expresses  it  mathematically 

/"Tf 

r0(t)  = / dr 


(■ 


dz  g,(t,r  z)  y(r  z) 


(3  3) 


This  expression  has  a parallel  in  classical  filtering  theory,  however,  it  can  be  deceptively  simple, 
so  a few  comments  are  appropriate 

First,  one  should  note  that  there  are  two  issues  in  specifying  the  beam  operation  the 
geometry  of  the  array,  or  aperture,  as  given  by  and  the  weighting  pattern  g(t,r  z)  One  can 
consider  expressing  these  in  terms  of  a function  with  finite  limits  similar  to  the  analysis  ot 
temporal  functions  Generally,  however,  it  is  useful  to  separate  the  two  The  array  geometry 
imposes  more  fundamental  constraints  while  the  shading,  or  weighting,  is  adjusted  within  these 
constraints  In  addition,  the  introduction  of  geometries  of  two  and  three  dimensions  can  lead 
to  some  very  subtle  considerations,  especially  with  regard  to  spatial  impulses,  or  singularity 
functions,  when  one  uses  a theory  dual  to  the  Fourier  transform  pair  relationship  between  the 
system  impulse  response  and  the  transfer  function 

In  this  context  it  is  worthwhile  to  introduce  a simple  example  to  illustrate  the  nature 
of  these  subtleties  Let  us  consider  a linear  array  oriented  along  the  zx  axis  We  represent  its 
response  as 


r(t)  = 


fl/2 

dr  I dfi  g(t,r.C)  y(r  £az  ) 

kn  x 


(3  4) 


If  we  want  to  use  infinite  limits  and  specify  the  array  response  in  terms  of  a single  function, 
we  have 


Kt)  = 


dz  g0Q  (t,r  z)  y(r  z) 


(3  5a) 


where 


goo<t’T  = 


g(t,r  zx)  uQ(zy)  uQ(zz),  |zxl  <L.  T0<r=STf 
0,  |zxl  > L,  or  r<T0  or  r>  Tf 


(3  5b) 


Observe  that  we  incorporate  the  dependence  along  the  Zy  and  zz  axes  via  the  use  of  the 
impulses  While  this  is  a very  straightforward  example,  similar  results  appear  in  more  complicated 
contexts,  particularly  when  one  uses  an  analysis  via  transform  methods  For  simplicity,  we  use 
an  integral  representations  as  illustrated  in  Eq.  3 4 


One  can  introduce  various  formalisms  to  incorporate  these  impulse  terms,  providing  one  is 
careful  about  his  description  of  the  aperture  and  its  response 

We  now  develop  the  concepts  of  spatial  filters  in  terms  of  their  wave  number  response 
This  type  of  approach  has  many  of  the  same  advantages  as  a frequency  domain  analysis  for 
purely  temporal  processes,  and  the  two  descriptions  comp'ement  one  another 


3.2  FREQUENCY  AND  WAVE  NUMBER  DESCRIPTIONS  OF  ARRAY  RESPONSES 


Representing  signals  and  filter  responses  in  the  frequency  domain  leads  to  convenient 
and  intuitive  methods  of  analyses  In  describing  the  operation  of  spatial  filters  on  signals,  one 
can  employ  a similar  analyses  with  comparable  benefits  In  this  section  we  introduce  the 
necessary  tools  along  lines  parallel  to  those  used  for  space  time  processes 

We  assume  that  at  each  point  z on  the  aperture  £2  the  filtering  operation  is  time 
invariant,  such  that 


g(t,T  Z)  = g(t-7.z)  = g(At  z),  At  = t - T 
The  temporal  frequency  response  at  a point  on  the  aperture  is  given  by 


(3  6) 


G(co 


:i)=[ 


g(At  z)e’^w^t  d(At) 


(3  7) 


fFor  discrete  arrays,  G(co  z ) is  the  transfer  function  of  the  i**1  sensor  element ) 


In  most  of  the  analyses  that  we  consider,  the  spatial  operation  is  of  fundamental 
importance  to  us,  we  generally  assume  th'at  the  temporal  frequency  w is  fixed  in  that  we  are 
concerned  with  a narrow  frequency  band  If  the  signals  involved  are  narrowband,  this 
analysis  suffices,  if  they  are  broadband,  one  needs  to  integrate  the  analysis  over  the  frequency 
band  of  interest  The  most  interesting  and  useful  function  is  the  frequency  wave  number 
description  which  characterizes  the  response  of  the  spatial  filter  to  a plane  wave  with  wave 
number  k and  temporal  frequency  co  We  define  this  function  to  be 


g(<*> 


k)  =f°j g(r  z)e‘J(WT‘-*£)  dTdz-  f 


G(co  z)e^— dz 


(3  8) 
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A beam  pattern  can  be  obtained  from  the  wave  number  response  by  fixing  the  magnitude  of  k, 
usually  at  co/c,  and  evaluating  as  a function  of  the  elevation  angle  0 and  the  azimuth  0 

Since  this  function  is  of  particular  importance  in  our  analyses,  we  introduce  several 
examples  Tne  simplest  wave  number  response  is  for  a linear  array  with  uniform  weighting 

l 

of i-  m magnitude  and  phasing  of  e “ 1 , or 


(U2  - 

g(w  k ) = I j- 
-L/2  L 


-J^raa6 


;,-'-a!1d«  = s1„c[(k-kT)-aai]. 


(3  9) 


where  kj  is  the  wave  number  associated  with  a particular  target  direction  This  wave  number 
response  is  unity  whenever  k has  the  same  projection  in  the  aa  direction  as  k-p  The  width  of 
the  main  lobe  along  this  direction  is  1/L  between  null  points 

If  one  introduces  a triangular  shading  with  the  same  phasing 


g (r  i)  = 


0, 


I Cl  <L/2 
l«l  > L/2 


(3  10) 


gr(co.k)  = 


e-*TV 


e^V 


(3  11) 


we  obtain  the  wave  number  response 


9(0)  k)  = smc2  Jck-kj)  * aa  i-J 


(3  12) 


This  beam  pattern  has  lower  sidelobes  but  twice  the  mam  lobe  width  of  the  uniformly 
weighted  aperture 

If  one  wants  to  take  the  inverse  transform  of  these  functions  to  produce  the 
weighting  pattern,  one  should  be  careful  because  of  the  aforementioned  impulse  terms 
For  example,  assume  k-p  = §y  and  aa  = ax  We  then  have 


g(w.k)  = smc 


a 


(3.13) 
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The  inverse  transform  is 


dk  n n 

G(co  k)e 


— u0(zy)u0(zx)u0(T),  1*1  < L/2 

(3  14) 

0,-  |*l  > L/2 

Other  shadings  would  lead  to  a d.fferent  frequency  wave  number  response  Extending  the 
array  into  the  other  coordinates  just  leads  to  a parallel  analysis  For  example,  a disc  of 

radius  ra  situated  in  the  x,y  plane  with  a radial  shading  f(r)  and  a phasing  e - generates 
a beam  pattern  given  by 


I S f[(z2  + z2 


« zx  + zy  = Ra 


f[(z2  + z2))e'JkT'ZeA‘Zd 


=f  \lrf  rd0f(r)eJ<-'kT)'r(COS(^)-x  + s,n^V 

'a  Jo 


jr[(k  -kT  ) cos (<f>)  + (k  -kT  ) sm(0)) 
rd#f(r)e  x y 


fRa  f2*  jr[(kx-kT  ) c< 

I dr  I rd0f(r)e  x 

Jo  n> 

/"R  a 

= 2<r  J f(-)  J0[r|k-kT|  xy)r  dr 


ESS 


lk-klrxy  = ((kx-kTx)2  + (ky-kTy)2] 1/2  (3  16) 

Consequently,  g(co  k)  is  a zero  order  Hankel,  or  Fourier-Bessel  transform  of  f(r)  with  respect 

to  the  term  |k-k|T  In  general,  for  circular  geometries  the  Fourier-Bessel  transforms  assume 
*xy 

the  role  of  classical  Fourier  methods  for  linear  arrays  In  the  special  case  of  a ring  array 
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the  pattern  leads  to 


g^:k)  = J0  (|k-kT|  xyRa) 


If  the  array  forms  a disc,  then 


f(r)  = 


r < R, 


t 


0, 


l(|k-kTlxyRa> 
l-'-flxy  R 


r > R, 


(3  I 


(3  1 


(This  is  the  familiar  Airy  disc  response  which  is  quite  useful  in  optical  signal  processing ) 
Ring  and  disc  arrays  can  be  usea  in  two  contexts  typically  These  correspond  to 


l--T|xy  = <kx2  + ky2)  1/2 
when  the  target  is  normal  to  the  array  front,  or 


(3  1 9a) 


1/2 


(3  1 


when  the  target  propagation  direction  is  parallel  to  the  array  surface 

We  now  investigate  how  we  can  describe  the  airay  response  in  terms  of  g(co  k). 
Consider  an  arbitrary  signal  field  y(t  z).  We  can  represent  this  function  in  terms  of  a fre- 
quency wave  number  transform 


y(t-z) = 


j(wt-k'z) 

Y(oo:k)e 


dw 

2tt 


where 


Y(w  k)  = 


y(t  z)e'J^wt"-’-/  dt  dz 


(3  2 


(3  2 


•ym1 


The  response  of  the  spatial  processor  becomes 


j(WT-k-Z)  dco  dk  \ 

Y(tok)e  — lg(l-rz)drdz 


or  in  the  frequency  domain 


R0(to)  = J • J g(u  k)  Y(to  k) 


We  see  that  the  aperture  collimates  the  signal  according  to  its  wave  number  with  a weighting 
g(oj  k)  Equivalently,  we  may  consider  that  when  we  have  a pure  plane  wave  with  fre- 
quency and  wave  number  fcp;  1 e , 

N+l 

Y(<o  k)  = u0(cu-<uT)  u0(k-kT)  (2ir)  (3.23) 

we  are  observing  it  through  a window ;<pw  k)  such  that  we  have 

R0(w)  = u0(cj-'.  ) g(wTkT)(2ir)  (3.24) 

It  we  wani  to  select  some  particular  region  of  wave  number  space,  as  we  do  when  detecting 
plane  wave  signals,  <j(oj  k)  should  be  as  narrow  as  possible  However,  just  as  in  temporal 
filtering,  this  introduces  an  attendant  sidelobe  problem,  and  much  of  optimum  array  theory 
essentially  involves  determining  the  best  trade-off  of  these  two  issues  according  to  a 
statistical  measure. 

We  have  not  introduced  the  most  general  type  of  spatial  processing  In  much  of 
optical  theory  lenses  are  considered  as  wave  number  filters  that  generate  a reradiated  field 
of  the  form 

^output  - 9(co  k)  ^mput  (3  25) 

which  is  completely  parallel  to  temporal  filtering.  With  the  exception  of  some  application 
to  multiple  beam  outputs,  we  do  not  need  to  use  this  more  general  formulation  for  our 
applications  On  closer  examination,  it  does,  however,  suggest  some  interesting  possibilities 
for  implementing  our  processors 
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3.3  FILTERING  OF  RANDOM  PROCESSES 


The  next  issue  of  concern  is  describing  how  random  processes  propagate  through 
these  apertures  Here,  the  frequency  wave  number  representation  of  the  ambient  noise  field 
is  particularly  relevant  We  confine  the  discussion  to  those  random  fields  which  lend  them- 
selves to  this  descnption  Using  Eqs  2 1 2 and  3 2 1 , we  have 


>•////■ 

<«)-  Up 


dY(to  k)e 


j(cot-k-z) 


d R0(o>)  = I I l g(u  k)  dY(u>  k) 


Our  harmonic  analysis  led  to  the  result  that  disjoint  regions  of  the  frequency  wave 
number  space  had  uncorrelated  increments  dY(u>  k)  The  output  spectrum  is,  therefore, 
given  by1 


g(«:k)  " P(cu:k) 


Tlus  particular  formula  is  extremely  important  m our  subsequent  analyses  It  is  parallel  to 
the  input-output  relation  for  temporal  spectra 

^output  = IH(c«»)|~  S,npU(  (cu)  (3  29) 

We  integrate  over  all  k space  because  our  aperture  acts  as  a weighted  collimator  of 
the  ambient  plane  waves  If  we  want  to  examine  a particular  region  of  frequency  wave 
number  space,  e g , to  make  an  estimate  of  Py(w  k),  the  above  formulae  implies  that  we 
look  at  this  region  thiough  an  aperture  weighting  of  |g(w  k)|-  Ideally,  we  would  like  to 
make  this  function  as  impulsive  as  possible  The  finite  aperture  limits  our  ability  to  ao  this 
We  are  again  led  to  a trade  off  b<  tween  sidelobe  level— which  causes  the  other  regions  of  the 
space  to  interfere-and  the  beamwidth  which  compromises  our  resolution 


1 The  vector  output,  or  multibeam,  generalization  of  this  is  gnen  by 

/*  f 

k)  P(u>  k)  9 (ut  k; rr 

o j - (2*r 
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One  can  consider  | gto>  k)|-  to  be  the  power  transform  function  for  the  aperture 
It  would  be  convenient  to  define  the  aperture  autocorrelation  function  as  the  inverse 
Fourier  transform,  or 


rG(<°  2)  = 


| g(co:k)|  ^ e&'2dk 


(3  30) 


When  one  accounts  for  difficulties  with  spatial  impulse  functions  in  inverse  trans- 
forming g (to  k),  it  is  not  surprising  that  these  are  compounded  when  dealing  with  Igloj  k)|- 
unless  one  is  dealing  with  linear  or  planar  geometries 


3.4  SENSOR,  OR  RECEIVER,  NOISE 

Any  real,  or  physical,  receiving  aperture  cannot  measure  the  incident  signal  field 
perfectly,  as  the  observation  operation  is  inherently  noisy  Sometimes  this  noise  may  be 
insignificant  compared  with  other  system  noises,  however,  it  does  ultimately  set  limits 
on  the  performance  of  the  receiving  aperture  and  on  any  subsequent  processing  This  sensor, 
or  receiver  noise  typically  may  manifest  itself  in  several  different  ways  The  electronics  of  the 
sensor  elements  and  their  associated  preamplifiers,  microseisms  for  seismic  systems,  or  flow 
noise  past  the  hydrophone  for  underwater  acoustics  are  possible  sources 

In  our  discussion,  this  component  of  the  noise  process,  denoted  w(t  z)  is  modeled 
as  being  additive  such  that  the  ambient  signal  field  plus  the  receiver,  or  sensor,  noise  is 
recorded  at  the  sensor  output  The  noise  w(t  z)  has  temporal  and  spatial  bandwidths  that  are 
much  larger  than  any  other  processes  of  interest  The  essential  aspect  of  the  model  is  that  the 
observation  noise  is  uncorrelated  among  sensor  locations  For  a continuous  aperture,  Si,  this 
implies  that  the  space/time  corielation  function  across  the  receiving  aperture  is  given  by 

V 

Kw(t.r  z,f)  = E[w(t  z)w*(t  f)J  = N0S0(t-T)8n(z-|)  (331) 

The  use  of  the  operator  8jj(z-J)  deserves  some  comment  The  SI  subscript  on  the  indicates 
that  its  sifting  property  as  an  identity  operator  is  defined  only  across  the  extent  of  the  aperture 
At  first  glance  it  would  seem  appropriate  to  model  the  white  noise  as  being  uncorrelated  across 
all  regions  of  the  soatial  domain  Unfortunately,  such  a model  can  introduce  fundamental 
difficulties,  and  divergent  results  for  an  otherwise  realistic  model  are  often  obtained  For 
example,  this  occurs  in  the  filtering  of  two  dimensional  white  noise  with  a linear  array 
For  a discrete  array,  the  covariance  between  elements  is  given  by 

Kw(t,TiZj,  Zj)  = E(w(t.Zj)  w*(7.Zj)]  = N^ft-T^y  (3  32) 
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where  the  itk  element  is  denoted  by  its  location  z(  (Often  the  spectral  level  N0  is  denoted 
by  an  effective  operating  temperature  Teff  such  that  one  has 


No  “ kBTeff 


(3  33) 
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where  kg  is  the  Boltzmann  constant,  I 38  X ! 0~"’  watts/Hz-°K  ) 

As  pointed  out  above,  this  white  noise  component  is  present  in  virtually  all  physically 
motivated  problems  Sometimes  it  can  be  realistically  neglected  as  other  effects  dominate  the 
system  performance  From  a theoretical  aspect,  however,  one  can  demonstrate  that  many 
detection  and  estimation  problems  are  singular,  i e , they  predict  perfect  performance  if  the 
white  noise  component  is  not  present  For  example,  in  many  of  the  array  processing  problems 
that  we  shall  study,  the  gain  produced  by  the  array  may  become  artificially  large  as  the  number 
or  density,  of  the  sensors  is  increased  Essentially,  results  which  are  singular  imply  a perfect 
measurement  of  the  ambient  field  and  then  a resulting  cancellation  process  or  a very  sensitive 
situation  where  very  precise  knowledge  of  the  system  parameters  is  required  Therefore  we 
consider  a moie  detailed  study  of  the  role  of  the  receiver,  or  sensor,  noise  in  our  analysis  of 
array  and  apertures  Unfortunately,  the  issues  are  not  nearly  as  apparent  as  they  are  when 
one  deals  with  temporal  processes  While  duals  of  temporal  process  results  can  be  used 
extensively,  some  aspects  of  spatial  processing  have  no  duals  as  they  are  inherently  coupled  to 
the  array  geometry 

In  this  section,  we  present  some  issues  regarding  receiver  noise  We  discuss  the  possibility 
of  using  frequency-wave  number  concepts  and  the  problem  of  equating  discrete  and  continuous 
array  performance  Finally,  we  discuss  a conservation  property  which  sets  the  minimum  output 
level  that  the  receiver  noise  can  have 

At  first  inspection  one  would  agree  that  receiver,  or  sensor,  noise  can  be  modeled  by 
using  a flat  frequency  wave  number  spectrum  analogous  to  that  which  is  done  for  purely 
temporal  processes.  Intuitive  as  this  approach  may  be,  it  is  fundamentally  incorrect,  except  in 
some  very  special,  albeit  important,  situations  The  basic  difficulty  is  that  the  noise  is  coupled 
to  the  geometry  of  the  array  To  illustrate  this,  we  examine  two  situations 

Assume  that  one  is  using  a line  array  which  is  operating  in  a two-dimensional 
environment  The  array  is  capable  of  discriminating,  or  filtering,  wave  numbers  projected 
along  the  array.  Consequently,  if  one  uses  a two-dimensional  flat  wave-number  spectrum, 
l e , two-dimensional  white  noise,  the  noise  power  which  propagate-  through  the  filter  along 
wave  numbers  with  the  same  projection  is  infinite.  Similarly,  one  can  use  a ring  array  as 
described  by  Eq  3.15  which  has  resolution  in  all  directions  Here  the  wave  number 
response  does  not  fall  off  rapidly  enough  as  | k|  -*■ <*>,  it  behaves  as  I / 1 k|  and  the  noise  power 
propagating  through  the  array  to  the  filter  output  is  likewise  infinite,  i e , 


N0JJdk  J02(|k-kT|xyRa) 


(3  34) 


tav- 
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Thus  one  cannot  routinely  extend  one’s  concept  of  temporal  white  noise  in  representing 
sensor  noise  This  is  unfortunate  since  it  introduces  some  difficulty  in  applying  trequency 
domain  concepts  to  the  design  of  aperture  response  weightings  In  some  situations,  particularly 
linear  or  rectangular  arrays,  one  can  consider  the  noise  to  have  a spectrum  which  is  flat  with 
respect  to  the  wave  number  as  projected  along  the  array  surface  Using  this  approach,  one  can 
proceed  in  the  study  of  linear  arrays  more  or  less  in  parallel  to  a temporal  domain  analysis 
For  more  complicated  arrays,  e g , crossed  arrays  or  ring  arrays,  one  needs  to  be  quite  careful 
about  the  effects  of  this  noise 

In  our  study  of  arrays,  we  generate  an  aperture  weighting  that  produces  a beam  pattern 
which  is  directed  at  a specified  wave  number,  yet  suppresses  the  background  noise  If  this 
background  noise  is  composed  of  just  sensor,  or  receiver  noise,  the  optimum  aperture 
weighting  to  minimize  the  noise  power  in  the  beam  output  is  to  phase  the  array  towards  the 
target  wave  number  and  use  a constant  amplitude  weighting.  This  is  the  spatial  analog  to 
matched  filtering  Verifying  this  result  is  quite  direct  If  we  wish  to  direct  a beam  at  wave 
number  k-j-,  we  require 


g(0J  kT)  = 


/< 


G(co  z)e J T - dz  = I , 


V to. 


(3  35) 


with  the  noise  power  response  given  by 


Sn0(“)=  /dSl  /c 

n n 

= N°/ 


S„  (w)=  J dz,  J dz2  CXco  z,)  G*(<o  z2)N„  ia(lri2) 
n n 


dz  |G(oo'zl“ 


(3.36) 


Straightforward  application  of  the  calculus  of  variations  minimizes  Eq.  3.36  subject  to  the 
constraint  of  Eq  3 35  This  yields 


G(co  z)  = — L e‘^-T"- 

An 


(3  37a) 
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where 


f 

t = I dz.  the  “s 


A j}  = J dz.  the  “area”  of  the  array 

n 


(3  37b) 
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and 
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n0(«)  ; 

sn  («)  = — (3  37c) 

0 Aa  r. 

!f  we  alter  the  aperture  response  so  that  some  other  noise  source  can  be  combated  we  increase  ? 

the  effect  of  the  white  noise  at  the  output  of  the  beam  It  is  when  these  effects  reach 

equilibrium  that  we  have  one  of  the  fundamental  tradeoffs  which  one  must  make  in  deter- 

mining  the  optimum  aperture  response  ’■ 

i 
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3.5  DISCRETE  ARRAYS  VERSUS  CONTINUOUS  APERTURES 


We  have  chosen  to  pursue  an  analysis  which  models  the  aperture  as  a continuum  In  most 
physical  systems,  the  sensor  elements  are  discrete  which  leads  to  matrix-vector  formulation 
We  have  gone  to  a continuous  aperture  formulation  for  the  following  reason 

A discrete  array  is  essentially  a spatial  sampler  As  such,  the  sampling  aspects  of  the 
analysis,  particularly  the  imbedding  of  the  geometry  in  vector  notation,  often  tend  to  cloud 
the  more  fundamental  issues  of  the  spatial  processing  Many  arrays  in  current  or  proposed 
systems  involve  interelement  spacings  which  are  so  small  that  one  is  well  above  the  spatial 
Nyquist  sampling  frequency  Here,  a denser  sampling  of  any  coherent,  or  spatially  bandhmited, 
part  of  the  signal  field  is  redundant  In  these  situations,  the  continuous  approach  is  more 
informative  regarding  the  actual  processing,  and  the  performance  of  the  physically  discrete 
system  is  closely  approximated  by  the  continuous  aperture.  Remember  here  that,  in  contrast 
to  temporal  processes,  one  is  always  working  with  bandhmited  signals  that  arise  through  propa- 
gation in  the  medium  They  are  strictly  bandhmited  due  to  the  fixed  upper  limit  of  2t r/A  for 
the  maximum  value  of  the  wave  number,  or  spatial  frequency  When  one  pursues  a temporal 
analysis,  the  continuous  representation  is  usually  more  natural,  although  the  implementation 
may  be  done  using  a sampled  system  with  digital  filters  This  does  not  imply  that  we  can 
neglect  these  questions  It  is  simply  our  assertion  that  we  feel  many  of  the  concepts  of 
optimum  processing  are  more  transparent  and  can  be  approached  more  directly  using  a 
continuous  analysis. 

There  are  two  issues  which  concern  us  in  comparing  discrete  and  continuous  arrays 
First,  if  we  have  a beam  pattern  which  we  generate  with  a continuous  model,  we  want  to 
determine  the  interelement  spacing  necessary  to  produce  approximately  the  same  beam 
pattern  by  using  a discrete  array  Second,  we  want  to  equate  the  effects  of  sensor,  or  receiver 
noise,  for  discrete  and  continuous  arrays  so  that  we  can  compare  their  performance  in  subse- 
quent sections  To  solve  both  of  these  problems  in  general  for  an  arbitrary  array  geometry 
is  quite  difficult  Significant  understanding  can  be  obtained,  however,  by  examining 
linear  arrays 

Let  us  assume  that  we  have  a linear  array  with  «n  aperture  weighting 
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G(t  z)  for  z = 2aa,  I C|  < L/2 , 


(3  38) 


which  produces  a beam  pattern  g(oc  k)  If  this  array  is  sampled  with  a spacing  AL  and 
the  resulting  beam  pattern  analyzed,  one  obtains 


L 

2AL 


gs(w  k) . 


G„(<o)e 


-jk-a  AL  n 


(3  39a) 


2AL 


where 


Gn(u)  = G(u  naaAL)  (3  39b) 

g(c o k)  and  gs(w  k)  can  be  related  via  a direct  parallel  to  the  temporal  sampling  theorem 
This  yields 


?s<«  y = zr  9(u  (3  40) 

n= 

Consequently,  we  have  the  spectrum  repeated  at  intervals  of  (2jr/AL)aa  Two  effects  are 
significant  here  If  g(u  k)  is  significant  for  k outside  the  region  | k-aa|  < — then  dis- 
tortion is  created  via  aliasing  a classical  problem  of  temporal  filters  Thus,  we  should 
consider  the  factors  that  govern  the  spatial  bandwidth  of  the  arrays  Conventionally,  the 
minimum  beamwidth  measured  in  radians/m  of  an  array  of  length  L is  on  the  order  of 
(2ir/L),  or  l/(L/X)  radians  The  implication  for  this  minimum  bandwidth  is  that  at  least  two 
samples  spaced  at  an  interval  of  AL  < L/2  would  suffice  The  difficulty  here  is  that  this 
response  would  be  reproduced  at  intervals  of  2jrn/AL  or  2ir/L  which  could  introduce 
significant  sidelobe  issues  With  larger  beamwidths  moic  samples  are  required,  but  the  side- 
lobes  remain  To  alleviate  this,  one  usually  shifts  these  sidelobes  out  beyond  the  region  of 
propagating  signals,  l e.,  beyond  the  wave  number  region  in  which  other  sources  could  enter 
In  effect,  one  is  controlling  the  beam  pattern  across  the  region  |kl  < 2ir/\,  not  simply  across 
the  main  beam  region  Examination  of  the  proof  of  the  sampling  theorem  shows  that  this  is 
more  of  the  essence  than  just  the  simple  prevention  of  aliasing  The  sampling  required  then  is 
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which  is  a classical  result  in  array  theory 

Related  to  this  issue  is  the  theory  of  superdirective  arrays,  defined  in  the  conventional 
sense  to  be  considerably  narrower  than  that  indicated  by  classical  theory  with  either  discrete 
or  continuous  signals  This  is  done  at  the  expense  of  creating  a large  sidelobe  structure  outside 
the  region  where  propagating  signals  can  arise,  i e , for  |k|  > 2tr/X.  In  this  case  significantly 
denser  samplings  will  be  required  to  approximate  a continuous  array  and  to  prevent  aliasing 
Alternatively,  one  could  start  from  a discrete  formulation  directly,  as  is  done  in  the 
Dolph-Techebychev  ( 1 6]  theory  When  more  than  one  coordinate,  or  dimension,  is  intro- 
duced, a number  of  issues  appear  which  do  not  have  analogs  as  developed  in  the  classical 
temporal  theory. 

Tile  second  issue  that  we  wish  to  consider  is  the  effect  of  sensor,  or  receiver,  noise 
In  particular,  ideal  sampling  of  such  a noise  field  by  a point  element  cannot  be  realized  The 
question  arises,  then,  of  how  we  can  compare  the  continuous  white  noise_  level  of  Eq  3.3 1 ,. 
for  an  array,  with  finite  extent  and  discrete  white  no.se  as  expressed  by  Eq  3 32,  for  a 
finite  aperture  element  array 

The  key  to  this  analysis  is  to  consider  that  a discrete  element  is  small  enough  that 
any  coherent  part  of  the  noise  field  is  observed  undistorted  while  the  incoherent,  or  white 
part,  is  averaged  over  an  effective  element  region  Essentially,  we  have  a spatial  parallel  of  a 
finite  time  averager  Let  us  consider  an  example  relating  the  equivalence  of  a continuous 
linear  array  and  a very  densely  spaced  linear  point  array 

We  assume  that  the  l1*1  element  of  the  discrete  linear  array  is  located  at  location 
zn  = (nAL,0,0)  and  the  interelement  spacing  is  AL,  as  shown  in  Figure  3-4  By  assumption, 
this  is  very  small  compared  to  any  spatial  correlations  of  the  coherent  parts  of  the  observeu 
field  to  the  noise  levels  We  consider  that  the  output  of  the  nth  sensor  is  even  by 

zx  +AL/2 

1 f " 

rn(t)=  7T  y(tz)dzx  (3  43a) 

J zY  -AL/2 
xn 

where 

zv  = nAL  (3  43b) 

*n 


K: 


I 

1 


j 


l 

i 


Jw 


56 


CONTINUOUS  LINEAR  ARRAY 
LENGTH  L 


DISCRETE  LINEAR  ARRAY 
LENGTH  L,  SPACING  AL 


Figuic  3*4.  Continuous  and  discrete  array 


If  we  consider  y(t  z)  to  be  composed  of  a coherent  part  s(t  z)  and  spatial  white  noise  w(t  z) 


with  level  N„  , we  have 


r„(t)  - s(t,zn)  + wn(t)  . 


where  s(t,jn)  is  the  approximation  for  the  coherent  part  of  the  signal  and  the  sensor  noise  is 


1 f"" 

"n(t)=zrj 


w(t:z)  dzx 


zv  -AL/2 
*n 


The  cross  covariance  between  elements  of  the  wn(t)  is  given  by 


E[wn(t)wm(r)]  • 


— 5nm8(t-r)  = N0d6nm6(t-r) 


N = . 

°d  AL  ’ 


i.e.,  the  effective  noise  density  for  the  discrete  element  is  smoothed  over  the  sensor  extent, 
so  that  the  larger  the  sensor  face  the  smaller  is  the  noise  level.  We  finally  point  out  that  this 
noise  is  not  isotropic  as  it  is  uncorrelated  from  sensor  to  sensor  no  matter  how  close  they  are 
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In  discussing  optimum  array  proc.  sing,  we  consider  the  comparison  of  discrete  and 
continuous  models  m more  detail  This  equivalence  can  be  made  more  specific  and  the 
validity  of  some  rather  qualitative  statements  in  replacing  a discrete  array  by  a continuous 
equivalent  is  somewhat  surprising 


3.6  KARHUNEN-LOEVE  EXPANSIONS  FOR  SPACE/TIME  PROCESSES 


In  the  design  of  array  processors  which  utilize  the  statistical  properties  of  the 
ambient  noise  field,  one  often  needs  to  visualize  the  signal,  which  is  observed  over  a segment 
of  time  and  across  an  aperture  in  terms  of  an  orthonormal  expansion  This  is  particularly 
true  for  detection  theory  issues  By  using  this  expansion,  one  can  design  and  specify  the 
resulting  performance  of  the  array  processor,  which  operates  on  a random  field  by  apply- 
mg  classical  detection  and  estimation  methods  to  the  generalized  Fourier  coefficients  of 
the  expansion  of  random  variables.  This  procedure  becomes  particularly  easy  when  the 
functions  in  the  orthonormal  expansion  are  chosen  such  that  the  Fourier  coefficients 
are  uncorrelated  This  is  the  essence  of  a Karhunen-Loeve  representation  [17]  whose 
use  is  well  established  m random  process  theory.  We  discuss  briefly  some  aspects  of 
this  representation  The  only  issue  which  is  somewhat  peculiar  to  our  problem  is 
that  we  need  to  incorporate  the  spatial  aspects,  this  generalization,  however,  is  straight- 
forward. 

At  this  point,  we  are  forced  to  make  a somewhat  artificial  distinction  between  con- 
tinuous apertues  and  discrete  arrays  With  a continuous  model  we  need  an  expansion  of  the 
processes  over  the  product  space  created  by  the  observation  time  and  the  receiving  aperture, 
while  with  discrete  array  a more  natural  approach  is  to  use  an  expansion  tor  vcclo,  random 
processes  where  each  component  represents  an  individual  sensor  output.  The  issues  in  the 
development  of  a vector  Karhunen-Loeve  expansion  are  discussed  by  Van  Trees  [17]  Since  we 
pursue  an  analysis  that  treats  the  array  as  having  a continuous  aperture,  we  use  the 
former  expansion  over  the  continuous  space  created  by  the  observation  time  and  the 
receiving  aperture.  Generally,  it  is  quite  easy  to  take  the  results  and  generate  the  parallel 
discrete  array  or  vector  formulation  We  emphasize  that  the  expansion  is  done  over  a 
specific  receiving  aperture  and  observation  tune  rather  than  over  all  space 

We  state,  without  proof,  the  folio  mg  theorem  (Ref  5,  6 and  1 7)! 

Assume 

1 that  y(t,z)  is  a space/time  random  process  defined  over  a region  which  includes 
the  interval  T0  < t < Tf  and  the  spatial  domain  Jit  and 

2 that  y(t,z)  has  finite  means  square  value,  i e , y(t  z)  e or 


E[|y(t,z)l2]  <~,  Vt,ie[T0,Tf],n 

lrnte  essential  issue  in  the  proof  is  the  validity  of  Mercers  theorem  where  the  expansion  region  is  the  spaceftime  field  rather 
than  just  a one-dimensional  closed  interval 
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We  can  then  represent  y(t  z)  in  an  orthogonal  expansion  with 


N 

X 

N-*»  n=I 


a)  y(t  z)  = hm  yn«n(t  z),  te[T0,Tf),.  zeCl, 


where  y, 


{ dt  / dz*n 


*(t  z)y(t.z)dz  (Generalized  Fourier 


coefficient) 


(t  z)0m(t  z)  dtdz  = Snm  (Orthogonal  Basis 


T0  13 


functions) 


c)  E[ynym»  ] ^n^nm 


(Uncorrelated 

Coefficients) 


if,  and  only  if,  the  Karhunen-Loeve  integral  equation 


r,  [ 

\„<>n tt  z)=  J dr  /df 

T„  V. 


dfKy  (t,r  z&tn(T:[) 


(3  46a) 


(3  46b) 


(3-471 


(3  48) 


(3  49) 


is  satisfied  This  is  a very  straightforward  extension  of  the  Karhunen-Loeve  theorem  for 
processes  defined  solely  over  a temporal  domain  or  index  set. 

In  general,  the  actual  determination  of  the  eigenvalues  and  eigenfunctions  is 
difficult.  For  the  case  ot  processes  defined  over  only  a temporal  domain,  several  basic 
methods  have  been  developed  for  determining  analytical  solutions  (see  ref.  3 and 
references  therein)  Except  for  some  special  cases,  one  usually  must  resort  to  computational 
procedures  to  implement  these  methods  The  extension  of  these  methods  to  space  time 
processes  would  require  that  (a)  Py(cJ  k)  is  a rational  function  of  the  components  of  k, 

(b)  there  is  a partial  differential  equation,  describing  the  generation  of  the  process,  i e , a 
distributed  system  model,  or  (c)  there  is  some  very  special  property  of  the  kernel  which 
leads  to  an  analytic  result 


59 


* ^ **  ^ ^ Arv.  <'-*--*-1  . »u •. '%  * . •\*v«\\\v,_v_v 

ls$ 


Even  if  one  of  these  situations  does  exist,  incorporating  the  geometry  of  the  array 
poses  some  very  fundamental  problems  The  overall  situation,  then,  is  that  direct 
implementation  of  this  theorem  is  by  no  means  a trivial  task  A similar  situation  exists 
for  discrete  arrays  with  many  elements  If  there  are  a large  number  of  significant  eigenvalues 
in  the  expansion,  this  implementation  rapidly  becomes  unmanageable  with  an  array 
containing  a large  number  of  elements  Consequently,  in  our  approach  and  discussion,  we 
use  some  relatively  simple  examples  which  can  be  approached  analytically  These  highlight 
the  major  issues  involved  and  lead  to  considerable  insight  into  the  behavior  of  more  com- 
plex problems  Fortunately,  in  many  cases  of  practical  interest,  we  can  determine  approxi- 
mate solutions  by  “factoring”  the  temporal  and  spatial  dependence.  This  leads  to  an 
analysis  which  is  done  in  the  temporal-frequency  domain  Consequently,  for  a narrow 
band  situation  we  essentially  have  a complete  analysis,  whereas  for  a broadband  situation 
we  must  consider  the  issue  of  integrating  our  result  over  the  band  of  interest  For  temporal 
processes  with  stationary  kernels  and  “long"  observation  intervals,  we  have  a well-known 
result  that  the  eigenfunctions  become  sines  and  cosines,  or  complex  exponentials,  and 
that  the  eigenvalues  have  the  same  distribution  as  the  spectrum  of  the  p-ocess  This  result 
allows  us  to  “factor”  the  temporal  issues  out  of  the  problem.1  The  essential  step  in  this 
factoring  is  to  assume  that  the  solutions  to  the  integral  equation  of  Eq.  3.49  have  the  form 

1 in^o4 

*i  = — e *m<nwo  iK  (3  SO) 

>F 

where  oj0  = 2tr/T  and  we  have  replaced  the  single  index  l with  the  double  index  n,m 
(Since  both  are  countable  schemes,  this  is  permissible  ) If  we  substitute  this  into  Eq  3 49 
we  obtain 


\n<nwo  F 

yjT 


jnw0t 

e VH1) 


f Ky(t-r  z,f)  eJn“°  4>m(nw0:pdrdJ 

a \T 


jOj(t-T) 
Sy(cj:z,£)e  dw_ 
2?r 


jnwQr 

e $m(nto0:£)drdf 


* The  term  “ factor  ” is  somewhat  misleading  since  the  temporal  dependence  is  sttll  imbedded  in  the  spatial  parameters 
eg,  the  magnitude  of  the  wax  e number  in  our  plane  hsic  model 
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(3  51) 


The  term  in  the  braces  of  the  last  expression  of  Eq  3.51  approaches  an  impulse  as 
T becomes  large  We,  therefore,  have 


1 Jn“ol 

Vn“o£)-e  $m<”M 

yT  (3  52) 

i Jn<V 

The  final  result  is  that  when  dealing  with  temporally  stationary  random  fields  observed  over 
a long  time  interval,  we  can  obtain  an  approximate  solution  to  the  Karhunen-Loeve  equation, 
as  expressed  by  Eq.  3 49,  by  factoring  the  complex  exponential  temporal  dependence  and 
solving  the  equation 


l 


Sy(n<o0:a,J)  <J>m(n<u0  f)df,  — i<  t <4  J.e  S2 


\m(nu0)  ^m(nu0:z)  = I Sy(nw0  z £)  4>m(nw0  f)d£,  z e 12  (3  53a) 

Jn 

In  this  for, nulation  we  can  solve  for  the  spatial  factor  with  the  frequency  of  the  complex 
exponential  assuming  the  role  of  a parameter  In  the  case  of  linear  arrays,  Eq  3 53  is 
identical  in  structure  to  the  equation  which  needs  to  be  solved  m the  Karhunen-Loeve 
expansion  for  temporal  processes  Consequently,  the  repertoire  of  methods  that  exist  for 
solving  this  equation  can  be  brought  to  bear  here.  For  example,  isotropic  noise  with  its 
spectrum  Sv  (co)  sine  [kn  |Az|]  corresponds  to  the  situation  with  bandlimited  noise  which 

y0  u 

has  prolate  spheroidial  waveforms  for  its  solution  Unfortunately,  the  temporal  situation,  in 
which  the  spectrum  is  rational  does  not  have  direct  application  since  it  implies  a wave  number 
function  of  infinite  extent  Since  all  of  our  propagating  signals  have  a bandlimited  wave 
number  spectrum  at  a given  temporal  frequency,  such  models  can  serve  only  as  approxima- 
tions Array  geometries  other  than  linear  ones  can  be  quite  difficult  to  analyze  and  usually 
lead  to  copious  use  of  the  higher  transcendental  functions  This  is  one  of  the  reasons  for 


using  relatively  simple  array  geometries,  for  we  can  extract  the  fundamental  concepts  under- 
lying the  array  processing  operation  without  becoming  emeshed  in  a complicated  analysis 
Even  at  this,  however,  we  must  develop  some  nontrivial  spatial  results  for  which  there  is  no 
temporal  parallel 

In  most  ot  the  problems  of  interest  to  us,  the  observed  process  r(t  z)  is  compose  1 of 
two  components  one  has  a finite  mean  square  value,  the  other  is  space/time  white  noise  The 
first  of  these  has  a well  defined  Karhunen-Loeve  expansion  while  the  v.  hite  noise  does  not 


M0  S(t-T)  «n(2-_»  = N0  2)  ?).  - ^ < t.  T <y, 


for  any  complete,  orthonormal  set  Consequently,  we  can  operationally  use  the  result  that 
white  noise  is  “white”  in  any  coordinate  representation,  m that  it  has  eigenvalues  of  N0  when 
projected  against  any  complete  orthonormal  representation  Expanding  a signal  with  an 
independent  additive  white  noise  component  we  have 


Kr(t,r  zj)  = Ky(t,r  ?J)  + N0  5(t-r)  Sn(z-f) 


OO 

= XI  (Xm  + No><W‘  S>  *m<T  P>  T<7> 


z,£e  ft  (3  53c) 


so  that  the  eigenfunctions  remain  as  they  are  tor  y(t  z)  alone,  while  the  eigenvalues  become 

r y 

X =>  A + 


Henceforth,  we  suppress  the  r and  y superscripts 

Finally,  we  note  in  the  case  of  what  are  termed  “separable  kernels”  we  can  solve 
Eq  3 53a  exactly  The  situation  of  interest  to  use  here  is  that  when  the  spectral  covariance 
function  is  composed  of  directional  noise  components,  i e., 


sy(wz,r)=]T^  S,j(0J)e 
i=l  j=l 


J(k,-z-kj-f) 


the  solution  is  also  composed  of  terms  e 1 " 1 e , 


‘Wlwo  2)  = Cmk <ntJo>e  7-e  n 
k=l 


Substituting  this  yields 


Zjk. -z 
Cmi(n<o0)e 


f ^52  5]  SU(n“o)eJ"'  “~J  ® ]CCm«(n"°>e"£  ~ 


to  1=1  j=l 


N . N 


= 52  eJ-»  - ^ Vn“o)  52  Cm£(ntoo>An  sincfi  IkjHtgl, 


i=l  j=l 


smcjj(k)  = • 


[ jJS'S 
1 * ‘ 


The  solution  to  this  is  given  by 


\n(nw0) 


I-  [S(ncj0)J  (sinc^CAk)  ] 


1 


smc|k|-k2] 


[smcjjfAk)]  - 


smcn[k2-k,] 


(3  57c) 


(3  57d) 


This  is  a homogeneous  matrix  equation  which  has  N solutions  at  any  given  frequency 
(nto0)  We  will  discuss  the  solution  to  it  at  various  points  where  we  encounter  noise  fields 
with  directional  components. 

We  now  consider  the  application  of  the  results  to  our  fundamental  problem  of  mterest- 
the  design  and  performance  analysis  of  optimum  processors 


4.  OPTIMUM  ARRAY  PROCESSING  FOR  PLANE 
WAVE  SIGNALS 


In  the  previous  sections  we  have  examined  the  representation  of  space/time  random 
signals  and  the  response  of  apertures  during  observation  Until  now,  the  two  most  important 
quantities  were  P(w  k),  the  distribution  of  a homogeneous  space/time  signal  in  temporal  and 
spatial  frequencies,  and  g(co  k),  the  response  of  an  aperture  to  a plane  wave  signal  with  fre- 
quencies co  and  k In  the  process  of  doing  this  we  have  now  developed  the  concepts  necessary 
for  describing  array  methods 

In  the  design  of  an  array  processor,  the  structure  of  the  desired  signal  and  the  ambient  and 
aperture  noises  all  have  a fundamental  influence.  The  classical  problem  in  array  processing 
consists  of  finding  the  optimum  beam  pattern  within  the  constraints  of  the  array  geometry  for 
observing  and/or  detecting  a plane  wave  signal  with  a wave  number  kp  Usually  this  is  done  in 
a somewhat  intuitive  manner  Beamwidth  and  sidelobe  level  must  be  traded  off  and  nulls  are 
positioned,  all  under  the  influence  of  the  constraints  imposed  by  the  array  geometry  and  the 
noise  field  encountered.  The  methods  of  optimum  array  design  essentially  address  the  problem  of 
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giving  a rationale  to  this  procedure  in  that  the  “optimum”  processing  method  is  specified  for 
a given  array  and  noise  field.  Naturally,  we  must  make  our  definition  of  optimality  precise, 
and  the  formulation  of  and  solution  for  the  “optimum”  array  processor  is  the  topic  of  this 
chapter. 

Like  many  problems  in  which  linear  processing,  quadratic  performance  measure,  and / 
or  Gaussian  statistics  are  introduced,  this  one  is  also  quite  robust  in  that  the  same  solution 
can  be  arrived  at  via  several  different  methods.  We  have  chosen  a formulation  which  empha- 
sizes the  linearity  and  quadratic  aspects  of  the  modeling.  Various  app-  aches  have  been 
unified  for  discrete  arrays  [18] 

We  focus  upon  processing  for  signals  which  are  plane  waves  At  the  end  of  this  section 
we  comment  on  array  processing  for  spatially  spread  signals,  1 e , the  signal  is  composed  of  a 
possibly  uncountable  number  of  plane  waves  Despite  their  importance  as  a more  realistic 
model  for  many  problems  and  the  numerous  parallels  to  temporal  random  signal  detection, 
the  literature  does  not  have  many  contributions  on  spatially  spread  signal  detection.  We 
shall  subsequently  discuss  some  aspects  of  these  problems. 

One  of  the  important  mathematical  entities  introduced  here  is  the  inverse  kernel 
This  is  essentially  the  continuous  analog  to  an  inverse  matrix  and  frequently  appears  m 
detection  and  estimation  problems.  Many  of  the  operations  appropriate  for  matrices  have 
parallels  in  the  continuous  case.  This  function  is  important  in  our  analyses  in  later  chapters 
Therefore  we  emphasize  its  importance  early.  Van  Trees  offers  a more  detailed  discussion  of 
this  function  [3], 

In  this  section  we  have  a very  concise  formulation  and  solution  for  finding  the 
optimum  beam  pattern,  or  array  processor.  This  solution  can  be  invoked  to  explain  how  large 
numbers  of  array  processors  operate,  from  superdirective  arrays,  to  endfire  gain,  to  null  place- 
ment procedures,  and  to  sidelobe  level  effects.  It  provides  a unified  approach  to  the  design  of 
beam  patterns  and  array  processing.  One  of  the  goals  of  our  later  chapters  is  to  develop  enough 
important  examples  so  that  the  mathematics  become  more  transparent  and  the  underlying 
theory  can  be  used  in  an  intuitive  and  practical  manner 


4.1  OPTIMUM  WAVE  VECTOR  RESPONSE  FUNCTIONS  AND  MINIMUM 
VARIANCE  ARRAY  PROCESSES 


Our  mathematical  formulation  for  determining  the  optimum  beam  pattern  reflects 
our  intuitive  concept  of  how  a good  array  processor  should  perform.  We  want  to  observe  a 
plane  wave  signal,  which  we  call  the  target,  with  temporal  frequency  and  spatial  wave 
number  ky.  To  do  this  we  want  to  direct  a beam  in  this  direction  and  obtain  a response  which 
has  a minimum  of  interference  caused  by  ambient  and  sensor  noise  We  may  want  to 
determine  the  signal  wave  shape  itself  or  we  may  simply  want  to  measure  its  statistics 
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Mathematically,  the  operation  of  directing  a beam  is  introduced  by  requiring  that 


g0tokT|kT)=  / GQ(w  z|k()^T  dz=  1, 

•n 


Vco 


(4. 1)1 


while  minimizing  the  effects  of  the  noise  leads  to  the  requirement  that  the  output  noise 
power 


2 

a 

o 


(w|kT)  = S0(w)  = 


G0(co.z|lcT)Sn(u>  z,f)G0*(wl  |*.T)  dzdj; 


(4  2) 


t 

! 

5 


be  a minimum  Note  that  a ^ (wlky),  the  minimum  noise,  is  a function  of  ky  and  that 
Sn(tJ.z,£)  is  the  temporal  frequency  spatial  correlation  function  observed  across  the  array 
aperture.  The  constraint  of  unity  response  is  a convenience.  We  could  require  only  a finite 
response;  but  the  result  is  deterministic,  so  we  could  arbitrarily  scale  it  back  to  unity. 

If  the  wave  number  versus  frequency  for  a plane  wave  signal  is  given  by  kyfto),  then 
requiring 


?o  [w  JST<«)  liiT<<u)]  = 1,  (4  3) 

will  lead  to  an  array  processor  which  observes  this  signal  undistorted  This  approach  has  been 
used  to  formulate  the  optimum  distortionless  filter  for  discrete  arrays 

We  are  presuming  that  g(co.klky)  can  be  designed  independently  for  all  cj,  i.e  , the 
solution  at  Uj  is  not  influenced  by  the  solution  at  wy  The  total  broadband  output  noise  for 
observing  a signal  is  given  by 


A pictorial  representation  of  our  formulation  is  indicated  in  Figure  4-1 
The  solution  for  the  optimum  beam  pattern  can  be  found  by  straightforward 
application  of  the  calculus  of  variations  where  we  impose  the  constraint  using  a Lagrange 
multiplier  (3(cu).  For  optimality  we  require  that  for  all  choices  of  8G0(<o.z|ky)  that 


f We  use  the  notation 

g(ui  k\kT) 

to  represent  a beam  pattern  in  £ fora  beam  "directed"  at  ware  number  ^,ie , 
3(us  kj-jkj-/  B I 
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This  implies  that 


/ Sn(«:z,» 


jkj’i 

G0*(<J_f]  k-p)df  = -(3(o))e  ,.  ze£2,Vto 


We  first  need  to  eliminate  the  Lagrangian  multiplier  (3(o>).  This  can  be  done  by  applying  the 
constraint  of  Eq  4.1,  and  we  have 


-f5(w)=l  / G0(ioi|JsT)Sn{toz,r>G0*(coi|kT)dzdf 


which  implies  that  0(oj)  is  a real  quantity 

Substituting  this  into  Eq.  4.6  yields  the  implicit  solution  for  G(oj:z.|k-p). 


/ Sn(ca.Z, 


,nG0*(u>i|kT)dr 


I f G0(io:z'|kT)  Sn(cc:i,r)G0»(«:r|iiT)di  dI' 

JnJn 


ze  SI,  Vco 


Eq.  4.6  specifies  the  optimum  beam  pattern  as  the  solution  to  an  integral  equation 
which  appears  often  in  solving  temporal  detection  and  estimation  problems.  Some  reflection 
should  indicate  that  we  have  a problem  almost  parallel  to  that  of  detecting  a known  signal  m 
colored  noise.  We  have  a known  signal  structure,  e^-T  plus  a background  of  spatial  colored 
noise,  n(t:£).  With  the  exception  of  a normalization  constant,  the  integral  equation  for  the 
optimum  array  weighting  G(wz  |fty)  and  the  optimum  correlation  detector  are  the  same. 

The  solution  for  G(w.z  |kj)  can  be  found  in  a number  of  ways.  First,  we  observe 
that  we  really  want  to  solve 


J Sn(cou£)G0*(«:Dd.£  =ekT‘Z>  *en’Vco 


and  then  scale  the  solution  in  amplitude  to  satisfy  the  constraint 

We  want  to  introduce  the  spatial  generalization  of  the  inverse  kernel.  The  defining 
property  of  this  function  >s  specified  by 


/ Sn(wz£')  Qn( 


,(«!'£)«'  = z.fefi, 


1 e.,  the  spatial  functional  analog  of  an  inverse  matrix.  We  also  have 


/Qn(<o:z,z')  Sn(w.z',f)  dz'  = 5n(z-_£), 


which  corresponds  to  the  result  that  right  and  left,  or  post  and  preinverses,  of  a matrix  are 
identical.  One  of  the  important  properties  of  the  inverse  kernel  which  we  use  relates  to  its 
expansion  in  terms  of  its  spatial  eigenvalues  and  eigenfunctions.  Following  the  temporal 
theory  [3]  for  example,  we  can  expand  both  Sn(to  z,f)  and  Qn(co:z,f ) in  terms  of  the  eigen- 
values and  eigenfunctions  found  in  Section  3 6 in  the  discussion  of  Karhunen-Loeve  expan- 
sions Mercer’s  theorem  asserts  that  we  have  the  expansions 


Sn(<o.z,j;)  = ^ Xm(<o)4.m(w.z)tf>m»(oj  f),  z,|  e H 


oo 

Qn(w2£>=  X]  xm  (w)*m("2)*m ("i),  hiesi  (4.12b) 


where  Xm(<o)  and  4‘m(cj:z)  are  the  solutions  to  Eq.  3 53a.  In  most  cases  n(t  z)  is  composed 
of  a correlated,  or  finite  mean  square,  component  and  a white  one.  For  this  we  have 


Xm(co)  = N0  + X'  (w) 
m 


where  \y  (to)  is  the  eigenvalue  of  the  correlated  component  We  then  have 
m 


Qn(oj:zJ-)  = ^](N0  + Xy  (w))'1  $m(co:z)  4»m(wj) 


= — (S(t-r)  8nU-p  - H(u:i,ni,  z,l  e a 


I The  function  is  an  impulse,  or  generalized,  function  which  is  defined  only  across  the  array  or  aperture  ft  Const 

quently,  no  sifting  property  outside  of  this  domain  is  implied.  The  same  function  was  used  in  discussing  sensor  noise. 
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where 


“ A*  («) 

H(«o:z,f)  = y.  <J>(ajz)4>(o;  f),  z,tf  (4.12e) 

No  + ^(“)  ' 


All  of  these  results  with  inverse  kernels  are  complete  duals  to  the  temporal  situation  therefore 
we  refer  to  ref.  [3}  by  Van  Trees  fora  more  complete  discussion 

Returning  to  the  explicit  solution  for  G(w.^ljop)  we  conjugate  Eq  4.9,  employ  the 
symmetry  properties  of  Sn(co.z,£)  and  then  post  multiply  by  Qn(to \z',z)  and  integrate  with 
respect  to_f . Using  Eq.  4. 1 1 this  yields 


/-)kT-z. 

e Qn(io:z^)  dz,  | < 


Apply  the  constraint  of  Eq.  4. 1 and  we  have 


and  the  final  solution  for  G0(wf|kT)  is 


G0(w-»  =' 


I e Qn(<o  z,£)  dz 

3 

f f -JhT'i!  JkTl’ 

/ / e Qn("i£>  <fe'd£ 


qeil  (4  15) 


The  optimum  beam  pattern  is  given  by 


90(<ok|iiT)  = 


f f 'JJSt  — Jk-f 

/ / e Qn(w:z,ne  aze 



f f *T’2'  , , JkT-r 

J I e Qn(«z.i  )e  dz.  d 
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Note  that  g^cck-plk-p)  is  unity  This  last  equation  is  particularly  important  in  our  subsequent 
analyses  where  we  use  it  in  contrast  to  some  of  the  more  conventional  beam  patterns  discussed 
earlier  Equations  4 15  and  4.16  constitute  the  solution  for  the  optimum  array  weighting  and 
the  optimum  beam  pattern.  The  fundamental  issue  in  determining  a solution  is  finding 
Qn(co  z,f).  In  general  this  is  a difficult  problem  Fortunately,  for  the  examples  of  interest 
we  can  determine  the  inverse  kernel. 

The  only  remaining  issue  is  to  find  an  expression  for  the  output  noise  power  Ogtcjlkj) 
of  the  optimum  beam  pattern  We  simply  substitute  our  solution  into  Eq  4 2 Using  the 
properties  of  the  inverse  kernel  we  obtain 


-I 


2 

a 

o 


-jkT‘z  jk-f 

e Qn(tc  z£)e  dz  d£ 


(4  17) 


We  should  contrast  this  to  the  output  of  an  arbitrary  beam  pattern 


onto  | kT) ! 


F// 

nJn 


G(cc:z|kx)  Sn(u>z,f)  G*(to:f|kT)  dzdjj 


(4  18) 


The  term  (<o|kT)  is  the  minimum  value  that  a2(cj  |kj)  can  obtain  over  all  cho-ces  of  G(cj:z|kT). 


4.2  HIGH  RESOLUTION  ESTIMATES  OF  THE  FREQUENCY  WAVE  NUMBER 
FUNCTION  P(co:k) 

One  of  the  most  difficult  aspects  in  applying  techmqueswhich  utilize  the  process  statistics 
is  determining  the  actual  statistics  themselves  There  are  a number  of  ways  to  do  this,  and  a 
large  number  involve  measuring  sample  covariances  of  the  observed  signals  In  array  processing 
work,  this  is  typically  done  as  a function  of  frequency.  At  each  of  the  anay  locations,  a data 
segment  is  first  transformed  and  then  used  as  in  the  calculation  of  the  estimate  of  the 
covariance,  l e.,  we  have 

1 N 

&co:z,p  = - ]T  Rn«o:z)  Rn*(to £ (4  , 9a) 

where 
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r(n+I)T=»AT 


r(t  z)  e‘JGJt  dt 


The  important  point  is  that  one  obtains  S(<o^,f)  only  across  the  aperture,  or  array  The  problem 

A 

now  is  to  convert  this  into  an  estimate  P(o;  k)  Since  the  aperture  is  always  finite  we  can  not 

A "A 

simply  transform  the  estimate  S(co  z£)  and  obtain  P(w  k),  no  matter  how  accurate  our  co- 
variance estimate  is.  This  has  led  in  the  temporal  theory  to  window  trade  offs  and  the  in- 
direct, or  Blackman-Tukey,  spectral  estimation  procedures  [19]  Several  authors  have 
brought  optimum  array  processing  to  bear,  and  it  has  even  had  its  impact  on  the  estimation 
of  purely  temporal  processes  [18,  20-221 

We  observe  that  if  we  were  to  measure  the  noise  power  originating  from  direction  ky  we 
would  direct  a beam  in  that  direction  and  try  to  suppress  all  other  contributions  of  the  noise 
field  as  best  that  we  could  This  is  just  what  the  optimum  beam  pattern  does  Now  if  we 
measure  the  output  noise  power,  a truly  ideal  beam  pattern  would  yield  only  the  noise  power 
stemming  from  direction  ky,  since  it  responds  to  signals  only  with  that  wave  number.  The  high 
resolution  method  simply  asserts  that  OgCculky)  is  the  best  we  can  do  m this  respect,  and  it 
represents  the  power  in  the  signal  field  with  wave  vector  with  the  remaindei  of  it  suppressed 
in  an  optimal  manner,  i e.,  it  is  an  estimate  of  Pn(co  ky)  Consequently,  to  find  this  high  reso- 
lution  estimate,  Pn(wkp)  we  find  Qn(cj  z,f)  from  the  sample  covariance  and  compute 
P(cj  kp)  as 


kp)  = 


■ii 


-iJsT‘1  A JJi'r'i 

Qn(wz,f)e  dzdf 


In  essence,  this  procedure  makes  all  the  optimum  trade  offs  in  selecting  the  best  estimate  for 
finding  $(«  kp)  Consequently,  in  relating  our  subsequent  work  to  using  i’fwrjs)  for  high 
resolution  estimates  for  Pfw  k)  we  should  examine  how  close  (w  kp)  corresponds  to 
P(oo.kp). 

One  of  the  less  satisfying  aspects  of  using  this  high  resolution  measurement  method 
is  the  difficulty  in  analyzing  its  estimation  accuracy  This  contrasts  to  the  classical  procedures 
where  the  effects  of  aperture  extent  and  shadings  employed  can  be  analyzed  in  terms  of 
confidence  intervals.  The  nonlinearity  introduced  by  the  inversion  operation  for  Qn(w  z,J) 
leads  to  complicated  analyses. 


4.3  FREQUENCY  WAVE  NUMBER  ANALYSIS  FOR  OPTIMUM  WAVE  NUMBER  RESPONSES 

Before  discussing  specific  geometries,  we  reexamine  briefly  our  analysis  of  optimum 
beam  patterns  in  terms  of  a frequency-wave  number  representation  for  the  signals  From  Eq  4.2 
we  want  to  minimize  the  output  noise  power  subject  to  the  unity  response  in  the  direction  of 


3 


the  target  If  we  have  a noise  field  consisting  of  a propagating  component  plus  sensor  noise, 
we  have 


Sn(co:z,f)=  J f Pn(co:k)eJ-  <-^)  — 5=^ 

•L  J,„  (2ir)N 


dk 

^ + No«> 


We  substitute  this  into  Eq.  4 2 and  obtain 


cl  (wlkT) 


■in 

No  / |G0(. 


|go(“-3s.|kT)rPn(“li)- 


dk 

(2ir)N 


+ N~  I l/"'0(o):z | k-j-) | 2 dz 
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(4  21) 


(4  22a) 


subject  to  the  constraint 

tI— = 1 


(4.22b) 


When  we  interpret  the  design  of  optimum  beam  patterns  m the  frequency  wave 
number  domain,  we  can  observe  that  we  want  to  minimize  the  overlap  of  the  frequency  wave 
number  function  and  the  beam  patterns  and  the  integral  of  the  magnitude  squared  of  the 
weighting  function  across  the  aperture  subject  to  the  target  constraint.  Generally,  minimizing 
the  two  terms  in  Eq.  4.22a  involves  a trade  off.  Because  of  the  constraints  of  the  array 
geometry,  minimizing  the  first  term  introduces  high  sidelobe  levels.  Smce  we  have 


l 


IGto  zlkT^dJg^'^r 

An 


(4  23)* 


sidelobes  which  are  in  excess  of  unity  increase  the  second  term  using  this  inequality.  Generally, 
any  increase  in  the  sidelobe  level  increases  the  contribution  of  the  aperture,  or  sensor  noise, 
to  the  output  noise  level.  As  a result  these  two  effects  must  be  traded  off  until  the  optimum 
compromise  for  minimizing  the  total  output  noise  level  is  obtained. 


1 It  would  be  very  usefltl  if  we  could  use  a general  form  of  Parseval’s  equality,  Le , 

Jic(wl)l2dl  = JJJ 


Unfortunately,  the  integral  on  the  right  fails  to  converge  m many  cases 


...  4.4  BEAM  PATTERNS  FOR  SPATIALLY  SPREAD  SIGNALS 

I 

In  many  physical  situations,  the  signals  which  are  the  object  of  the  optimum  array 
processing,  or  beamformer,  are  not  true  plane  waves  in  that  their  frequency  wave  number 
function  is  not  impulsive  This  may  occur  in  a multipath  environment  in  which  local 
spreading  is  due  to  scattering,  and  gross  spreading  is  due  to  separate  travel  paths  Figure  4-2 
illustrates  a possible  propagation  situation  leading  to  a spatially  spread  signal  situation  The 
most  evident  effect  observable  at  the  receiver  is  that  the  magnitude  of  the  signal  covan  e 
decays  as  a function  of  the  separation  of  array  location 

We  model  the  signal  as  a general  space-time  process  s(t  z)  with  a spectral  covariance 
Ss(<o:z,f).  In  the  special  case  of  a plane  wave  structure,  we  have 


jkT-(zi) 

Ss(toz,f)  = S0(w)e  (4  24) 

with  an  impulsive  frequency  wave  number  function  In  many  situations  concerning  spatially 
spread  signals,  the  frequency  wave  number  function  occupies  a narrowband  about  k-p.  The 
principal  difference  m analyzing  spatially  spread  signals  involves  their  representation  across 
the  receiving  aperture  In  the  plane  wave  situation,  the  statistical  representation  in  terms  of  a 
Karhunen-Loeve  expansion  is  one  dimensional,  i.e  , only  one  eigenfunction  is  required  in  the 
representation.  In  the  spatially  spread  situation,  more  than  one  eigenfunction  is  required  For 
example,  with  gross  multipath  spreading  we  need  at  least  one  eigenfunction  per  path,  while 
with  local  scattering  there  is  an  upper  bound,  which  is  determined  by  the  area  of  the  array  and 
extent  of  the  spreading,  on  the  number  of  eigenfunctions  required  to  represent  effectively  the 
signal  across  the  aperture  In  most  cases  of  interest  the  number  required  is  small  so  that  we  have 
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(4  25a) 


and 


M 

Ss(wz,{)=  ^2  Xm(nw0)4.m(nw0z)  4>^<ncoop  (4  25b) 

m=l 


with 


E[Sn,(n<*>0)s*k(&*>o)]  - Am(nco0)  5ng 


(4  25c) 
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Figure  4-2.  An  example  of  spatial  spreading  due  to  global  path  structure  and  local 
scattering  of  thermal  structure. 
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The  structuring  of  the  design  of  the  optimum  beam  pattern  is  somewhat  more  subtle 
than  that  for  a single  plane  wave  in  that  we  no  longer  want  simply  to  direct  a beam  at  a 
specific  wave  number.  In  our  formulation,  we  generalize  this  concept  by  designing  a set  of 
beams,  one  for  each  of  the  M eigenfunctions  in  the  spatiai  representation  of  the  signal  across 
the  aperture  For  each  of  these  beams  we  require  that  the  response  due  to  the  ambient  noise 
be  a minimum  under  the  constraint  that  there  is  an  undistorted  response  to  the  eigenfunctions 
for  all  possible  signals.  Just  as  we  noted  earlier,  we  could  alternatively  formulate  the  opti- 
mization in  terms  of  a detection  theory  problem  with  the  resulting  analysis  leading  to  the  same 
spatial  structure  that  follows  here 
We  have 


o.2(cj)  = 


LL 


G(<o  z|0j)  Sn(co  z,ff)  G(oj  z|0j)  dzdf 


(4  26a) 


nJn 


with 


[ V 

I G(<o  z|0j)  > 

•'O  1 


sm(<o)  $m(“ d£=  sft*#,) 


(4  26b) 


SI  m=l 


for  all 


S(col0,)| 


M 

i=l 


The  optimization  procedure  is  quite  similar  to  that  which  we  did  earlier  so  we  just  indicate 
the  results.  We  define 


<j>(w  z)  = 


4>l(w.g) 


4-M(a;z) 


y ze  SI 


(4.27a) 


£0(coz)  = 


Glc  G0(w  z|0,) 


GMq  G0(«  zl0MJ 


(4  27b) 
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The  optimum  set  of  beam  patterns  are  given  by 

G0(io:J)=  AJ*(w)  / <t>^7ojz)  Qn(to:z,f)  dz feJ2 

^0(w:k)  = A'^to)  I *(cuz)  Qn(eo.z,f)eJ“  - dz  df 
T2 

The  mean  square  output  level  due  to  the  signal  component  is  given  by 
e["  /|Gm(cj:z)  S(cj  z)  dz|  2]  = \m(u) 


while  the  output  due  to  the  noise  components  are  given  m matrix  form  by 
o2  (<o)  = E G0(co  i)N(ui)  d^^  G0(w  f)  N(oi:£)  djj 

-f  %(“)  1 

|A11«o)Ajj(to)| 


In  summary,  the  signal  components  each  pass  through  the  filter  undistorted,  i e.,  they 
are  unbiased.  These  components  have  a diagonal  covariance  matrix  as  given  by  Eq.  4.29a  The 
noise  components  are  correlated  for  each  beam  and  their  noise  covariance  is  given  by  Eq.  4 29b 
Various  means  of  combining  the  signal  can  be  considered  according  to  the  criterion  imposed. 
We  consider  this  further  in  Section  8. 
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PART  3 - PROCESSING  METHODS 


5.  LINEAR  ARRAYS 


The  most  important  aperture  is  the  linear,  or  line,  array  It  is  the  simplest  structure, 
yet  the  analysis  of  .t  yields  many  significant  results  which  more  complicated  structures 
change  only  in  detail,  but  not  in  substance.  Much  of  the  analysis  of  this  array  parallels  that 
of  temporal  processes  observed  over  a finite  time  duration  Often  this  is  a distinct  advantage 
since  only  dual  results  need  be  developed  This  can  also  prove  to  be  somewhat  of  a liability  as 
several  important  aspects  of  the  analysis  of  array  processing  do  not  appear  until  more 
complicated  geometries  are  introduced.  For  our  purposes  we  consider  conventional  and 
optimuri  array  patterns  when  operating  m different  types  of  noise  fields  Then  we  consider 
some  general  feature  of  the  processing,  such  as  wave  number  analysis,  superdirectivity,  and 
null  placement  where  we  use  the  linear  array  as  a discussion  vehicle 

As  introduced  in  Section  3,  the  receiving  aperture  consists  of  a line  of  length  L, 
centered  at  the  origin  (for  convenience  only),  and  oriented  with  a tangent  vector  in  the 
direction  specified  by  the  unit  vector  §a. 

« = {z  = Caa,  |Cl  < L/2^  (5.1) 

as  illustrated  in  Figure  5-1 . 

The  beam  pattern  of  the  conventional  array  directed  in  wave  number  direction  k-j- 
is  found  to  be  (see  Eq.  3 •)) 


i [m 

g(w.k)=—  I < 

L <Ll2 


,i(k-kT)-aa« 


dS 


= sinc|(k-kT)-aa~J 


(5.2) 


When  we  apply  a triangular  shading  of  the  form 


I SI  < L/2 


(5  3) 
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Figure  5-1.  Linear  array  of  length  L,  orientation 

we  obtain  a beam  pattern,  or  wave  number  response,  of  the  form  (see  Eq  3 12) 

j(u  k)  = sine2  |(k-kT)-aa  (5  4) 

Both  of  these  are  standard  array  patterns.  Others  can  be  developed,  for  example,  the  split 
array  illustrated  in  Figure  3-3  We  are  primarily  interested  in  optimum  array  patterns  and 
their  performance  o2(co  k-p)  for  various  noise  Fields  We  contrast  these  optimum  patterns  and 
performance  with  those  obtained  conventionally  to  determine  how  they  operate  in  obtaining 
their  improved  performance  We  also  compare  the  performance  to  determine  whether  the 
additional  eftort  in  optional  processing  is  worth  while 

5.1  LINEAR  ARRAYS  AND  A SINGLE  DIRECTIONAL  NOISE  SOURCE 


The  simplest  nontrivial  problem  which  one  can  analyze  for  an  optimum  processor 
is  the  linear  array  combatting  a directional  noise  Field  as  shown  in  Figure  5-2.  This  problem,  in 
spite  of  its  apparent  simplicity,  yields  a surprising  number  of  useful  results  and  an  insight  into 
the  general  aspects  of  optimum  array  processing. 

We  have  a directional  noise  source  of  power  level  located  at  k„  - (co/c)  an  = 

(2ir/\)gn.  Consequently,  the  colored  component  of  the  noise  is  a plane  wave 


'>c<U>  = no('-~) 


i v j- 


From  our  earlier  analysis  in  Section  2,  the  temporal  frequency  spatial  covariance  function 
of  this  component  is  given  by 

S_  Mz,f)  = S.  (co)e  , (5  6) 

IlC no 

while  the  corresponding  frequency  wave  number  function  is  given  by 

Pn  (wk)  = Sn  (<u)  5(k-kn)  (5  7) 

“c  o 

Wc  also  assume  that  a white  noise  component  due  to  sensor  noise  is  present  In  light 
of  our  earlier  comments  regarding  this  class  of  noise  process,  we  have 


Sw(<o.z,f)  = N0  8n  (zi) 


(5.8) 


The  total  noise  field  is  specified  by 

Sn«oz,f)  = SnoMe  +N06n(z^) 


C9) 
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We  are  interested  in  determining  the  optimum  processor  and  its  performance  for 
detecting  a signal  at  k-j-  The  first  uep  in  designing  this  optimal  processor  is  to  find  the 
inverse  kernel.  For  directional  sources  this  is  particularly  simple  because  the  discrete  spatial 
frequencies  of  the  noise  lead  to  a separable  kernel  with  a finite  number  of  eigenfunctions 
The  Karhunen-Loeve  integral  equation  is  given  by 
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#*»)  = *««,),  ISj  I < L/2 


(5  10) 


The  single  non-tnvial  solution  is 
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(S  11a) 
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The  inverse  kernel  is  therefore  given  by 
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From  Eq.  4 1 we  have  that  the  aperture  pattern,  or  weighting,  is  given  by 
G(<o  z|kj 
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Similarly,  the  beam  pattern  can  be  found  using  Eq  4 1 2 This  yields 
g(co.k|j£y) 

sine  [(k-kT)-aa  ISn-(-^+-No  sine  [(k„-kT)-aa  ±J  sme  [(k^-a,  |) 
LSn  (co) 


^°n  K ; r il 

-Ls;(J)+N0-sinc?[(w^i] 


(5  14) 


The  operation  of  the  optimum  array  can  be  interpreted  as  a null  placement  operation 
which  depends  upon  the  relative  strengths  of  the  directional  and  white  noise  If 


LS„o«o) 


<<  1 


then 


g(to  k|  kT)  = sine  Jfk-kjVaj,  y 


(5  15) 


This  is  the  conventional  beam  pattern  which  has  a uniform  weighting  and  is  optimal  when 
only  white  noise  is  present.  Conversely  when 


LSn  (to) 


>>1 


then 


9(to.k|kT)  = 


r Li 

sine  (k-k j)  -aa  — - sine 

(-n'-T> -a  j]  Slnc  i— -ni’-a  7] 

1 - sinc^ 

(fcn"-Ti’-a  J 

(5  16) 
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In  the  direction  of  noise  source,  this  beam  pattern  has  the  response  level 


g(^knliT>  ~0’  (5  17) 

that  is,  a null  is  placed  there  Placing  this  null  is  done  at  the  expense  of  increasing  the  output 
noise  due  to  the  white  noise  component  Specifically,  in  this  situation  the  white  noise  is 
increased  by  a factor  of  approximately 


l/(l-sinc2  (kn-kT)\aay]j2 


over  the  minimum  level  obtained  with  a uniformly  weighted  array  If  the  noise  and  target  are 
widely  separated,  i e , 


(kn-kT)-aay  >>2* 


the  null  placement  has  no  effect  However,  if 


(kn-kT)-2aT  - 0, 

then  the  white  noise  response  level  is  increased  by  a factor  of  approximately 


[l  L l"4 

3(kn-kT)'Ja-j  >• 


which  is  a very  dramatic  increase.  Consequently,  in  designing  an-ay  weightings  with  the  target 
and  noise  close  together,  one  follows  a null  placement  procedure  up  to  the  point  that  response 
due  to  the  white  noise  becomes  dominant. 

Figure  5-3  illustrates  a situation  in  which  the  directional  noise  is  dominant  At  separa- 
tions of  30°  the  beam  pattern  resembles  that  of  a uniformly  weighted  array.  At  10°  and  3° 
the  pattern  has  very  strong  nulls  near  the  noise  term  Finally,  at  1°  the  null  can  no  longer  be 
maintained  for  the  white  noise  has  become  the  dominant  consideration  due  to  the  large  side- 
lobes  appearing  in  maintaining  the  null  Note  that  in  the  optimum  beam  pattern  response,  the 
peak  does  not  fall  right  at  the  target  location  when  the  noise  direction  is  close  to  that  of  the 
target  Effectively,  it  is  displaced  over  to  facilitate  placing  the  null.  We  shall  return  to  this  in 
subsequent  discussions.  The  performance  or  noise  power  output  o2(co|kj)  for  this  array  and 
noise  field  can  be  found  using  Eqs  4 17  and  5.12.  We  have 


83 


Optimum  array  beam  pattern  for  a directional  plus  white  noise  environment  — rr 1 20 
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We  see  that  all  spatial  aspects  of  the  problem  are  incorporated  in  the  parameter 


A 

p - sine 


j(kn-kT)-aa  y 


while  the  strength  of  the  interfering  noise  appears  in  the  ratio 


LSn  (co) 
“o 


In  Figure  5-4  we  have  plotted  the  output  noise  power  o';  (co  | k-j-)  versus  p relative  to  a 
reference  level  of  N^L 

In  interpreting  this  figure  we  first  note  that  when  the  target  and  noise  are  widely 
separated,  i.e.,  p 3 0 


2 "o 

oMkT)  « — 


while  when  they  are  coincident,  i.e.,  p =*  1 


o^(«|k, 


No  //  LS"o(w)V 


In  terms  of  the  performance  using  optimum  array  design,  the  most  important  parameter  is 
the  resolution,  that  is,  how  close  to  the  directional  noise  source  a target  can  be  before  the 
noise  power  increases  significantly.  Alternatively,  when  using  “high  resolution”  methods  for 
noise  field  measurements,  one  wants  to  know  the  resolution  of  these  methods 
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Since  the  spatial  aspects  are  incorporated  in  the  parameter  p,  the  resolution  is 
determined  by  this  parameter.  However,  the  value  of  p at  which  the  interference  becomes 


significant  is  dependent  upon  the  parameter  L S_  (a>)/N„,  as  stronger  interference  degrades 

• I1Q  u 


the  performance  faster. 

For  discussion  purposes  let  us  define  the  value  of  p which  leads  to  a 3dB  performance 
loss  by  P3 or  a 3dB  increase  above  the  background  white  noise.  We  have 


LSn  (co) 
o 


1+- 


(’  -P1  3d b) 


= 2 


(5  24a) 


P 3dB  = 


LS„  (u) 
“o 


N„ 


LS.  (oj) 


(5  24b) 


N„ 


This  value  ranges  from  one,  as  a 3dB  reduction  never  occurs  for  noises  with  a level  less  than 
(L/N0)‘*.  to  l/\/T for  strong  interference  Consequently,  even  for  strong  interference  one 
always  can  obtain  a performance  with  3dB  of  the  limiting  white  noise  level  providing  the 
spatial  geometry  yields  a p less  than  1/y /~T.  Equivalently,  the  resolution  can  be  defined  by 


sine  (krkn)'%~  > >/v^" 


(5.25a) 


The  response  to  a directional  signal  will  be  less  than  3dB  above  background  outside  of  this 
region.  Solving  this  yields 


(krkn)-iiay  <1.39 


(5  25b) 


For  broadside  geometries 


2 cosl 
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(5  26a) 


87 


S2S<3cyifl'r««.'  id  V.V  ■«'  I‘V  ■'  * 


r<;.j 


m 

$ 


KM 


’■A 


Mi 


where  is  the  length  of  the  array  in  wavelengths.  Figure  5-5  illustrates  a typical  broadside  array 
configuration  for  directional  noise  and  target. 
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ANGLE 


TARGET 
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Figure  5-5.  Broadside  array  configuration  for  directional  noise  and  target 
For  an  array  larger  than  one  half  a wavelength,  or  > 1/2,  we  have 


.444  j 25.4° 

0 =*  — - — rad  

S L» 


This  should  be  the  “critical”  separation  angle  for  a linear  array  using  “high  resolution” 
methods  in  the  broadside  case.  (Note  that  we  assumed  no  other  sources  to  be  present  in 
the  noise  field ) 

As  the  array  is  “beamed”  off  broadside,  the  effective  resolution  becomes  smaller. 
Let  us  consider  the  endfire  situation  as  represented  in  Figure  5-6.  Here  we  have 
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Figure  5-6.  Endfire  array  configuration  for  directional  noise  and  target. 


For  large  arrays  with  > 1 this  becomes 

0S  “ .942VTX  rad  = 53.6°/v/T^  (5  29) 

Observe  that  the  effective  resolution  of  an  array  changes  from  Ly1  dependence  at  broadside 
to  one  of  Ly1^2  for  an  endfire  configuration 

Before  illustrating  the  performance  of  the  optimum  beam  pattern,  we  discuss  the 
performance  of  some  suboptimum  beam  patterns.  In  this  context,  we  consider  two  commonly 
used  array  shadings,  uniform,  or  conventional,  and  triangular;  and  then  we  contrast  their 
performance  to  that  of  the  optimal  array  processor  just  analyzed. 

For  a target  with  wave  number  lsj,  the  conventional  or  uniformly  weighted,  linear 
array  has  a weighting 


Gc(co  X)=j-  e'jkT'aa*.  |S|  >k 


As  discussed  earlier  in  Section  3,  this  weighting  produces  a beam  pattern 
9c(w:k)  = sine  |(k-kT)-ga  -yj  = p 


(5.30) 


(531) 


The  mean  square  noise  output,  or  the  performance  for  a single  directional  signal  and  white 
noise  environment  follows  directly  from  Eq  4 18.  We  have 


2("l  kT)  = S„o  (<o)  sire2  jf^-kx)-^  -yj  +-^ 
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A triangular  shaded  linear  array  directed  to  a target  at  wave  number  kT  has  a weighting 

^ 2 r 2 ... 


Gj(uz)  = y-  [''Y  |{|]e  -T_a,.  |8|<i 


A similar  analysis  shows  that  the  beam  pattern  generated  is  given  by 
gj(<o  k)  = sine2  [(k-kg.)-^  -^J 

while  the  mean  square  noise  output,  or  the  performance,  versus  this  noise  field  is 
2.  .....  . a r.  LI  4/No\ 


(“I*T>  = Sn0<w)  s,nc4  [^n^T)'aa  j]  + ^(~j 
m L Sn  (w) 

o no  a r LI  4 

= — n sme' ’[(kn-kT^a  +T  (S  35) 


In  Figures  5-7  through  5-9  we  have  plotted  the  performance,  or  noise  output  for  the 
optimum,  uniform,  and  triangular  weighted  arrays  for  van'  ms  ratios  of  the  directional  signal 
power  to  background  white  noise  level  We  see  that  the  optimum  array  weighting  combines  the 
desirable  features  of  both  the  uniform  and  triangular  arrays  Its  enhanced  performance  is 
particularly  significant  when  the  ratio  of  directional  to  white  background  noise  is  high  In 
this  situation  optimum  design  of  the  array  weighting  is  worthwhile  providing  the  results  are 
not  unduly  sensitive  to  the  model  assumptions 

The  next  observation  concerns  the  sidelobe  structure  of  the  noise  power  response 

versus  kT.  When  the  directional  noise  is  high,  i.e.,  L S (w)/N.  > > 1 , then  the  sidelobe 

o w 

level  of  o 2 (oi|  kT)  for  a uniformly  weighted  array  goes  as  [2/(2n+I ) ir]2,  or  from  1 3 dB 
down,  while  n2  (<o|  kT)  for  a triangular  weighted  array  goes  as  (2/(2n+l  )jr)  4,  or  from  26 
dB  down 
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Performance  of  optimum,  uniform,  and  triangular  linear  arrays  vs  target  position 
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With  the  optimum  array  weighting,  Eq.  5.20  implies  that  the  maximum  sidelobe  level  of 
o2(ai|kT)  goes  as 


LSn  (<o) 
o 


\ [<2n+iH ) 


1.05/ 


m 


This  indicates  that  the  sidelobe  level  is  inversely  dependent  upon  the  directional  noise. 

For  example,  in  Figure  5-7  the  sidelobe  level  is  down  approximately  21  dB.  Naturally,  there 
is  a premium  to  pay  for  this.  As  we  have  previously  discussed,  the  optimum  array  shading  is 
a null  placement  operation  in  the  situation  of  strong  directional  noise.  Unfortunately,  this 
becomes  a very  sensitive  situation  requiring  precise  knowledge  of  the  direction  of  the  noise 
which  is  typically  not  available  in  an  experimental  situation.  This  sensitivity  to  direction  is 
quite  evident  when  one  examines  the  beam  pattern  As  indicated  in  Figure  5-3  the  beam 
pattern  has  a very  high  sidelobe  structure  such  that  a small  amount  of  error  between  the  actual 
noise  direction  and  the  nulflocation  would  lead  to  a very  high  response  to  this  noise. 
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5.2  LINEAR  ARRAYS  WITH  TWO-DIRECTIONAL  NOISE  SOURCES 


The  optimum  array  achieves  an  enhanced  processing  gam  by  effectively  displacing 
its  beam  pattern  to  null  the  noise  If  the  target  is  located  between  two  noise  sources  this  pro- 
cedure is  not  viable,  for  one  would  have  a strong  response  to  the  directional  source  due  to 
stdelobes  on  the  other  side  of  the  beam  In  these  situations,  optimum  processing  can  still 
improve  the  array  performance  although,  depending  upon  the  geometry,  not  as  much  as  in 
the  single  noise  source.  A particular  example  of  two-directional  noise  source  field,  as  shown 
in  Figure  5-10,  introduces  important  issues  in  the  limited  “high  resolution”  mapping  methods 
and  superdirectivity  We  consider  the  issue  of  determining  when  two  noise  sources  are 
rcsoluable,  or,  from  an  optimum  beam  pattern  viewpoint,  we  determine  if  we  can  create  a 
beam  that  can  observe  a signal  located  between  two  noise  sources  This  problem  is  a very 
special  case  since  we  can  obtain  closed-form  expressions  for  a noise  field  composed  of  an 
arbitrary  number  of  sources  Insight,  however,  is  lost  in  the  matrix  manipulations  which  are 
necessary.  Moreover,  this  problem  is  simple  enough  that  we  can  obtain  considerable  physical 
insight  by  carrying  out  a detailed  analysis  if  we  assume  that  the  noise  sources  have  the  same 
temporal  spectra,  then  we  have  a spectral  covariance  function  of  the  form 
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Sn  (cu:k,f)  = 


Vw> 


( -JiKftei:) 

\c  ••  i 


(5  36) 


Substituting  this  into  the  spatial  Karhunen-Loeve  Eq  3.53,  we  can  obtain  a matrix  equation 
since  the  kernel  is  separable,  i.e  , it  contains  a finite  number  of  nonzero  eigenvalues  For  a 
linear  array  of  length  L,  oriented  with  array  tangent  aa,  we  obtain  the  eigenvalues  and 
eigenfunctions 


=L 


\(W) 


(1  + Ptb.j-Js.2)) 


(5.37a) 
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(5.37b) 
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Figure  5-10.  Geometry  for  two  directional  noise  source  field. 


where 


i«l  < L/2 


(5.37e) 


(5  370 


We  observe  two  features  characteristic  of  the  eigenvalues  and  eigenfunctions  of  a 
directional  noise  field.  The  eigenfunctions  consist  of  a linear  combination  of  the  functions 

JV.»2 

e 1 while  the  eigenvalues  are  determined  by  how  similar  these  functions  are  across  the 
aperture.  In  this  example,  if  the  signals  are  orthogonal,  i.e.,  ) = 0,  then  the  eigenvalues 

are  equal  and  reflect  the  noise  power  of  ~ Sn  (co)  from  the  individual  directions.  If  they  are 
very  similar,  i.e.,  p(}i.2— lii ) “ 1 , then  there  is  one  dominant  eigenvalue  of  L Sn  (co)  which 
represents  the  situation  that  almost  all  the  noise  power  originates  in  a direction  which  is  the 
average  of  the  two  source  locations. 
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We  can  determine  the  optimum  beam  pattern  and  its  output  noise  power  using  these 
expressions  and  assuming  that  sensor  noise  with  level  N0  is  present.  Since  the  total  noise  power 
is  also  the  term  that  normalizes  the  beam  pattern  we  determine  it  first  We  have  after  some 
manipulation 


tM<o|kT)= 


-}kj‘z  jkT-_f 

e Q(o>z,£)e  - dz^dj 


L l 2 (l+(3[l+p(k,-k2)]) 


(p(k]-kT)  + p(k2-kT))2 


4 ' 

g(w : k |kT)  = |o2( w | kT)j jpftkx)  (5.3 

4 TTJTi'-fpdt'^ji  W^i)  + ^2)Hp<Krlii)  + Mkrk2)) 


T l HU  -Pfkr^)))  0’atil > - *»2»  <p(*t*i > -p(kr*2»  , 


P = LSno(<o)/2N0 

While  complicated,  these  expressions  can  be  interpreted.  If  p(kj-k2)  “ 0,  then  the  beam 
pattern  has  nulls  of  approxunately 

P(ki-hT) r ! 5 

1 +010  - P2(JSj -kx)  — PZ(!S2~^T)1 

at  each  kt,  i.e.,  stronger  noise  sources  produce  deeper  nulls.  If  we  have  p(kj-k2)“  1,  the 
beam  pattern  and  response  resemble  that  for  a single  noise  source. 

In  Figures  5-1 1 and  S-I2,  we  have  illustrated  these  results  for  an  array  and  noise  source 
geometry.  We  can  observe  that  for  2°  and  6°  separations,  the  noise  power  output  resembles 
that  of  a single  source.  Here  p(kj-k2)  “ .99  and  .94  respectively.  At  20°  and  60°  separations, 
the  effects  of  the  two  noise  sources  are  resolvable  with  each  source  producing  a response 
approximately  equal  to  (3.  If  we  examine  the  conventional  beam  pattern  response  in  Figures  5-1 3 
and  5-14,  for  comparison,  we  find  much  the  same  behavior  for  the  peaks;  however,  the  sidelobe 
level  is  very  high  and  could  easily  be  attributed  to  other  sources  of  a lower  power. 
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The  patterns  illustrated  in  Figures  5-1 5 through  5-18  demonstrate  how  the  output  noise 
power  is  minimized.  For  a 2°  separation  the  beam  pattern  differs  substantially  from  the  conven- 
tional one  only  when  0=60  Even  a minor  reduction  of  the  beam  pattern  in  the  noise  directions 
is  worthwhile  when  the  directional  noise  level  is  high  Note  that  a higher  sidelobe  level  is  intro- 
duced to  accomplish  this  As  the  separation  between  the  two  noise  sources  increases,  all  the 
beam  patterns  differ  from  the  conventional  beam  pattern  The  pattern  for  0=60  has  very  high 
sidelobe  in  excess  of  the  main  lobe  for  separations  of  6°  and  20°  This  behavior  contrasts  with 
that  for  a single  source  In  that  situation  the  closer  the  separation,  the  higher  the  sidelobe  In 
this  case,  the  beam  pattern  must  minimize  the  effect  of  two  sources,  consequently,  putting  the 
target  on  the  skirt  of  a mam  lobe  is  not  a viable  alternative  At  wide  separations  the  large  side- 
lobe disappears,  and  only  the  field  with  a very  high  noise  source  places  a deep  null  in  the  noise 
direction  In  all  cases  we  can  see  the  gross  behavior  which  governs  the  beam  pattern  for  a 
directional  noise  field  We  try  to  place  nulls  in  the  direction  of  the  noise  sources  while  not  dis- 
turbing the  nulls  at  other  noise  sources  and  minimizing  the  total  sidelobe  level  so  as  to  reduce 
the  effects  of  the  sensor  noise 

We  can  interpret  the  output  noise  power  in  more  detail  by  manipulating  Eq  5.38 
This  also  gives  some  insight  into  the  optimum  resolution  of  the  two  noise  sources  We  can 


express  a^(cj|kj)  as 
o2(<o|kT) 


N0  010+0)  (P2(irlST>+P2(k2-3iT))-2PP(lil-k2)P(lii-liT)P(t2-kT)] 
L (/J+l)2-j52p2(krk2) 


The  second  component  in  the  brackets  represents  an  increase  due  to  the  directional  noises, 
which  is  analogous  to  the  term 

sme2  |(kn  -i.T) ' 3a  ^)L  Sn0(tJWNo 


l+LS„o«o)/N0 

for  the  single  directional  noise  case  (see  Eq.  5 20).  It  is  a quadratic  function  of  the  terms 
p(k i and  p( k2-ky );  consequently,  the  increase  in  output  noise  follows  elliptic  contours 
versus  these  parameters  These  ellipes  have  a major  axis  representing  the  geometry  for  maximum 
interference  and  a minor  axis  where  it  is  a minimum.  The  major  axis  is  specified  by 

POSj-kx)  = Pfkj-kr)  (5.41a) 
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Figure  5-15.  Optimum  beam  patterns  for  two-directional  noise  sources  at  0=90' 
0n=±l.o'>; 0-f=9O,  d>T=0°;and  (3=60.0,  6.0, 0.6. 
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ar-d  the  minor  axis  is 


P(k  1 -llT)  = -p(k2-iiT> 


(5.41  b) 


Since  p is  an  even  function  of  Jj,  the  maximum  noise  interference  occurs  when 

£)-£t  = Ji2'-T  (5.42a) 


or 


k]-kT  - <k2-k-j-) 


(5.42b) 


The  first  occurs  when  the  noise  sources  are  close  together  and  the  target  is  identical  to  the 
case  of  a single  directional  noise  source  twice  the  level.  In  the  second,  the  target  is  situated 
between  the  two  noise  sources,  which  is  the  geometry  that  we  illustrated  earlier. 

We  define  the  sources  to  be  resolvable  if 


2 

o 

o 


H£])  > 


(5.43) 


In  the  geometry  analyzed  previously,  this  requires 

1 +0(1  — (3)  sinc^  (2jrL^sin$n)  > 2sinc^  (rrL^sin^j) 

•[1+0(1-  sinc(2irLx  sin0n))J  (5.44) 

In  Figure  5-19  we  have  plotted  the  0 required  for  resolving  tne  two  sources  versus  the  separa- 
tion angle  between  them  for  = 3 For  example  once  they  are  closer  than  <pn  =*  3°,  a 0 of  at 
least  600  is  required. 
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Figure  5-19.  Required  SNR  for  separating  two-directional  noise  sources. 
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5.3  WAVE  NUMBER  ANALYSIS  FOR  LINEAR  ARRAY 


Before  considering  analyses  for  more  noise  fields,  we  examine  a frequency  wave  number 
approach  to  optimal  beam  pattern  design.  We  have  found  that  the  quantity  o^(w|kTl  specified 
the  noise  power  at  the  beamformer  output.  From  Section  4.2  (Eqs.  4.22a,  4.22b),  we  minimize 


2 

a 

o 


(w|kT)  = 


| g0<w:&  |St>1  2Pn<":fc> + No 


L 


|G0(co:z|kT)|  dz 


(5.45a) 


with 


So(«:kTlkT)=  1 (5.45b) 

For  linear  arrays  we  can  obtain  a particularly  useful  and  intuitive  approach  using  wave-number 
concepts.  This  was  apparently  used  first  by  Woodward  [23J  and  later  extended  for  analyzing 
superdirective  arrays  by  Francia  (24J . 


& 
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r > v . . - ? > 


Let  us  define  Gg(w:C)  to  be  a weighting  across  the  array  and 


ka=JS.*aa 


(5  46) 


We  have  for  the  beam  pattern 


g(co:k) = 


jk-aa2 

Gj(<j:2)e  d£  = gg(co:ka) 


(5  47) 


where  gg  is  the  Fourier  transform  with  respect  to  the  spatial  variable  of  Gj.  Consequently, 
only  the  projection  of  ka=k-aa  is  relevant  in  determining  the  beam  pattern  output.  This  is 
illustrated  in  Figure  5-20. 


Figure  5-20.  Planes  of  constant  response  for  a linear  array  with  orientation  ag  and  length  L. 
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Figure  5-21.  Optimum  beam  pattern  with  only  sensor  noise  present. 


this  one  would  cancel  the  response  in  the  noise  direction  with  a beam  of  opposite  sign.  We 
cannot  do  this,  however,  in  an  arbitrary  manner.  Since  array  has  finite  length,  we  cannot 
realize  an  arbitrary  <j(co:k).  The  sampling  theorem  says  that  on  the  average  we  can  specify  one 
sampled  value  per  distai.ee  of  2ir/L  independent  of  the  other  samples  and  still  realize  the 
function  with  a finite  length  array.  When  we  specify  these  points  as  g |co.n  ^ j,  we  can 
determine  any  intermediate  values  by  the  interpolation  formula  associated  with  the  sampling 
theorem,  i.e.. 


CO 

g«o:k|kT)=  Y) 
n=-°° 


g^cu:n-^|  sine  jfk-jsj)*; 


with 


g(oj:kT|  kT)=  ! 


(5.52) 
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If  we  want  to  place  a null  in  this  beam  pattern  we  must  adjust  the  values 

g (to  n =2-|  To  place  a null  for  a single  noise  source,  two  situations  can  occur  When  the  noise 
' *-  / / ■)  \ 
source  is  outside  the  mam  beam  region  the  null  can  be  placed  by  adjusting  the  g(co  n •=— 

/ 'a  \ *-  / 

nearest  it  so  as  to  cancel  the  main  beam  response  Only  a modest  value  of  g (to  n —■] is  needed, 
since  the  mam  beam  response  is  quite  low  Consequently,  the  added  Side  lobe  response  is  small 
If  the  noise  source  is  inside  the  main  beam,  one  must  cancel  a very  large  mam  beam  response 
value.  This  requires  a large  value  of  the  adjacent  g(to  n)  and  as  a result  very  iatge  sidelobes  are 
introduced.  Both  of  these  situations  are  illustrated  in  Figure  5-22.  The  amount  of  additional 
sidelobe  response  is  the  value  expressed  in  Eq.  5.51  m terms  of  the  effects  background  noise 
introduced.  When  noise  single  source  is  near  the  target,  the  response  function  appears  to  be 
displaced.  The  amount  of  displacement  is  governed  by  the  ratio  of  directional  to  white  noise 
and  the  proximity  between  the  target  and  directional  noise  source  This  is  the  dominant  effect 
with  a single  directional  noise  source 

When  several  directional  noise  sources  are  present,  the  situation  becomes  more 
complicated.  In  adjusting  the  coefficients  g^co:n  we  must  consider  not  only  the  side- 
lobe  level  but  also  the  mutual  interference  created  by  placing  the  various  nulls  If  the  sources 
are  widely  separated  relative  to  the  resolution  of  the  array,  the  interference  is  minimal  If  they 
are  not,  we  must  determine  the  exact  trade  offs  using  our  optimum  array  design  procedures 
For  ware  number  fields  which  are  not  directional,  the  imtuitive  appeal  of  Eq.  5 5 1 
is  still  useful  Basically  we  attempt  to  minimize  the  overlap  of  the  beam  pattern  and  the  wave 
number  function  for  the  noise  field  without  undue  effects  from  sensor  noise  As  we  will 
subsequently  discuss,  this  minimization  of  the  overlap  can  often  be  exploded  in  a manner 
unique  to  spatial  processing.  Essentially,  for  propagating  signals,  the  maximum  wave  number 
for  a given  to  is  <o/c,  or  2ir/A,  and  all  propagating  noise  could  be  eliminated  if  the  sidelobes 
were  confined  to  a region 


In  principle  this  can  be  done,  however,  it  requires  extremely  large  sidelobes  for  ka  outside 
this  region.  Since  the  output  due  to  sensor  noise  still  is  affected  in  this  region,  the  response 
becomes  extremely  noisy. 


m 


m*!-riri7vr.rjf±r. .'  -J ')- 1 


S.4  LINEAR  ARRAYS  IN  ISOTROPIC  (OR  OMNIDIRECTIONAL)  NOISE 


A common  model  for  ambient  noise  of  deep  ocean  noise  is  the  isotropic  noise  field 
distribution  discussed  in  Section  2.2.  The  spectral  covariance  function  for  the  model  is 


Sn(<o  Az)  = S„o(<o)  sine  ^ lAz  ij 


(5  53) 


Many  beam  patterns  are  designed,  and  the  figures  of  merit  attached  to  their  performance  are 
discussed  in  term  - of  their  effectiveness  in  combatting  this  particular  noise  Unfortunately  as 
we  will  see,  using  this  noise  field  alone  can  be  misleading  for  it  is  closely  relate)  to  super- 
directivity  and  singularity  issues  it  In  this  section,  we  briefly  discuss  how  such  noise  fields 
can  be  incorporated  in  optimum  array  design  by  using  results  established  in  communication 
theory  and  electromagnetics.  We  also  apply  our  frequeney-wave-number  analysis  to  discuss 
endfire  gain  in  linear  arrays.  Tins  same  analysis  extends  to  the  nonuniform  distributions  of 
noise  which  we  discussed  in  Section  2.2. 

We  consider  the  noise  field  to  be  composed  of  an  isotropic  and  a white  noise  com- 
ponent such  that 


Sn(«:Az)  = S„o(u)  sine  | Az  | j + NQ  Sjj  (Az) 


(5.54) 


The  frequency  wave  number  function  in  three-dimensional  space  for  isotropic  noise  consists 
of  a sphere  of  impulses;  however,  when  one  integrates  all  those  parts  of  the  sphere  which  have 
the  same  wave  number  projection  k"aa,  a rectangular  function  of  X*a  results.  Alternatively, 
the  one-dimensional  transform  of  Eq.  5.54  with  respect  to  the  spatial  frequency  as  observed 
by  the  array,  i.e.,  ka,  yields  the  distribution  as  shown  in  Figure  5-23. 1 

For  a linear  array,  the  integral  equation  that  specifies  the  optimum  array  weighting 
(see  Eq.  4.6)  is 


/L/2 

S„  (oj)  sinc( — j G0(w:82l isf)  ^2  + No  Go'CJ:Cll  -T^ 
n 0 ' c ' 


-L/2 


(5.55) 


jkT-aJi 

= -(3{<i})e  l«jl<L 


:ka)  = J [sn0(£j)  sinc  (j  E)  + No5,C)]  L‘j 


Jka8d«, 


' Note  that  the  linear  array  allows  introducing  tne  sensor  noise  in  term  of  an  effective  wave  number  k 
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Figure  5-23.  Distribution  of  noise  power  versus  ka,  the  effective  wave  number  across 
the  array. 


where  (J(cj)  is  a normalizing  constant  for  unity  gain  at  k=k  j-  The  solution  of  this  integral 
equation  can  be  expressed  in  terms  of  the  spatial  eigenvalues  and  eigenfunctions  of  the  homo- 
geneous integral  equation  associated  with  the  temporal  frequency-spatial  covariance  of  the 
isotropic  noise.  For  this  particular  kernel,  the  prolate  spheriodal  wave  functions  specify  the 
basis  for  the  solutions.  When  one  maps  the  parameters  of  the  prolate  spheroidal  wave  functions 
to  the  array  design  problem  considered  here,  there  are  2L^  + 1 significant  eigenvalues.  As  an 
indication  of  this  behavior.  Figure  5-24  illustrates  the  six  most  significant  eigenvalues  as  a 
function  of  array  length  L.  For  2L^  > 1 the  eigenfunctions  resemble  sinusoids  where  the 
spatial  frequency  separations  ts-j-  . 

If  we  denote  the  eigenvalues  as  7n(u>)  and  the  eigenfunctions  as  ll>n(co:C),  the  optimum 
array  beam  pattern  is,  using  Eqs.  4.1 2d  and  4.15, 

g0(«.k|  kT)  = <?  (w|  kT) . 


-Jk'aaSl  fhl2 

e 0n(°J-sl)d®l 

•'-L/2 


f1#  * JkT‘^aS2  //  \ 

1/  *n<w:*2>e  dS2/(Vt',>  + No) 

•'-L/2  / 


(5.56a; 
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where  the  performance  as  encompassed  by  the  output  noise  power  is 


o (co  | k-j-) : 


t / 

, n=i  -L/2 


,, ,,  i -%-5aS  1 2 ’ 

L'2  le  «n(tu:e)dS 


(5  56b) 


When  the  array  length  exceeds  unity,  the  eigenfunctions  approach  a sinusoid 
harmonically  related  at  a fundamental  2it/L.  Switching  to  a complex  exponential  represen- 
tation for  the  sines  and  cosines  we  then  have 


a (w  Kj) 2 


rn(u)  + N0 


No  / X Sno(C0)\ 

T KlTT 

No/  X V*>\ 

TKin ]’ 


4[« 


1 + LXPka(“  kTaf 


(5.57) 


|kT‘§a|  <2irA 

lkT'2al  “ 2tr/X  (endfire) 


* o 
L ' 


l-T'^al  b®yond  region  of  propagating  targets 


This  represents  the  approximate  performance  of  both  the  conventional  and  the  optimum 
arrays.  It  is  only  in  regions  where  the  isotropic  noise  is  dominant  that  the  optimum  array  can 
achieve  somewhat  superior  performance  by  displacing  the  main  lobe  into  a region  with 
k-aa  > ??  , where  only  the  white  noise  is  present.  This  is  an  example  of  how  a superdirective 

A 

beam  is  formed. 

In  Figure  5-25  we  have  illustrated  the  performance  found  by  solving  Eq.  5.55.  We  have 
also  indicated  the  results  indicated  by  the  approach  discussed  above  for  > 1 . The  arrays  show 
a gradual  decrease  noise  power  output  at  endfire.  This  decrease  can  be  interpreted  as  putting 
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the  array  wave  number  response  in  a region  where  the  frequency  wave  number  spectrum  of  the 
noise  is  lower  The  conventional  array  performs  almost  exactly  as  predicted  by  the  large  array 
approximation,  while  the  optimum  array  displays  an  improved  performance  for  the  high 
isotropic  noise  case.  This  can  be  attributed  to  superdirective  effects. 

The  presence  of  high  end  fire  gains  has  been  observed  by  several  people,  while  the  dis- 
cussion by  Vanderkulk  [25]  is  best  known,  the  first  observation  was  by  Schellkenoff  [26] 

A report  by  Schulteiss  [27]  is  notable  because  it  has  several  sets  of  figures  Its  equivalent  in 
signal  design  for  temporal  signals  apparently  has  not  been  discussed. 

If  we  examine  Figures  5-26  and  5-27,  the  optimum  beam  patterns  for  operating  in  an 
isotropic  noise  environment,  we  observe  that  when  the  target  is  broadside  the  beam  pattern  has 
a conventional  shape.  In  the  case  of  very  high  isotropic  noise,  the  beam  is  slightly  displaced  into 
a region  of  nonpropagating  or  sensor  noise  only.  Again  we  observe  that  optimum  processing 
offers  substantial  improvement  when  one  can  place  part  of  the  beam  pattern  in  a wave  number 
region  where  there  is  a lower  noise  level  and  yet  satisfy  the  constraint  of  having  a unity  response 
in  the  target  direction. 

We  can  interpret  high  endfire  gain  for  isotropic  noise  directly  in  terms  of  three-dimensional 
frequency  wave  number  concepts.  This  interpretation  also  predicts  the  occurrence  of  similar 
behavior  for  other  types  of  noise  fields.  In  Figure  5-28  we  have  indicated  the  frequency  »vave 
number  function  of  isotropic  noise  as  distributed  on  a sphere  of  radius  2ir fk.  According  to  our 
intuitive  wave  number  analysis,  there  is  a strip  of  width  2 X 2vJL  representing  the  main  lobe  of 
the  beam  in  the  target  direction.  For  broadside  beams  this  strip  intersects  a ring-shaped  region 
whose  approximate  area  is  ^ X lit  X ?2,  as  indicated  jjy  target  location  on  Figure  5-28  At 
endfire  the  region  has  a shape  of  a cap  whose  area  is  ir  1 2?  ^ Q j" , as  indicated  by  location  B. 


rlil 
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The  intersected  area  of  the  endfire  region  is  half  that  of  the  broadside  region.  The  gain  can  be 
additionally  enhanced  by  displacing  the  lobes  outside  the  sphere  and  increasing  its  gain.  If  the 
noise  does  not  have  the  uniform  distribution  of  isotropic  noise,  the  intersected  regions  must  be 
weighted  appropriately.  For  example,  if  there  is  a high  noise  component  normal  to  the  array,  the 
endfire  performance  will  be  enhanced.  Such  would  be  the  case  for  a vertical  array  in  a noise  field 
such  as  that  shown  in  Figure  2-8.  Conversely,  the  high  layer  noise  of  Figure  2-4  would  reduce 
endfire  effects  of  the  array  performance.  Some  typical  results  for  noise  fields  with  nonuniform 
distributions  are  indicated  in  Figures  5-29  through  5-32. 
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Figure  5-31 . Output  noise  power  for  a 3\  vertical  linear  array  in  high-surface  noise 
environment  (see  Figure  2-4). 


5,0  SUPERDIRECTIVE  LINEAR  ARRAYS 


Superdirective  arrays  have  a long  and  rather  interesting  history  The  work  on  them  lias 
been  motivated  by  the  apparent  ability  to  achieve  very  small  beamwidths  with  finite  length 
apertures.  One  of  the  first  to  recognize  their  existence  and  properties  was  Schelkunoff  [ 26] 
in  his  algebraic  theory  of  linear  array's  They  have  also  been  extensively  discussed  in  electro- 
magnetics, where  the  problems  inherent  in  their  realization  have  been  investigated  Several 
synthesis  procedures  are  related  to  superdirective  arrays,  in  particular  those  of  Pritchard  [ 28] 
and  some  aspects  of  Dolpli-Tchebychev  arrays  They  are  also  intimately  related  to  singular  detec- 
tion theory  which  in  itself  has  seen  a long  period  of  development  [3]  Quantitative  results  are 
hard  to  come  by,  however,  the  general  principles  about  superdirectivity  and  its  relation  to  singular 
detection  can  be  understood  quite  intuitively 

For  purposes  of  illustration  we  use  a linear  array,  although  there  exist  examples  where 
other  types  of  arrays,  e g.,  circular,  where  similar  approaches  have  been  used  If  we  have  an  array 
of  length  L o-iented  with  an  array  tangent  aa,  its  frequency  wave  number  response,  in  principle, 
extends  to  all  spatial  frequencies  If  one  considers  the  response  to  plane  wave  signals  in  a medium 
with  finite  propagation  velocity,  then  the  incident  wave  numbers,  or  spatial  frequencies,  are 
cj  lit 

limited  in  magnitude  to  |k|  < — = — The  response  outside  this  range  is  theoretically  arbitrary, 
c X 

For  a superdirective  array,  the  wave-number  response  m this  region  is  manipulated  such  that 
within  the  range  |k|  < — it  has  a very  narrow  beamwidth 

We  have  pointed  out  the  difficulty  with  this  in  a practical  system,  the  narrow  beamwidths 

2tr 

are  created  by  generating  extremely  large  sidelobes  in  the  region  |k|  > — . If  one  investigates  the 
beam  pattern,  or  wave  number  response,  by  plotting  it  versus  0,  these  sidelobes  never  appear. 

Sensor  noise,  however,  does  not  have  a propagating  structure  but  has  an  equivalent  wideband 
wave  number  spectrum  which  can  extend  well  beyond  the  limit  on  propagating  limit  of  |k|  < 2 tt/\ 
Consequently,  it  enters  through  in  these  large  sidelobes  as  we  can  see  from  the  second  term  in 
Eq  5.51. 

We  can  indicate  some  superdirective  effects  for  a half  wave  length  array  by  designing  a beam 
which  has  half  the  width  of  a conventional  pattern  We  use  the  procedure  of  Section  5.2  and  let 
N0  *♦  0,  or  0 -*  <®.  This  yields 

q(<o.k|kT)=  C Jsmc^(k-kT)-aa  yj 

4 (S,"C  [f--T --a4)  • 3a  l]  + sinc  jj^-"-T  + -a  ' 3a  j])]’ 
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rare  the  output  noise  level  of  the  superdirective  array  to  that  of  a conventional 
the  ratio  of  the  ratio  is 

.2 
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The  superdirective  array  is  much  more  sensitive  to  sensor  noise 

The  optimum  design  of  superdirective  arrays  is  not,  surprisingly,  closely  related  to  the 
optimum  processing  methods  that  we  have  been  considering  Also,  this  illustrates  the  connection 
between  superdirectivity  and  problems  m singular  detection  theory  An  ideal  array  would  have 
a unity  response.1 


g j(wka)  = 


*<ka-Ta>-  l*al  <*/* 

9S(“  ka>-  lkal  >2b/X 


For  an  array  of  length  L we  want  the  transform  of  this  response  to  have  the  form 


Ge(cu  8)  = 


Gg(w:C)  = e 


finite  , 181 « L/2 


181  > L/2 


( 

+ ,u 


-jka8  dka 


Gg(w:8)  = e Ta  + Gg(cu:8)  * jj5(8)  - -j  sine 


This  yields 


j J d82  sine  Jjy  (£i  " g8(w  *2* = e 


L/2  < 8,  < L/2  (5.62) 


K 


This  is  the  integral  equation  which  we  would  obtain  if  we  posed  the  optimum  detection  prob- 
lem of  detecting  a plane  wave  signal  in  the  presence  of  isotropic  noise  alone,  with  no  white  noise 
present  Such  gn  integral  equation  is  a Fredholm  equation  of  the  first  kind,  and  the  general 
solution  has  singularity  functions  appearing  at  the  end  points  of  the  observation  interval  for 
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‘As  noted  earlier  this  value  of  unity  near  the  target  direction  must  yield  a nonzero  finite  measure  when  integrated  in  a 
region  "near"  it 
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temporal  processes  or  at  the  aperture  extremities  for  spatial  processes  One  can  specify  a forma! 
solution  in  terms  of  the  solutions  of  the  homogeneous  equation— the  prolate  sphoidial  wave 
functions;  however  the  resulting  series  does  not  converge,  in  general,  to  well-defined  functions 

The  essential  points  to  be  made  regarding  superdirective  arrays  are 

1 . They  are  primarily  a mathematical  issue  since  the  effects  of  inter-element  coupling 
and  the  sensitive  adjustments  required  cannot  be  controlled  precisely  enough  to  implement 
them  practically 

2.  Their  sensitivity  is  closely  coupled  to  singular  detection  and  is  of  interest  as  a 
limiting  situation  when  the  white  noise  in  a system  becomes  small 

3.  They  are  only  effective  across  narrow  frequency  bands  since  the  aperture  phasing  is 
coupled  to  the  wavelength  specification. 

4 Their  response  should  be  determined  for  all  wave  numbers,  not  simply  as  a function 
of  the  incident  beanng  angle 

1 6.  OPTIMUM  ARRAY  PROCESSING  FOR  DIRECTIONAL  NOISE  FIELDS 


Noise  fields  that  consist  of  a finite  number  of  directional  sources  plus  a white  component 
can  be  analysed  exactly.  The  final  results  require  knowledge  of  sincj^(k),  the  conventional  wave 
number  response  of  the  array,  the  correlation,  if  any,  among  the  noise  sources,  and  a matrix 
inversion  Since  we  have  observed  in  Section  3.6  that  directional  noise  field  has  a finite  number 
of  nonzero  eigenvalues,  we  should  be  able  to  specify  the  solution  in  terms  of  a set  on 
linear  equations 

The  noise  field  has  the  temporal  frequency  spatial  covariance  function 
Xr^>  x?"'  j-i’i  _J!Sj 

sn <“:£P  = 2L/ 2-te  S<j(W) * + N0  6^  (z-f),  (6.1) 

i=l  j=l 

where  Sjj(cj)  represents  the  correlation  among  the  various  directional  components  Note  that  if 
we  require  the  directional  noise  to  form  a homogeneous  field,  then  homogeneity  requires  that 
plane  waves  with  different  wave  numbers  be  uncorrelated.  For  convenience  in  the  subsequent 
discussion,  we  define  [S(co)J  to  be  the  matrix  whose  elements  are  S,j(co).  in  addition,  we  define 
the  propagation  vector  e n to  be 
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where  the  kj  are  the  wave  numbers  for  the  various  noise  components.  Consequently,  we  can 
express  Eq.  6.2  as 


'i! in-^ 


JiSn'i 


Sn(cj:z,r)  = e (S(m)1  e +N0«j2(z-f) 


(6.3) 


To  determine  the  wave  number  response  and  the  output  power,  we  must  find  the  inverse 
kernel  across  the  aperture  fl.  We  appeal  to  a functional  analog  of  the  following  matrix 
inverions  lemma 


[uTv  + A]->  = A"*  - A'^u^  [£+vA"*u^)‘*  v A'* , 


where  A is  N X N,  and  u and  v are  M X N.  We  identify  A as  N0  6jj  (z-J),  while  u and  v are 


replaced  by,  respectively,  n and  (S(oi)]e  -.  The  spaces  of  concern  are  S2  X J2  and 


n X Rjq.  We  have  for  the  inverse  kernel 


Qn(wz|)  = — (z-£) 

”o 


-iV2f 


JL  " [In-a„!SM1  ( ,]  ISO 
N0  L N0  } N, 


[S(w)l  ‘iiSn’i 


(6.4) 


where  we  define 


1 sincn(krk2)  sincjjikj-kg) 

sincn(k2-!si)  1 

sincn^3'!ii) 


(6.5) 
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(The  term  [pn]  is  identical  to  the  term  (smcn(8)]  used  in  Eq  3 57c  ) This  matrix  reflects  the 
mutual  interference  of  the  various  noise  sources  as  observed  across  the  aperture 

The  output  noise  power  as  a function  of  the  target  dnection  follows  from  Eq.  4 17  as 


smcj^k-kj) 

PN<!s)  - sincj2(k-k2) 

_sincfi(k-kN) 

This  can  be  simplified  considerably;  the  final  expression  is 


<7*{M|kT)=  ^ [l  -TrWl  p(kT)  pf(kT)  (IN+[0J  [Pnll'1)]'1,;  (6.7 


Ajj[S(oj)] 


This  is  the  same  form  we  had  for  the  single  noise  source  case.  The  limiting  performance,  or 
lowest  noise  level,  is  given  by  N0/Ajj,  while  any  increase  in  the  noise  output  level  enters  in  the 
form  of 


2 VAn 

Co(oj|kT)=  -y— - ,, 


“ = Tr  ([01  PN%>  eJ,<hT>  (In  + [01  [Pnl)'1 ) ,. 
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For  a conventional  wave  vector  response,  we  obtain 


oj(cj|kT)  = ^ 1 + Tr  [|/5]  p(kT)  £T(kT)j| 


(6.9) 


We  can  observe  that  the  mutual  interference  is  introduced  via  the  term  [p^]  in  the  optimum 
beam  pattern  output,  whereas  it  is  not  included  in  the  conventional  processing 

The  optimum  beam  pattern  follows  directly  from  Eq.  6.4  using  Eq  4 1 6 This  yields 


g(w:k|kT) 


= ancn(k-kT)  - Tr  |t(5l£N(kT)  P^)  [iN  + [g] 

1 -Tr{(/5]  pN(kx)pT(^T)  [IN  + [0]  [pnr'| 


(6  10) 


We  have  a conventional  beam  pattern  and  a term  which  introduces  nulls  in  the  direction  of 
the  vanous  noise  components. 

We  can  completely  analyze  the  optimum  array  performance  for  a source  of  arbitrary 
geometry  operating  in  a directional  noise  field.  Since  we  have  introduced  most  of  the  intuitive 
concepts  in  our  discussion  of  linear  arrays,  particularly  null  placement,  we  need  not  pursue 
this  further.  The  major  remainmg  issue  is  the  effect  of  various  array  geometries 


7.  ARRAY  GEOMETRIES  IN  HIGHER  DIMENSIONS 


The  analysis  of  arrays  with  geometries  that  are  more  complex  than  linear  generally 
becomes  quite  complicated.  Since  our  primary  consideration  is  the  development  of  the  major 
concepts  that  are  involved  in  optimal  array  processing  techniques,  we  have  focused  our  attention 
upon  linear  arrays.  The  more  tedious  mathematical  issues  did  not  appear,  yet  we  were  able  to 
analyze  substantive  examples  illustrating  the  issues  under  consideration.  In  many  respects,  this 
is  deceptively  simple,  in  that  added  dimensions  introduce  unique  problems  and  possibilities 
which  have  no  direct  temporal  equivalent.  This  does  not  imply,  however,  that  the  introduction 
of  more  complex  geometries  changes  the  essential  concepts  we  have  developed  Higher 
dimensional  problems  are  sufficiently  complicated  such  that  any  insights  that  can  be  drawn 
from  simpler  problems  are  valuable 

In  a discussion  of  optimum  processing  with  more  general  array  geometries,  we  cannot 
pursue  as  extensive  an  analysis  as  we  did  with  linear  arrays.  There  are  simply  too  many 
combinations  of  array  geometries  and  noise  fields  which  rapidly  reach  a point  of  diminishing 
return  once  a few  representative  examples  have  been  introduced  We  discuss  some  of  the  more 
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common  geometries  and  then  concentrate  upon  some  representative  examples  of  optimal 
processing  and  the  limitations  imposed  upon  it  by  vanous  noise  fields 

To  keep  issues  in  perspective  we  point  out  that  the  only  noise  fields  which  can 
reasonably  be  analyzed  are  directional  ones,  and  even  these  become  quite  complex  for  any 
more  than  two  sources  Exact  analysis  of  fields  with  noise  sources  distributed  in  spatial  wave 
number  is  impossible  even  for  such  a simple  structure  as  isotropic  noise  The  more  complex 
geometries  do  not  change  these  results.  We  have  already  done  a complete  analysis  for  the 
directional  noise  field  and  the  only  new  issue  is  computing  the  sincj^Ck)  function  The  exact 
analysis  for  distributed  noise  sources  remains  intractable,  and  we  must  still  resort  to  approxi- 
mate and/or  asymptotic  analysis 

In  the  analysis  of  linear  arrays  we  found  that  the  important  wave  number  was  k*aa— the 
component  of  the  wave  number  as  projected  upon  the  array  tangent  We  are  similarly  concerned 
here  with  projections  on  a set  of  lines  for  crossed  arrays,  > a '5  for  planar  arrays.  Just  as  for 
linear  arrays,  these  projected  wave  numbers  can  be  quit  . calculating  wave  number 

responses 

When  dealing  with  white  noise  and  frequency  wave  number  issues  one  must  be  careful  As 
before,  we  regard  this  class  of  noise  as  nonpropagating  in  origin  and  independent  from  point  to 
point  on  the  array  This  leads  to  a wave  number  representation  which  is  flat  where  the  array  is 
capable  of  filtering  and  no  definition  for  the  entire  wave  number  space  is  implied  If  the  completely 
flat  response  for  all  wave  number  space  is  assumed,  we  are  led  to  infinite  noise  power  propagating 
through  our  spatial  filters.  This  is  not  a very  satisfactory  model  of  a real  array  Essentially,  we  want 
the  Parseval  relation 


N°/'/ l9(w  '|2  ^ = N°  / G|(w  z)|2  dz 


In  general  this  is  no]  true,  and  only  for  those  special  cases  in  which  we  confine  our  integration  to 
projected  wave  number  components  can  we  wnte  relations  of  this  form. 

There  are  three  situations  of  general  interest  in  arrays  with  spatial  extent  of  more  than  one 
direction  • the  target  and  any  directional  noise  sources  are  near  broadside,  the  target  is  broadside 
and  significant  components  of  the  noise  field  propagate  across  the  array,  both  the  target  and  noise 
propagate  across  the  array.  Each  introduces  its  particular  set  of  considerations.  Several  typical 
cases  are  illustrated  in  Figure  7-1 . 
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Figure  7-1.  Typical  array,  target  and  noise  geometries. 

7.1  PLANAR  ARRAYS 

Planar  arrays  are  the  simplest  two-dimensional  army  structures  Our  analysis  of  planar 
arrays  is  essentially  the  same  as  for  linear  arrays  Figure  7-2  illustrates  the  array  geometry  under 
consideration,  for  convenience,  we  orient  the  array  with  a unit  normal  vector  a^  and  tangent 
vectors  aj  and  g2  (these  are  not  necessarily  orthogonal). 

We  can  specify  the  array  weighting  pattern  as  a two-dimensional  function 

L,  L2 

G(<o:z)  = GfuCja,,  C2a2),  jS,|  <-J-,|«2|<_  (7.1) 

We  can  attach  all  the  common  weightings  or  shadings  as  we  did  with  the  linear  array 
for  each  coordinate  8 1 and  82,  or  we  can  specify  more  complex  shadings  with  a function 
which  need  not  factor.  We  consider  some  examples.  For  a target  at  Kx*2.1  = h-y-ai  = 0 with 
conventional  shading,  G(w  £|  ,82)  = 1 /Lj  L2,  or 
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Figure  7-2.  Lj  x L2  planar  array  with  orientation  aj  x a2  oan. 


As  we  have  seen,  sincjj(k)  is  important  in  our  discussion  of  directional  noise  analyses  and 
its  form  is  sketched  in  Figure  7-3.  Figure  7-4  is  the  beam  pattern  when  the  array  response  is 
presented  in  dimensionless  units.  If  we  examine  the  array  pattern  in  Figure  7-3,  we  see  that  it  has 
aresolution  of  2-(2ir/L|)  in  the  aj  coordinate  direction  and  2-(27r/L2)  m thea2  directior  The 
resolution  cell  has  an  area  (2jt/Lj  ) X ( 2tt/ l-2 ) m wave  number  space.  The  distribution  of  this  area 
is  determined  by  the  orientation  vectors  §j  and  a2 

A triangular  shading  in  each  direction  produces  a beam  pattern  of  the  form 


2/  LA  7/  l2\ 

g(w  k)  = sinczUc-a|  — J sine'1 1 k'a2  —I 
We  can  also  split  the  array  such  that 

g(wk)=  |Fj[(k-Ak)-a,]  -F,  Kk+AlO-a,]} 

* {f2  [(k-Ak)-a2J  -F2(-(k+Ak)-a2l|, 
where  F|  and  F2  are  the  beam  patterns  on  the  array  halves 


(7.3) 


(7.4) 
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This  beam  pattern  typically  has  a saddle  point  structure  near  the  origin  in  wave  number 
space  where  the  cross  defining  the  critical  point  locus  of  the  saddle  is  along  the  vectors  aj  and 
a2.  We  see  that  in  this  situation  the  generation  of  single  beam  is  not  a desirable  processing  method 
We  need  to  generate  two  beam  outputs,  one  for  each  coordinate  direction  in  order  to  obtain  our 
desired  response  for  nulling  at  a particular  point 

The  analysis  for  optimal  processing  in  the  presence  of  directional  noise  is  substantially  the 
same  as  for  a linear  array.  The  performance  may  be  taken  directly  from  Eq  6.8  as 

2 N0  1 + (L,L2/N0)S  (co) 

°0(“|-t)=  l^t; 2 — (75) 

1 2 1 + (L  j L2/N0)  SnQ( u)  ( 1 -smcn  (k-j-kn)) 

The  beam  pattern  also  follows  directly  from  Eq.  6.10 


g(u  k |kq-l  = 


sincj^fk-kf)  -smcjj(kfkn)  Sincjj(k-kn) 


1 - sincfi  (k-j-kn) 


(7.6) 


In  Figures  7-5  through  7-10  we  indicate  the  performance  for  some  typical  noise 
geometries  We  observe  the  same  effects  that  entered  in  our  study  of  linear  arrays  Note  the  high 
sidelobe  levels  in  the  conventional  versus  the  optimum  processing  The  results  when  two  direc- 
tional sources  are  in  the  noise  field  are  parallel  to  those  for  linear  arrays. 

For  analyzing  isotropic  noise  we  encounter  the  same  type  of  difficulties  as  in  the  linear 
case  if  we  pursue  an  ^xact  analysis.  We  can,  however,  pursue  as  approximate  analysis.  First,  we 
need  to  find  the  asymptotic  distribution  of  the  eigenvalues  for  the  isotropic  noise  kernel  expanded 
over  the  aperture.  The  eigenfunctions  will  be  spatial  sinusoids.  Asymptotically,  the  eigenvalues 
approach  a skewed  version  of  the  two-dimensional  wave  number  transform  as  the  array  extent 
increases.  (We  assume  that  the  array  is  located  in  the  (zx,  zy)  plane  ) We  have  from  our  previous 
results 
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Computing  the  noise  output  as  we  go  from  broadside  to  endfire  is  more  difficult  for  there  is  not 
a monotomc  increase  as  is  the  case  in  the  linear  arrays  If  L > 1 and  L2\>  1 the  eigenfunctions 
will  approach  spatial  sinusoids  with  an  eigenfunction  distribution  given  by  the  above  two  dimen- 
sional wave  number  function  An  approximation  performance  is 


°o(“l-T)=E^  (N0  + P<^kTa) 


(7  8) 


_No 

L)L2 


No 

LjL  2 


s„  («) 


{m 


-1/2 


i 2tt 

• N < T 


>r 


Tins  diverges  as  kT  approaches  endfire  or  2ir/\,  so  it  is  hard  to  predict  the  endfire  performance 
It  is  easy  to  observe  that  our  approximation  can  be  refined  with  two  considerations  If  we  inte- 
grate over  the  finite  beamwidth,  we  average  over  the  integrable  singu  anty  in  the  wave  number 
function  In  addition,  at  endfire  half  of  the  beam  response  is  in  the  white  noise  region  which  has  a 
lower  noise  level  Consequently,  we  have  that  the  output  noise  increases  until  we  are  beyond 


. - 2ir  2 ir  2r r / I \ 


at  winch  point  it  levels  off  and  then  begins  to  decrease  The  optimum  processing  can  increase 
these  effects  somewhat  by  appropriate  showing  of  the  beam  Some  representative  performances 
are  illustrated  in  Figures  7-1 1 and  7-12. 

One  can  introduce  various  noise  field  geometries  representing  surface,  bottom  or  layer 
noise.  The  results  of  our  discussion  on  noise  representation  indicate  how  the  two  dimensional 
wave  number  function  can  be  found,  from  which  the  approximate  performance  can  be 
determined  Alternatively,  we  can  pursue  an  anal>  sis  in  three  dimensional  wave  number 
space  with  the  noise  field  distributed  on  a sphere  as  we  did  in  our  discussion  for  linear  arrays. 
We  are  led  to  the  same  results  in  either  method. 


7.2  CROSSED  ARRAYS 


The  importance  of  crossed  arrays  lies  in  their  capabilities  of  providing  spatial  information, 
or  resolution,  along  two  coordinates  with  fewer  sensors  than  a planar  array  While  in  a planar  array 
the  number  of  sensors  increases  as  L j L2  of  the  crossed  array  sensors  increases  as  Lj  + L2  In  most, 
but  not  all,  cases  crossed  arrays  perform  virtually  as  well  as  their  planar  counterparts  when  their 
outputs  are  processed  correctly.  This  is  very  significant  in  reducing  computational  demands 

One  of  the  most  important  facets  of  the  analysis  of  crossed  arrays  is  their  tutorial  value 
With  a minimum  of  complexity,  they  introduce  several  significant  concepts  in  the  processing  of 
array  data.  In  particular,  the  importance  of  frequency  wave  number  concepts  as  a means  of 
obtaining  approximate  expressions  for  the  system  performance  For  our  analysis  we  consider  two 
noise  fields  — directional  and  isotropic. 

We  consider  an  array  structure  as  illustrated  in  Figure  7-13.  The  class  of  beam  patterns 
which  we  can  generate  is  given  by 


,L,/2 

g(u)  k)  = / G|(<j: 

4,/2 


jk-ai«i  fl 
■«l)c  dC,  +J 


jk-a2f2 

<182 


g|(«.k*aj)+  g2(co:k-a2) 


Figure  7-13.  Crossed  array -legs  of  length  L|,L.2>  oriented  in  direction  a j,  ^ crossed  at 
center  of  legs. 
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This  is  a much  more  restrictive  form  than  that  which  we  had  for  the  planar  array  We  can- 
not generate  resolution  cells,  1 e , regions  where  there  is  a main  beam  concentration  and 
outside  of  which  the  beam  is  approximately  zero  With  a crossed  array  we  can  generate  only 
combinations  of  resolution  strips  In  Figures  7-14  and  7-15  we  have  plotted  the  conventional 
beam  pattern  for  a typical  crossed  array  configuration  This  beam  pattern  is  given  by 


g(to  k|kT)  = 


l1+l2 


Lf 

l2 

h 

sine 

(k-kT)-a,  — 

+ sine 

l,+l2 

(k-kT)-a2  — 

(7  10) 


Conventional  beamforming  can  introduce  ambiguity  effects  if  a noise  source  appeared  in  one  of 

the  “sidelobe  strips,”  one  should  weigh  the  use  of  the  other  line  of  the  cross  to  a degree  greater 

than  that  of  the  conventional  weighting  of  Lj/Lj+L-)  Optimal  processing  achieves  a substantial 

amount  of  its  improved  performance  and  sidelobe  suppression  by  introducing  this  implicitly 

rather  than  by  using  esoteric  beam  patterns 

We  consider  the  analysis  of  a directional  noise  source  of  level  S„  (to).  From  our 

"o 

previous  example,  we  can  specify  the  performance  and  optimal  beam  pattern  by  inspection 
For  the  performance  we  have 
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where 


a (co  |kj) 
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An  = Li + l2 

L,  / LA  L2  / 

smc«®=  r^r2  s,ncf2i  T/+r^r;  smcfa-2 

This  is  plotted  in  Figures  7-16  through  7-18  At  first  inspection,  it  may  appear  that  the  above 
does  not  perform  as  well  as  a simple  linear  array  when  the  noise  appears  in  one  of  the  sidelobe 
strips  However,  a comparison  of  this  situation  with  the  corresponding  linear  array  yields 
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The  largest  increase  occurs  when  the  noise  appears  in  a null  of  the  crossed  array  length  with, 
i e.,  it  is  the  main  beam  of  one  length  and  the  null  of  the  other, 
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It  is  useful  to  demonstrate  how  the  array  obtains  its  increased  performance.  The  optimum 
beam  pattern  is  computed  straightforward,  and  the  next  result  is 
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We  consider  a situation  (kn-.kx)'a2  — 0 and  (kn-kx)-a]  > ir  The  resulting  beam  pattern 
generated  is  illustrated  m Figure  7-19  It  simply  has  a null  in  the  crossed  array  beam  pattern 
In  tile  special  case  of  a dominant  directional  noise  source  and  L|  =L->  we  have 


(L,+L2)Sn  (co) 


» 1 


1 

sincjj  (krkn)  =*  - 

i and 

| 

\ G0(co:k|kx)  =■  j (sincfi(k-kT)  - i sincn(k-kn))  (7.16) 


This  places  a null  at  k=kn  while  maintaining  unity  response  in  the  target  direction  A natural 
question  which  arises  is  why  doesn’t  the  array  simply  ignore  that  section  along  which  there 
is  no  resolvabihty,  i.e  , assign  a weighting  of  zero  along  that  cross  If  this  were  done,  we 
i would  have  an  array  output  noise  of  approximately  NQ/L  while  if  the  optimum  response  is 

used,  the  output  noise  is  2 dB  less,  or  4/3  N0/2Lj.  Both  weightings  suppress  the  highly 
directional  source  by  nulling  it,  however,  in  the  optimum  array  a smaller  amount  of  white 
noise  enters  through  its  reduced  sidelobe  level  Only  if  the  noise  field  has  a strip  of  sources  with 
(kq— kn)'»2  — 0 does  the  beam  pattern  approach  that  of  a linear  array  A representative  beam 
pattern  for  a noise  field  approximating  this  class  is  illustrated  in  Figure  7-20  Note  that  theie  is 
effectively  no  resolution  along  ky  direction 

For  crossed  arrays  in  directional  noise  fields  the  optimum  patterns  place  nulls  at  the 
wave  numbers  of  the  directional  noise  sources,  while  simultaneously  minimizing  the  mean 
square  amplitude  of  sidelobes  While  the  conventional  array  is  designed  independent  of  the  noise 
field,  the  optimum  patterns  indicate  an  intuitively  appealing  design  to  exploit  the  structure  of 
the  noise  fields  encountered  These  optimum  patterns  are  not  too  complicated  and  it  is  only 
when  the  noise  and  targets  are  close  together  that  they  become  complicated,  introducing  the 
features  we  encountered  earlier  in  the  discussion  of  linear  arrays 
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The  design  of  crossed  arrays  for  isotropic  noise  cannot  be  pursued  analytically  without 
considerable  amount  of  mathematical  detail  The  crossed  geometry  still  reduces  to  the  prolate 
spheroidal  wave  functions  that  we  encountered  in  linear  arrays  We  can,  however,  pursue  an 
approximate  analysis  which  encompasses  the  basic  attributes  of  the  processing 

To  analyze  the  performance  of  crossed  arrays  m noise  fields  which  have  a smoothly 
varying  wave  number  function  for  Ik  l > 27t/X,  we  must  isolate  a set  of  temporal  random  processes 
which  characterize  the  array  outputs  We  can  do  this  quite  easily  using  wave  number  > i epts 
Each  leg  of  the  array  forms  a set  of  orthogonal  beam  outputs  Xm(w)  and  Yn(w)  at  intervals  of 
27r/L,  m wave  number  space,  as  illustrated  in  Figure  7-21  For  convenience,  we  subsequently 
assume  the  array  crosses  are  perpendicular  The  spectra  of  each  of  the  outputs  Xm(w)  and  Yn(w) 
is  given  approximately  by 


with 


Sx  (w>s 
xm 


1 C ’ ' ' B 2ir  dk,  N, 

E7#  Pn(c  m-a|+k2a2)_^r 


T—  (7.17a) 


Sv  (u)  = 


1 / ' ' ' / , 2w  \ dk,  N0 

r2J  + (717b> 


[|m|]  < L,/X  = L,/\ 
[|n|]  < L2A  = L2/X 

a, -a2  = 0 
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Figure  7-21 . Wave  number  analysis  for  crossed  array  performance. 

where  Pn(u  k)  is  the  two-dimensional  wavenumber  function  Consequently,  we  have  that 

there  are  2(Li  + L2  ) random  processes  The  random  processes  xm  (t)  and  x (t)  are  un- 
' A A 12 

correlated  since  the  beams  are  orthogonal;  however,  xn(t)  and  ym(t)  are  correlated  since 
the  beams  overlap  at  the  crossing  point  We  have  for  their  cross  spectral  density 


n , , 1 J 2r  2ir  \ 

sW")=^r2  plf:wr;zi+nr2y 


so  that  as  the  array  length  L]  and  L2  increase,  the  effect  of  the  correlation  decreases 


Figure  7-2 1 Wave  number  filters  for  crossed  array  performance 
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where  a and  p are  adjusted  to  weight  each  array  so  as  to  minimize  total  noise  output  of  the 
beam  We  have 


(7.20) 

Optimizing  a and  p yields 


aLl 

«L]+0L2 


S*  (w)-S. 


xmy, 


(w) 


Sx  (w)  + Sv  (w)-2S 


xmyn 


(w) 


PL2  SXn«u)-SX[nyn«o) 

«L[+/SL2  sx  (u)  + Sv  (u)  - 2SX  v (w) 
m yn  xmyn 


This  yields  a performance  of 


Sx^^-S^co) 

a (wkt)  =* 

o -1  Sx  (co)  + S (co)  - 2SX  v (co) 


(7.21a) 


(7  21b) 


(7.22) 


where  the  various  spectra  are  computed  via  Eq.  7. 1 7 and  7 1 8 


V V -■  -*  VVV  VV  V .*.■  V „■  .. V.  -r. 


As  an  example,  we  consider  the  special  case  of  three-dimensional  isotropic  noise  of 
Example  2.2a  with  a crossed  array.  We  have 


sx  (“) ' 
*i 


(|  vw)+N°) /L' 

No lh< 


'■Kl.tl 


U,  >|-| 


(7  23a) 


Sv(oi)  = 


Substituting  this  into  Equation  7 22  we  obtain  for  the  special  case  of  Lj=L2  with  P » 1 


2 

O (tJ  ky)  = 


% (l  + p0  \ 

2L|  \ir/(Li/X)/ 


(7.24) 


where 


ASn  (co) 
no 


2N„ 
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The  improvement  o»er  a simple  linear  array  is 


We  have  assumed  that  Lj/A>  1,  otherwise,  we  would  have  need  to  consider  superdirective 
effects  immediately.  In  approximating  the  smoothness  of  P(te  k)  we  must  have 


or 


-^-<1 

(irL/4A) 


If  this  is  not  satisfied,  then  the  improvement  becomes  1/2.  Finally  if 


then  we  should  be  careful  of  superdirective  effects,  in  summary,  Figure  7-22  represents  an  example 

2ir  2 

of  a superdirective  beam  pattern  when  kj- = , 0, 0 , while  plots  for  0o(oi|k-p)  may  be  found  in 

Figure  7-23.  Observe  that  in  Figure  7-23  at  0p  = 0°  it  is  easier  to  place  beam  pattern  volume  outside 
|k|  > 2ir/A  than  it  is  at  = 45° , therefore  the  endfire  affects  are  more  significant 

We  can  also  use  the  above  analysis  to  derive  approximate  expressions  for  the  spatial 
eigenvalues  and  eigenfunctions  of  the  crossed  array  when  operating  in  a noise  field  whose 
frequency  wave  number  function  vanes  slowly  over  the  resolution  width  maximum  (L/A)-' 
of  the  array.  The  results  are  indicative  of  what  is  obtained  in  the  analysis  of  higher  dimensional 
geometries.  We  emphasize  that  our  analysis  is  approximate.  It  combines  a classical  temporal 
analysis  with  the  features  of  the  crossed  array.  To  determine  the  eigenvalues  and  eigenfunctions 
we  observe  that  at  the  frequency  <o  the  random  variables  Xj(cj)  and  yj(io)  lead  to  a spectral 
covanance  matrix  of  the  form 
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-Sx.nM  0 


0 


[S(co)]  = 


Sx,y1(“) 


sx,y/“) 


Sy-M(w)  0 


0 V"* 


where  N = Lj/X,  M = L2/X  and  the  elements  Sx.(w)  are  computed  as  indicated  in  our  previous 


discussion  The  generation  of  a set  of  random  variables  which  have  a diagonal  covariance  matrix 
is  a classical  orthogonalization  problem  in  matrix  algebra  [31 . One  needs  to  solve  for  the  eigen- 
values and  eigenvectors  of  the  homogeneous  equation 


[[S(<4]-  Xk(<j)  l]  4>k(cj)  = 0 


The  eigenvalues  of  the  process  at  frequency  cj  are  given  by  the  set  jx^Cco)  ^ (N+M+ 1 ) j 


v / 

The  spatial  eigenfunctions  can  be  expressed  m terms  of  the  elements  of  the  eigenvectors.  If 
we  separate  them  such  that  the  first  2N  + 1 components  are  associated  with  the  x or  Lj  leg, 
while  the  remaining  2M  + 1 are  with  the  y or  L2  leg,  we  have 


ck, 


K-mM 


Ck.  M<"> 
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and  the  eigenfunction  ^(co  z)  associated  with  Aj,(oj)  is 


(aZ  ck,,(“)e' 

1-N 


j(2tt/L|  )Cai  -z 


0k(co.z)  = 


M 


“IX, 

. i=-M 


(to)e 


j(2ir/L-,)ia-,-z 


|z  = Sa,  |C|  < L,/2 1 

(7  27b) 

|z  = %2  |«|<L2/2l 


where  a is  a normalization  constant  Note  that  the  basic  approach  is  to  isolate  a set  of 
random  variables  from  which  we  could  reconstruct  the  observed  signal  and  then  to  apply  an 
orthogonalization  procedure. 

We  have  assumed  that  our  array  is  crossed  at  the  center  of  each  leg  This  simplified 
many  of  our  results  since  there  was  no  relative  phasing  between  the  center  of  each  leg  We 
now  briefly  consider  the  situation  in  which  the  center  is  separated.  We  assume  a geometry  as 
illustrated  in  Figure  7-24.  The  class  of  beam  patterns  which  can  be  generated  is  given  by 


jk*Rj  jk‘R2 

g(oj:k)=  gj<<o:k-aj)  e + g2(<ok*a2)e 


(7  28a) 


where 


g1(w:k,)  = 


jk]S] 

G](<o  kj)  e d£| 


(7  28b) 


g2(w:k.2)  = 


jk2S2 

G2(w  2j)  e d^2 


(7.28c) 


with,  respectively,  Gj(co  E ( ) and  G2(cj  S2) as  the  weighting  patterns  for  the  first  and  second 
array  legs. 
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In  our  analysis  we  are  concerned  with  the  interference  relationship  between  the  two 
legs.  We  have 


2 2 2 
|g(cok)|  = lgi(tok-a[)|  + I g2<“  ira2)| 


■jk-(R,-R2) 
+ 2Re  e 


gj(w  k*aj)  g2  (w 


k-a2) 


(7.29) 


For  example,  a conventional  beam  pattern  with  only  a constant  magnitude  weighting  and  a 
phasing  directed  to  k=0, 1 e , k-p=0,  has  a power  response  of 


2 / L1  o L2\ 

|g(co:k)|  = a2sinc2  [k'aj  — + 02sinc2  |k-a2  — I 

/ LA  / l2\ 

+ 2ap  cos  (k‘(R|-R2))  sine  I k-a j — I sine  lk-a2  —J  (7  30) 


We  see  that  the  term 


►X 


v/, 

H 


cos[k-(R,R2)] 


\ 


can  alternatively  cause  the  responses  to  interfer  either  constructively  or  destructively  The  crossed 
array  with  separated  centers  produces  a larger  array  which  leads  to  the  narrower  beam  pattern  ndges 
Since  the  array  is  not  filled,  1 e , continuously  connected,  large  sidelobes  result  giving  the  indicated 
interference  pattern  This  is  a general  result  when  one  is  dealing  with  sparse  arrays  A typical 
situation  is  sketched  in  Figure  7-25  We  can  observe  the  general  effects  for  the  two  types  of  noise 
fields  that  we  have  considered  For  directional  noises  the  relevant  term  is  the  wave  number  differ- 
ence between  the  target  and  the  noise  ky  - kn  If  this  difference  vector  falls  within  one  of  the 
strips  but  not  near  the  origin,  we  have  essentially  the  same  problem  as  discussed  earlier  The  noise 
and  target  cannot  be  distinguished  along  that  particular  strip,  i e , it  is  within  a sidelobe  strip  so 
that  the  strip  of  the  second  resolves  the  two  Our  analysis  would  be  substantially  the  same  as  before 
If,  however,  the  difference  vector  is  near  the  origin,  we  may  be  able  to  use  the  array  separation  to 
our  advantage  to  achieve  an  enhanced  resolution  compared  to  that  of  the  centered  arrays,  if  the 
difference  vector  falls  at  a minimum  of  the  interference  pattern  the  separation  can  be  significant 
This  effect  is  clearly  most  pronounced  when  the  noise  difference  vector  is  parallel  to  the  array 
center  separation  vector  One  still  has  all  the  difficulties  of  ambiguous  effects  as  we  shall  discuss 
for  clustered  arrays  This  effec^  can  not  be  used,  however,  if  the  two  vectors  are  orthogonal,  as 
indicated  in  Figure  7-2 o 

The  results  for  noise  spectra  with  smooth  frequency  wave  number  functions  are  straight- 
forward if  there  is  no  separation  we  have  a power  spectrum  response 

2 2 * 

|g|(<ok-a])|  +|g2(wk-a2)|  + 2Re  1 g,(co  k-a])g2(w  k-a2)] 

while  if  there  is  a separation  we  have 

2 2 T J!S‘(El-B2) 

|g|(w.k’aj)|  + lg2(<o  k-a2)|  +2Rele  gj(w 

jk-(Ri-R7) 

If  the  space  factor  e reduces  the  overlap  volume  of  the  last  term  on  the  right, 

then  there  is  a smaller  amount  of  noise  propagating  through  the  array.  This  is  generally  the 
case  Physically,  with  the  greater  sepaiation,  the  noise  field  decorrelates  Some  typical 
optimum  beam  pattern  responses  for  separated  center  arrays  are  shown  in  Figures  7-27  and 
7-28  These  should  be  compared  to  the  original  centered  array  results  (See  Figures  7-16 
and  7-17.) 
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Figure  7-25.  Conventional  beam  pattern  for  a crossed  airay  with  separated  centers. 
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7.3  CIRCULAR  ARRAYS 


The  analysis  of  circular  arrays  is  similar  in  principle  to  that  of  linear  arrays,  however,, 
the  transcendental  functions  appropriate  for  their  analysis  are  not  familiar  to  many  Conse- 
quently, the  salient  features  of  circular  array  properties  are  often  clouded  due  to  the  appearance 
of  these  functions  Before  discussing  how  the  response  patterns  can  be  controlled  to  conhat 
various  noise  fields,  we  examine  the  class  of  beam  patterns  that  can  be  generated  by  typi-al 
circular  array  structures 

Ring  Arrays 

A ring  array  is  the  simplest  circular  geometry,  as  indicated  in  Figure  7-29  For  this  array 


"0 

■A 


jk-Rar(«) 

G(co  0)  e Rd0 


jkfR  cos(0-^k) 

G(w0)  e Rd0,. 


kr=  V(^l)Z  + (^'22)2  (7.311 

0k  = tan'1  (k-a,,  k-a2)  (7.31( 

We  expand  G(co  0),  the  aperture  weighting,  in  order  to  identify  the  structure  of  the  beam 
pattern 


E °»<' 


if' 

2tr  J 


-jn0 

gn(w.0)e  d 0 


Substituting  this  into  Eq  7 31a  for  gtw  k)  yields 
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v2'  p"  j(k,R  005(0-^1.)  + n0) 

o.y  = 2_j  Gn(<o)R  / e 

n=-oo  0 

jn(0k  - P!f  j(krR  sm0'  + n0’) 

= R /Gn(u)c  2 I e ‘ 

n-oo  •'o 

= 2irR  2_,  Gn(w)e  v 2/JnfkrR) 


We  see,  for  example,  that  for  a uniform  weighting  of  ( 2jtR)“  * we  have 
9(co.k)  = J0(krR), 

which  for  a ring  array  plays  a role  similar  to  the  smc  ^ka  ~j  for  linear  arrays  The  functions 

J0(x)  and  sincfx)  are  sketched  in  Figure  7-30.  We  point  out  that  J0(x)  has  its  first  zero  at 
2 4 and  a first  sidelobe  level  of -.402  occumng  at  3.8,  as  compared  to  tr,  and  -.212  at  4 7 
for  sinefx) 

The  effects  of  the  higher  order  terms  can  be  quite  complex  since  they  must  be  added 
in  phase  The  weighting  for  each  is  Jn(krR),  and  the  first  five  members  of  the  series  are 
indicated  in  Figure  7-31.  To  illustrate  the  effects  of  the  various  terms  we  consider  the  n and 
-n  terms  in  ihe  senes.  We  have 

2*R  J„(krR)  • 2Re[Gn(<o)ej,lWk~’r/2)]  (7  34) 

This  produces  a sinusoidal  pattern  versus  0k.  Tile  total  beam  pattern  consists  of  a complex 
linear  combination  of  these  terms  The  weightings  have  an  amphtude  dependence  upon  k_R 
determined  by  the  term 


Consequently,  the  important  terms  within  the  region  of  propagating  signals  are  determined 
by  the  magnitude  of  Jn(2trR/X),  which  decreases  for  increasing  n,  the  order  of  the  Bessei 
function. 
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One  can  intuitively  see  how  superdirective  effects  can  be  produced  at  a specific  kr^, 
by  choosing  the  coefficients  needed  to  produce  an  impulsive  function  versus  0^,  1 e , one  can 
ignore  the  dependence  upon  the  wave  number  magnitude.  We  have 


J 

2irRJn(krTR)e 


— -i^kT 
2G„(co)  = e 


so  that 


g(to  k) 


-Jn(«k^kT' 
e ' 


(7  35a) 


(7  35b) 


which  for  kr  = kr^  approximates  an  impulse  at  As  the  value  of  kr  is  changed  the  exact 
phasing  necessary  for  the  construction  of  the  impulse  altered  which  leads  to  extremely  large 
sidelobe  levels.  If  the  white  noise  level  is  finite  or  there  are  signals  propagating  with  a differ- 
ent wave  number  magnitude,  this  can  introduce  a serious  degradation  in  performance.  An 
illustration  of  these  effects  is  given  in  Figure  7-32  We  observe,  consequently,  that  for  any 
particular  kr,  we  can  achieve  a superdirective  beam  pattern,  however,  we  encounter  a sensitiv- 
ity problem  in  doing  so  In  any  finite  element  array,  we  cannot  achieve  the  impulsive  beam 
pattern;  however,  the  basic  arguments  are  similar.  In  reference  13,  a detailed  analysis  of 
circular  arrays  with  a finite  number  of  elements  has  been  done.  Consideration  was  given  to 
the  effects  of  sensitivity 

There  are  several  useful  relations  dealing  with  integral  representations,  recurrences 
and  series  expansions  for  Bessel  functions  Reference  1 1 is  one  of  the  readily  available 
sources  which  has  tabulated  many  of  them  We  can  illustrate  the  use  of  several  of  them  by 
considering  an  example  where  we  phase  the  array  to  a wave  number  with  finite  kp. 

For  a signal  with  wave  number  kp  we  need  a phasing  given  by 


1 -Jk-r  i -jkr  R cos(i ) 


■ c ‘ = e 

2trR  2irR 


T _ 


= G («-*),• 


where 


1 M 

krT=  l(kT-a|)2  + (kT-a2)2] 


(7.36) 
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Substituting  this  mto  Eq  7 31  yields 


g(co.k)  = J0[(k-kT)rR],. 


(k-kT)r  = t((k-kT)-a, )2  + ((k-kT)-a2)2] 1/2  (7  3' 

i e 3 we  simply  steer  the  array  to  ky  as  expected  For  example,  when  both  the  target  and 
wave  number  are  in  the  same  plane  as  the  array  with  wave  number  magnitude  of  27 r/X,  as 
illustrated  tn  Figure  7-1  b,  we  have 


g(oj  k)  = J 


(47rR  /A0\\ 

T~  sm\T// 


where  A 0 is  the  separation  angle  between  the  target  wave  number  and  h Naturally,  there  is 
no  change  as  the  target  moves  around  at  edge  because  of  the  circular  symmetry 

If  we  pursue  the  Founer  analysis,  we  arrive  at  the  same  result  in  a somewhat  more 
circuitous  manner  The  analysis  is,  however,  illustrative  of  the  use  of  the  relations  needed 
in  the  analysis  of  more  complex  problems  We  have 


i f* 

3"(“)=  5 J 


-jkrTR  cos(<H>kT)  -jn0 


-jn(v§) 


Substituting  this  into  Eq.  7 33  yields 


« jn(0k-0kT) 

g(wk>=  2_j  e Jn<krTR>  Jn<krR> 


One  now  uses  the  relation  [29,  Eq  9.1  79] 


Jn<W)P  = £ Jn+m(u)Jm(v)e,m“ 
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Figure  7-33.  Relation  of  terms  in  equation  7 40 

where  the  terms  W,u,v,  and  afi  are  geometrically  related  as  indicated  in  Figure  7-33.  This  yields 

g(cok)  = J0[(k-kT)rRj  (7  4!) 

As  pointed  out  above,  this  approach  is  not  as  direct,  however,  it  does  indicate  the  usefulness 
of  one  of  the  product  theorems  for  Bessel  functions.  Generally,  often  complicated  expressions 
can  be  simplified  considerably  by  using  these  relationships  or,  more  important,  there  is  a more 
straightforward  approach  to  the  analysis 

We  now  consider  the  effects  of  staving  a circular  array  whose  geometry  is  illustrated  in 
Figure  7-34  We  assume  that  the  staving  operation  multiplies  the  shading  G(co  <t>)  by  a gatmg 
function  Gst(0)  as  illustrated  in  Figure  7-35,  where  N is  the  number  of  staves,  <j>$  is  their  angular 
width,  and  <t>Q  is  the  origin  The  Fourier  senes  associated  with  this  function  is  given  by 
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Figure  7-35  Gating  function  for  representing  staved  arrays 
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The  effect  of  the  staving  can  be  represented  as  a gating  of  the  original  pattern,  i.e  , 


Gs(co  0)  — G^iQ)  G(o;  i 


Consequently,  the  coefficients  {Gst(n)  j and  |G(w  n) } , which  represent  the  Fourier  series  of 
the  shading  on  the  circle,  must  be  convolved  to  obtain  the  final  set  of  coefficients  for  the  array 
If  the  weighting  factor  has  significant  harmonics  greater  than  ^ , then  the  coefficients  are  aliased 
Generally  this  is  not  the  case,  the  only  effect  of  the  staving  is  to  cause  the  Fourier  coefficients 
to  be  repeated  at  intervals  of  N The  net  beam  pattern  is  then  (assume  <f>Q  = 0) 


GO  OO  i,  ^ 

g(to:k)  * N0wR  ^ G(co:m+n)el(m+nN>  V k 2/Jm+nN(krR) 

n =-“  m=-<” 

/nN0w\ 

• sinc(  ~ 1 (7.44) 

which  introduces  the  same  type  of  grating  structure  that  we  observe  for  clustered  hnear  arrays 
with  a periodic  repetition  interval  If  a highly  directional  source  can  propagate  through  one  of 
the  ambiguous  sidelobes  we  will  experience  deterioration  in  performance  If  not.  we  would  not 
be  bothered  except  for  the  increase  in  white  noise  due  to  the  increased  array  area.  The  simplest 
example  of  these  effects  can  be  demonstrated  when  the  array  is  conventionally  steered.  In  this 
case,  we  have 

G(cu  n)= — 1 — 6„o  , (7  45) 

N«SR  nU 

the  resulting  pattern  is  given  by 


jnN(0ic+7r/2)  , nN^s 

e K JnN(krR) smc  — 


Exact  statements  regarding  the  number  of  significant  sidelobes  are  difficult  to  make.  The 
question  becomes,  for  those  values  of  n where  sinc(nN0  s jj)  is  significant,  1 e , n < 2ir/0  s N, 
how  many  functions  Jnjj(krR)  have  maximum  values  in  the  wave  number  regio-i  of  interest 
This  question  requires  a rather  lengthy  examination  of  the  behavior  of  Bessel  functions. 

In  examining  the  design  of  optimum  beam  patterns  for  circular  arrays  we  confine  our 
detailed  analysis  to  directional  noise  fields  and  to  thost  which  are  isotropic  in  the  plane 
containing  the  array  The  general  concepts  illustrated  within  the  context  of  hnear  arrays 
are  still  applicable  and  are  quite  useful  in  analyzing  more  comphcated  noise  field  structures. 


As  with  crossed  arrays,  we  assume  that  there  is  sensor  noise  present  This  noise  has  a 
correlation  given  by 

E[w(t  tf>i)w (t.02)]  = N0S(t-r)5(p]-^2)  = NQ6(t-r)5 1 -02J  (7.47) 

Note  that  the  noise  output  generated  by  such  a process  is  given  by 


) = N° f 


a (<o)  = N0  / |g(co  0)|  Rd 0 

w 


(7  48) 


For  a conventionally  steered  array,  this  yields 


2 

o (w) 
w 


= No/ 


'2ir  jkpjR  cos(^x) 


2rrR 


Rd0  = No/2nR  = No/Af2  (7.49) 


Note  that  we  cannot  model  the  noise  with  a two-dimensional  white  noise  spectrum 
of  the  form 


E(w(t  z)w(r  f)l  = N0S(t-r)«(z-f) 


l i 50a) 


Pw(cok)  = N0 


(7  50b) 


The  reason  for  this  is  that  although  the  circular  array  has  resolution  capabilities  in  all  directions, 
its  wave  number  response  decays  asymptotically  only  as  |k|~^  When  integrated  across  the  entire 
two-dimensional  wave  vector  domain,  this  decay  is  not  sufficiently  fast  in  the  noise  power  output 
of  the  array  so  that  any  nonzero  aperture  weighting  is  necessarily  infinite 

In  Section  6 we  analyzed  optimum  beam  pattern  design  for  directional  noise  field 
with  an  arbitrary  number  of  components  and  arrays  of  a general  structure  We  can  easdy 
specialize  these  results  for  circular  arrays.  For  example,  for  a single  noise  source  we  have 


2irRSn(<o) 


}Jo([k-kT  ] R)-  - 


g(co.k  kT)  = 


27rRS  (to)  ^-T~-nlr^  Jo^--kn] rR) 


2irRSn(o>) 

2*RS°M"  J!(!-T'-n'rR) 
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°o(OJ^^=  jS 


2irRSn(o>)  , 2 \" 

N l+^T  ('-'o^T-WrR)) 


2irRSn(to) 


2 xz 

J (x)  =»  1 - — 
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I x | « 1 


sinc2(x)  “ 1- 


|x|«l. 


a circular  array  has  a 3/2  high  resolution  improvement  over  a linear  array  of  comparable 
dimensions  operating  at  broadside  In  Figures  7-36  and  7-37  typical  noise  power  outputs  for  a 
circular  array  for  directional  noise  fields  incident  normally  and  tangentially  are  indicated 

For  the  circular  array  we  can  analyze  the  optimum  processor  for  a noise  field  which 
has  a two  dimensional  isotropic  structure.  Our  representation  of  two-dimensional  noise 
processes  in  Section  2 and  some  useful  relations  among  Bessel  functions  are  the  key  factors 
which  produce  a relatively  simple  result  for  the  optimum  beam  pattern  and  its  noise 
power  output  From  Section  2,  Eq.  2 70,  we  have 


S"(“  'ZJ)=^[ 


p0(w  kr)J0(krjz|)krdkr  = — f0(o>  |z|) 
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for  the  structure  of  an  isotropic  two  dimensional  noise  field  From  Section  4,  Eq  4 9, 
we  want  to  solve  the  integral  equation 


n(w:  |z  - { | ) G(cj'£)d£  + N0G(w  z)  = X(w)e  zefi,. 


where  X(w)  is  chosen  such  that  the  beam  pattern  has  unity  response  in  the  target  direction.  For 
the  circular  array  geometry  this  becomes 


I s'(“l2R“"(^)l) 


G(co.0 2)  Rdtf>2  + 


jkj  R cos 

= X(w)e  f O<0]<2tr  (7  55) 

The  spatial  eigenfunctions  for  this  array  and  noise  field  are  ot  the  form  . This  suggests  using 
an  eigenfunction  approach,  which  for  this  case  is  simply  a Fourier  series  expansion  of  G(co  0),  1 e., 
we  expand  G(w  <j>)  m the  form 


OO 

G(to.0)=  Gn(oo)eJ" 


For  the  left  hand  side  of  Equation  7.55  we  have 


j ^ Sn^  ^Rsin  Gn(o,)eJn02Rd02  + No  J G^ule'"01 

= £ Gn(w)| J Sn^:|2Rsin^^e'Jn0Rd«  + No  e^'  (7  57) 


For  the  right  hand  side  of  Eq  7 55  we  use  the  Bessel  function  expansion  from  Ref  25 


This  is  the  continuous  analog  of  Gaardcr’s  observation  for  the  eigenvectors  for  discrete  arrays  / 131 


-'m  RS5«W.V»  t>:V.  .1 


gjZ  cos<p  = ^ ^ Jn(Z)eJn^+*/2> 


From  this  we  obtain 


X(u,)e'jkTrR  C0S(^)=  X(«)  2 J„(kT  R)eJn(’r/2-^)eJn01 


We  match  Fourier  coefficients  to  solve  forGn(oi),  or 


= X(ca)  ]T 


Jn(kTrR)eJn(!r/2^) 


J Sn(u:  1 2R  sin  ^e-Jn0Rd0  + Nq 


jnd) 

eJ  , (7.60) 


where  X(oj)  is  a normalization  factor.  The  term  in  the  denominator  can  be  simplified 
appreciably  by  use  of  the  wave  number  representation  of  Eq.  7.53.  Substituting  for 

Sn  (w:  |2R  sin  | ) yields 

/ Sn(w:  1 2R  sin  Rd«  = — f”  C p0(w:kr)J0(kr2R  sin(£) ) krRe'Jr 

•'0  ' ' (2 *kJq  Jq  '2' 


Again  one  of  the  Bessel  function  relationships  is  useful.  We  have 


i[  J°(2Zsin(?)) 


ejl^d0  = J2(Z); 
n 


therefore, 


£ Sn^r 1 2R sin  p0(«:kr)J2(krR)krd 


191 


In  order  to  compute  the  output  noise  we  need  from  Eq  4 1 5 


-jkT  R cosisHJj) 

o^(to  |kf)  = I / e r G(w:0)Rd0  (7.64) 


Using  tlie  Fourier  expansion  ofCKoi:^)  and  Eq  7.58  for  the  exponential  we  derive 
r / 7 


a\ to  k'j.)  = 2irR  ^ ^ 


J^(kTrR) 


if  P0<“:kr)JJkrR)Mkr+N0 


0^(10  kj)  = 


N^irR 

J2(kT  R) 
n *r 


PoCcu^JJ^krRJkrdkj+l 


A number  of  special  cases  can  be  considered  at  this  point.  For  example,  for  a target  normal  to 
array,  k[  is  zero  so  we  have 

°2(“'-T)=  5 j[  Po(oJ:kr)JJkrR)Vlkr+  >)  <7  »> 

\ e.,  only  the  term  n=0  is  significant  in  the  output  noise  power.  For  the  case  of  ring  noise,  from 
Eq  2 79,  we  have 


P0(“:k r> = T so<“)uo<kr-ko> 
Kr 


This  produces  an  output  power  of 


o^(oj.kj)  = 


N^irR 
J2(kT  R) 

n 1 r 


n=-“  \ N0RS0(u)  Jn(k0R)  + 


'••t.'-irsa ^i^w3vtf>s7s  t^CTirs  1 


Most  pfw:kr)  of  interest  are  zero  for  kj  > 2rr/X.  Since  Jn(x)  -*■  0 for  x < 2ir/X  with  n 
increasing,  there  are  a finite  number  of  terms  which  are  significant  This  is  analogous  to  linear 
arrays  when  there  were  L/X  + 1 significant  eigenvalues  Unfortunately,  one  must  resort  to 
numerical  procedures  to  evaluate  the  expressions,  just  as  in  the  linear  arras  Some  typical 
situations  are  indicated  in  Figures  7-38  and  7-39  As  the  final  topic  in  our  analysis  of  circular 
array  structures,  we  briefly  consider  some  aspects  of  disc  arrays  The  geometry  of  a disc  array 
is  the  same  as  a circular,  or  ring,  array  except  that  the  entire  area  is  filled  The  beam  pattern  of 
a conventionally  steered  disc  array  can  easily  be  found 


grc(co.k|kT)=  - — Qe 

An 


j(k-kT)-z 


dz 


This  is  often  termed  an  “Airy  Disc"  from  its  use  in  Fourier  optics  130) . In  Figure  7-40  we  have 
compared  the  beam  patterns  of  linear,  circular,  and  disc  arrays  with  the  same  length  and/or 
diameters.  We  observe  that  (2J  j(x)/x)2  has  its  first  zero  at  x = 3.8  and  a first  sidelobe  level  of 
0.0044  at  x = 5.2.  The  behavior  near  the  origin  is 

2J,(x)2  x2 

— st  1 - — >•  1x1  « 1 (7.69) 

Consequently,  its  resolution  is  approximately  the  same  as  a ring  array. 

In  analyzing  the  response  to  various  noise  fields  we  observe  that  directional  fields  are 
again  a special  case,  as  those  results  shown  in  Section  6,  and  that  analytic  results  for  isotropic 
fields  are  not  tractable  with  the  results  taking  the  form  of  Bessel  function  senes. 
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8.  ARRAY  PROCESSING  SYSTEMS 


In  many  situations  one  must  consider  systems  composed  of  arrays,  1 e.,  an  array  of 
arrays.  There  are  many  practical  reasons  that  lead  one  to  this  situation.  For  example,  one  often 
has  a discrete  set  of  sensors  the  apertures  of  which  have  a finite  set  of  geometries,  so  that  each 
sensor  may  have  directional  properties  in  itself.  If  the  individual  wave  number  responses  are 
more  or  less  omnidirectional  for  frequencies  of  interest  and  their  separations  are  large,  then 
the  theory  of  discrete  arrays  is  probably  most  appropnate.  Typically,  however,  the  sensors  are 
not  omnidirectional,  the  separations  are  modest,  and  we  need  to  consider  the  system  as  a whole 
The  physical  motivations  for  using  clusters  of  arrays  are  that  a complete  coverage  is 
either  too  expensive  or  impossible.  Further  economical  motivation  is  that  one  may  be  forced 
simply  by  high  data  rales  to  “cluster”  arrays  because  the  sensor  outputs  can  not  be  processed 
individually,  particularly  if  real-time  processing  is  important 

Usuaily  clustered  array  systems  have  a suboptimum  performance  so  it  is  important  to 
understand  when  the  losses  enter  and  to  compare  their  performance  to  an  optimized  system  In 
this  section  we  consider  the  analysis  of  suen  systems  and  consider  some  representative  examples. 
Figure  8-1  illustrates  the  general  structure  of  the  systems  being  considered. 
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Figure  8-1.  Clustered  array  structures. 

Several  possibilities  and  combinations  exist  in  the  (resign  of  these  systems  One  can 
conventionally  phase,  or  beamform,  with  or  without  shad.ng  each  subarray,  or  one  can  generate 
an  optimum  beam  as  we  have  done  throughout  the  text.  In  addition,  one  can  perform  direct  de- 
lay and  sum  beamforming,  or  one  can  introduce  the  theory  of  an  optimal  discrete  array  to  do 


.■'jmm&zttaMia 


in 


the  beam-combining  Consequently,  there  are  four  combinations  of  processing  that  we 
consider  We  have  tabulated  and  commented  upon  these  combinations  in  the  following  table 

Processing  combination  for  clustered  arrays 

Combinations  of  Beams 


Individual  Beams 

Conventional  Combination 

Optimum  Combination 

Conventional 

Arrays 

Conventional  array  theory  of 
sonar  texts 

General  area  of  application  of 
discrete  array  theory 

Optimum 

Arrays 

Possible  application  for  non- 
plane (distributed)  wave  targets 

Least  amount  of  performance 
loss  for  these  suboptimum 
systems 

There  are  also  two  considerations  which  commonly  simplify  the  details  of  the  analysis  con- 
siderably First,  each  of  the  individual  arrays  has  the  same  pattern,  except  for  a phase  offset  due 
to  the  displacement  of  their  centers,  in  the  second  arrays  are  regularly  spaced 

Since  we  have  considered  the  synthesis  of  individual  arrays,  we  concentrate  upon  the 
problem  of  combining  them.  'Ve  define  the  beam  pattern  of  each  array  for  a target  of  k-p,  as 
Gj(cj:k  |k-p).  The  output  of  each  array  is  given  by 


Jfj Gjfwklk- 


klkT)  dR  (w.k)  + dW;(to) 


Consequently,  the  spectral  noise  covariance  matrix  for  the  beams  is 


i(cjjk  kT)Pn(u,k)Gj  (w:k|kT) 


d K [ 2 

— +N06ij  J |Gj(co:z)|  dz 


We  constrain  the  individual  array  responses  to  satisfy 
G,(oj:kylkj)  — 1 ,. 

i.e.,  we  have  incorporated  the  steering  delay  in  the  individual  arrays.  There  is  no  loss  of 
generality  here  since  the  combiner  can  assign  an  arbitrary  weighting  to  compensate  for  this. 

The  beam  combiner  calculates  a weighted  sum  of  the  individual  array  outputs  In  suboptimum 
processing,  the  weights  are  preassigned  and  do  not  directly  depend  upon  the  noise  spectral  covariance 
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matrix  of  the  individual  array  outputs.  In  opamum  processing,  these  weights  are  optimized  to  achieve  *|l 
the  minimum  noise  output  power  From  discrete  array  theory  this  is  given  by  [ 30,  3 1 ) 


C,(co)  = 


|(‘olkT]t  [Sn(to)r'  eJkTZ| 


where 


IMkx]*  [Snico)]-1  E(to|kx) 
E[w  1 kT]  = [e'J-TZ  > , ,.  e"J-Tzn] 1 


and  z,  is  the  location  of  the  center  of  the  Nth  array 
In  either  the  optimum  or  conventional  case 


N 

dY0(<o)  = ]T  CjMdY^co), 

Pf 


HI 

g(<o.e|eT)  = y^  Cj(<o)  g,(<o  £|kx) 


(8.3) 


(8.4a) 


where  we  impose  the  constraint  of  unity  response  in  the  target  direction, 

N N 

Z Cj(co)Gi(w:K,kT)  = ^ Cj(w)«  1 
i=l  1^1 

The  beam  pattern  for  the  total  system  is 

M 

T, 

1=1 

and  the  performance  can  be  expressed  as  the  inverse  of  the  noise  output  power,  or 


(8.4b) 
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For  the  optimum  array  configuration  this  becomes 


Mw|kT)=  lSn(w)]“!  E[colkT]  j (8  7) 

which  is  identical  to  that  derived  for  discrete  element  arrays  This  is  as  far  as  one  can  go  with 
a general  formulation  for  an  arbitrary  collection  of  beams 

To  proceed  a step  further,  we  require  that  the  individual  arrays  have  a common  aperture 
weighting  except  for  a linear  phasing  for  its  arising  from  their  location,  or 


9,(w'K|kT)  = gcom(oi  k|kT)eJ(-  "t)  (8  8) 

Here  we  have  assumed  that  each  array  is  centered  at  z}  and  there  exists  the  possibility  of  only 
steering  the  individual  beams,  i.eM  the  weighting  is  constant  but  one  adds  a linear  phase  shift  to 
account  for  each  array  location 
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g0(w«  kT)  = 9Com(“  ^Cj(u>)e 
i=l 

= ?(w  klkT)gcomb(k-kT) 
where  Ocom  is  tiie  aperture  of  the  common  array 
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(8.9b) 
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Consequently,  we  find  that  the  final  beam  pattern  can  be  expressed  as  a product  of  the  beam 
pattern  for  the  individual  arrays  and  a beam  factor  for  the  combination  of  them  as  distributed 
in  space.  Note  that  the  C;(<o)  can  be  designed  either  conventionally  or  in  an  optimum  fashion. 

The  performance  is  expressed  in  basically  the  same  form  as  above  with  somewhat  minor 
simplifications. 


IV  IV 
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When  the  combmation  beam  is  designed  in  an  optimum  manner,  we  obtain 


0 (w|JiT) 


/ N N 
\ i=l  j=l 


[Sn..(co)l 

n« 


-1  ^T’Ji 
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= (ET(ca|kT)  [Sn(m)r>  Etwlkj))-'  (8.1 1) 

Essentially,  we  are  examinmg  the  covariance  for  a noise  field  which  has  passed  through  the  wave- 
number  filter  g(w:k  |kT).  There  is  an  extensive  literature  on  the  theory  of  optimum  discrete 
arrays.  Sometimes  the  relative  spacing  between  the  clusters  is  so  small  that  the  continuous 
array  theory  is  a more  convenient  analysis  tool  as,  for  example,  when  the  individual  array 
responses  are  omnidirectional  and  their  separations  quite  close  In  other  cases  the  spacing  is 
wide  and  the  continuous  theory  is  no  longer  appropriate  for  analyzing  the  design  and  operation 
of  the  beam  combiner  Essentially,  one  should  say  that  the  element  spacing  is  sparse 

The  analysis  of  sparse  arrays  is  quite  difficult  and  the  dominant  concern  is  the  sidelobe 
structure  of  the  beam.  There  are  two  special  cases,  however,  which  can  be  analyzed. 

Fortunately,  they  are  of  practical  interest  and  are  indicative  of  the  general  results  encountered. 
In  the  first  of  these  we  have  the  clusters  spaced  at  regular  intervals  m a lattice,  whereas  in  the 
second  the  individual  arrays  are  spaced  at  regular  intervals  on  rings.  The  analysis  of  lattice 
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transforms  has  been  studied  extensively  in  spectrographic  methods  for  the  analysis  of  crystal 
structures.  In  the  single  dimensioned  distributions,  the  theory  of  sampled  data  systems  and 
Z transform  methods  are  quite  useful.  The  analysis  of  ring  structures  is  more  complex  and  it 
is  easy  to  become  mixed  in  a maze  of  the  Bessel  functions  which  appear.  Fortunately,  there 
are  some  reasonable  simple  examples  which  indicate  the  basic  theory  involved 

A complete  discussion  of  lattice  array  would  take  us  too  far  astray  at  this  point  In 
particular,  a large  literature  exists  dealing  with  optimum  processing  for  discrete  arrays.  Most 
of  it  does  not  convey  the  intuitive  insight  that  we  have  obtained  for  continuous  arrays  We 
consider  line  clusters  which  have  a distribution  of  the  form  and  which  are  conventionally 
combined 

We  have  assumed  an  even  number  of  elements  The  modification  for  an  odd  number 
with  a cluster  located  at  the  origin  is  straightforward.  Substituting  into  Equation  8.9c  yields 
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Kk-aa)!(2n+l)/2]L 


This  is  in  a form  that  is  useful  for  defining  a term  which  puts  our  analysis  in  the  context  of 
Z transforms.  One  identifies  Z(k)  as 
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Figure  8-2.  Clustered  line  array. 


This  leads  to  the  expression 


c fkl  -j-'sa(N'l/2)L 

Gcomb(k>-e 


N-l 

Vc  N(«)e*‘8*Ln 

4-J  n-x 

n=0  2 


N-l 

= Z (y-(N-l/2)  V c N(w)Z(k)n  (8.14) 

JLmJ  n-~ 
n=0  2 

Tables  of  these  have  been  extensively  tabulated  in  conjunction  with  the  use  in  digital  filtering 
methods.  It  is  also  the  basic  formulation  used  by  Schelkunoff  in  his  studies  of  linear  arrays. 
First  we  note  that  Z repeats  itself  as  a function  of  k We  have 

Z(k)  = z(k  + n y 2a)  = z{k+n  f X'O  (815) 


202 


Consequently,  if  L/2X  is  greater  than  unity  then  the  beam  pattern  repeats  itself  within  the 
2jt 

region  of  |k|  < — or  where  propagating  noise  appears,  1 e.,  there  is  an  ambituity  in  the 

A 

beam  pattern.  We  examine  the  consequences  of  this  in  subsequent  discussions.  Basically,  the 
array  is  undersampled  at  this  point. 

In  order  to  make  our  discussion  more  concrete  we  consider  an  example;  let 


C»(“)  = N 


We  then  have 


Gcomb®  = (^)>HN-'/2)£  £ Z«n 


ZOO- (N— 1/2)  1 -Z(k)N  _ J.  Z(k)N''2  - Z(k )~N/2 
N 1-Z(k)  =N  Z(k)*/2  -Z(k)_*/2 


c jfk-oc,)  N -~J J sine  [(k-x.)  y] 


This  is  the  space  factor  for  a linear  array  weighting  of  the  individual  outputs,  as  illustrated 

in  Figure8-3.  The  main  lobe  of  the  pattern  array  factor  has  a width  of  • — — , i.e.,  it  is  still  determined 

NL 

by  the  total  array  length.  However,  there  exists  the  possibility  of  ambiguous  effects  due  to 
the  sidclob-  at  2jt/L.  If  this  is  outside  2rr/X,  then  only  white  noise  enters  through  this  side  lobe 
To  consider  the  effect  in  more  detail,  let  the  individual  arrays  be  phased  conventionally  and 
assume  that  their  length  is 
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In  Figure  8-4  we  have  sketched  the  product  G(cj.k|k-j-)  Ggom^k-kj)  which  assumes  the 
target  is  broadside  such  that  kj'«a  = 0.  We  make  a comparison  for  this  case  N = 4 There 
are  2k-l  side  lobes  of  level 


sine 


n = ± I , . . . , ±k 


One  consequently  obtains  a mam  lobe  of  width  2(ir/NL)  The  effect  of  the  sparse  sampling 
is  to  introduce  these  side  lobes  If  the  background  noise  is  over  regions  which  are 
significant,  then  one  has  a performance  aetenoration  due  to  the  sparse  sampling.  For  example, 
a directional  noise  sourse  could  enter  through  one  of  them.  If  the  noise  is  uniform,  or 

< 1/2  with  small  elements,  then  the  array  is  ot  sparse  and  one  does  not  observe  a reduction 
in  performance. 

While  we  have  chosen  a uniform  weighting  of  the  individual  array  outputs,  one  can 
use  any  type  of  weighting  desired.  Tables  of  Z transforms  are  useful  in  this  respect.  For  the 
more  common  shadings  one  observes  a spreading  of  the  width  2(n/NL)  and  a reduction  of 
the  minor  side  lobe  levels  between  the  beams. 

The  modification  for  the  situation  in  which  each  array  is  designed  optimally  is  straight- 
forward. One  still  observes  a modulation  of  the  main  beam  by  the  space  factor  However,  the 
more  important  issue  concerns  combining  the  beams  optimally,  which  is  the  basic  issue  in 
the  design  of  discrete  arrays. 

The  design  of  optimum  discrete  sparsely  spaced  discrete  arrays  is,  ip  gervnl,  best 
done  by  computational  methods  as  analytical  results  are  qu.te  tedious  to  obtain.  There  is, 
however,  one  example  which  does  not  involve  a great  deal  of  tedium  and  illustrates  some 
of  the  issues  which  appear  in  the  design  of  sparse  arrays.  As  one  might  expect,  the  noise 
field  consists  of  directional  signal  plus  a white  component.  We  assume  that  the  noise  has  a 
wave  number  kn  and  a level  Sn  («)  and  that  the  white  component  has  level  NQ.  From 
Equation  8.2  we  have 
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i<W":«|Ht>; 
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+VOJ)|GcomfoJhl-n)|2eJinUr2i)  =N°/' 


2 dk 

IGcom^yM  i <819> 

(2rr)N 
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Conventional 
Array  of  Length 


.vAiKfcV'*  n.s  vsa’Jvv 


We  define 


«-N  0f 


iGcomi“  bm  — 


* = V"),G“m<"  kn|kT)l 


the  matrix  whose  elements  are  defined  above  can  be  represented  as 


[S]  =«In  + (3 


e^-n " z 1 


£j-n"zN  I ej-n'-N 


The  inverse  of  the  matrix  can  be  found  using  the  inversion  formula  of  Section  4 


JV*n|  jej-n’zN  I 


The  optimum  noise  power  output  then  becomes 


oWtM- 

o Vet  l+/3N/a 

' i=I 


,K!sr!b>*3 


We  observe  that  the  space  factor  for  a uniformly  weighted  summing  of  the  individual  array 
output  enters  the  calculations.  The  same  comments  regarding  closed  form  expressions  for 
lattice  space  arrays  and  the  use  of  Z transform  theory  enter  here. 

If  we  consider  the  same  linear  regularly  spaced  distribution  discussed  above  with 
individual  arrays  of  extent  L„,  we  have 
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0 * Sn  (w)  sinc^ 
Ilo 


IV^  1 

Gcomb®  =l/ 


I 2. 


L 

sine 

^ 7 

sine 

L 1 

k-«a  yj 

2 LoN 


/L0N\  , Lv  / smc2(Ak-<yji) 

k** 


Ak  = krkn  (8.24e) 

This  is  very  similar  to  the  results  that  we  obtained  for  a continuous  array  of  length  L0N.  In 
Equation  8.24  we  need  only  identify  the  directional  noise  level  as  0 and  the  array  length  as 
L„N.  The  important  term  is  the  space  factor  which  is  identified  as  p.  This  causes  deterioration 
in  the  performance  at  wave  number  intervals  of 


(kT*n>«i  2 = v’ 


i.e.,  when  the  target  and  noise  are  separated  by  wave  number  components  of  2ir/L  as  projected 
upon  the  array.  At  these  points,  the  performance  becomes 


2 L0N  V*  - ,n, 

ao(“l*T>=TT  jl+  IT  Sno<“)Sinc2  K 


l0=k  n = ±l”  •’±k'1 
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Consequently,  if  > 1 , spatial  aliasing  enters  and  degrades  the  performance  The  only  way 
to  combat  this  is  for  the  individual  arrays  to  have  enough  resolution  to  eliminate  these  effects. 
In  Figure  8-5,  we  have  sketched  the  performance  for  some  representative  values  of  the 
parameters. 

This  concludes  our  discussion  on  clustered  arrays  The  principle  effects  are  space 
factors  and  spatial  aliasing  If  these  effects  are  not  significant,  then  clustering  represents  a 
possible  processing  method  which  has  distinct  computational  advantages,  especially  for 
optimum  processing.  We  have  confined  our  attention  to  linear  combinations.  While  the 
clustering  is  suboptimum  for  spatially  coherent  signals,  it  may  become  a superior  method  if 
the  signals  are  not  plane  waves.  Here  quadratic  operations  upon  the  clusters  may  become 
desirable. 


PREFORMED  BEAMS 

One  of  the  most  time  consuming  aspects  of  processing  array  data  for  changing 
noise  fields  is  in  calculating  and  forming  the  spatial,  or  beamforming,  operations  One  method 
of  reducing  the  amount  of  computation  is  to  preform  the  beams  and  operate  solely  upon  the 
beam  outputs.  F_.  —Tays  which  can  spatially  sample  the  ambient  field  on  the  lattice  dis- 
tribution of  points,  the  computational  advantages  are  even  more  attractive  since  fast  Fourier 
transform  methods  can  be  introduced  in  order  to  form  the  multiple  beam  outputs. 

Consider  an  array  processing  system  as  illustrated  in  Figure  8-6.  In  genera],  one  has 


HL 


r(r:f)gi(t-T:pdfdr 


The  combiner  performs  a linear  combination  upon  the  beam  outputs  so  as  to  direct  a sum 
beam  in  the  target  direction  kq,  with 


^ Cj(to:kT)  q;(oj:kT)  = 1 Vw 


For  conventional  processing  the  beam  spacing  is  quite  close  m wave  numbers,  so  that  the 
individual  beam  outputs  are  simply  scanned  to  look  in  the  respective  direction.  There  are  two 
areas  where  more  general  combinations  of  beams  are  useful:  in  null  placement  methods, 
which  are  closely  related  to  the  synthesis  of  optimum  pattern  for  directional  noises,  and  in 
the  approximate  synthesis  of  patterns  in  general.  As  the  two  are  closely  related,  we  consider 
them  simultaneously. 


209 


^ tw:kT)/Ajj 


U 
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r(T?)9i(t-T&<lJ  = yi(,)_ 

BEAM  # 1 


<{T-t)  92(l-r  f)  d?d r * y2(t) 
BEAM  #2 


Lf 


r(f.f)  sN(t-r.pdfdT=yN(t)_ 
8EAM#N 


y(t) 


The  basic  issue  is  the  selection  of  a set  of  coefficients  jCj(o>.kT)j 


N 

1=1 


g(u.k|kT)  = ^Cj(6J:kT)  g,-(«:k) 
i=l 


(8.29a) 


N 

G(co:z  |kT)  - } ' Cj(w:kj)  Gj(co.z) 
i=l 


(8.29b! 


such  that 


o2(w:kj)  =11  |g(w:k|kT)|  2 P(u:k)  -^=-  + N 0 / |G(«z|kT)|2  dz 
JJ  (2tt)n  Jn 


(8  30) 


is  minimized  subject  to  the  constraint  of  8.28.*  This  becomes 
N N 


WkT)  = TT.^WCti^W  f f fn(u:Sji(“:kW“:!i)"% 
M H JJ  J '27f>N 

+ N0  / Gj(w:z)Gf(w:z)dz 

Jn  f 


(8.3 1 a: 


We  are  temporarily  suppressing  the  \kj,  conditional 


and 


N 

^C,(<o:kT)  g;(«  kT)  = 
i=l 


(8.31b) 


To  optimize  the  choice  of  the  coefficients 


|cj((o.k-p)| 


i=! 


we  differentiate  incorporating 


the  constraint  by  Lagrange  multipliers.  By  differentiating  the  performance  with  the  constraint 
introduced  we  have  for  the  k**1  coefficient 


2 ^ ’ Sy(cj)  Cj*(cj)  + Xgj(co:kj)  = 0, 
j=l 


(8.32a) 


where  the  elements  of  the  matrix  [S^j(cj)]  are 


!ij(“)  =Jj PnCwtklg/ut^gKwik)  + NQJ  Gj(w.z)Gj(m.z)  dz 

mi: 


(8  32b) 


The  optimal  choice  of 


is  then 


C*(<o  kT)  = - — (S(oj)I  j]gj  («:kT) 


(8.33) 


i=l 


In  order  to  determine  the  Lagrangian  multiplier  X we  impose  the  constraint  equation  to 
obtain 


N N 


1 = ' | XaXa[S(‘J)ij!gi(“:-T)gj(‘j:-T)=-7  Tr([S(m)]pG(co  kx)),  (8  34a) 

i=l  j=l 


where 
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I 
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£o(  u:kT)  = 


g j(u:kj)|2 


9](o;  kj)g2(w:ilx^ 


g2(“  kT)g*(cj:kx)  | g2(tokT)|  2 

g,(cj:kT)g’(60  kT) 


(8.34b) 


As  a result,  we  have 


Cj(w  kj)  = 


N 4 

g.  (w:kx)  [S(to)]:1 
i=l 


Tr[pg(to:kT)  [S(w)l  '] 


The  system  performance  is  given  by 

jpg(<o:kT)  (S(co)l 


o2(co:kT)  = ^Tr  |pg(cj:kT)  (S(co)l  'jj 


(8  35) 


(8.36) 


Several  special  cases  are  of  interest  here  For  directional  noise  fields  with  a white 
background  component,  the  matrix  [S(w)l  is  evaluated  easily 


m r 

Sm<“)5i(“:!'m)9j*(“  !im)  + No/  C'i«o:z)G; 
m=l  TJ 


:z)Gf(co:z)  dz  (8.37) 


M 

E 

m=l 

This  is  further  simplified  if 

gm(«:kT)  = 0 m = 2JM 


i.e.,  all  but  one  of  the  preformed  beams  have  noils  in  the  target  direction.  The  performance 
becomes 


o2(ukT)  = [S(u)l  1 


with  a total  beam  pattern  of  the  form 

N 


(t  11  . 

g»  (oj:kT)y^tS(oj)l  jl1  gj(co:k)  (lS(w)]j|  J ' 
. 1 i=l  ' 


(8.38) 


This  is  shown  to  represent  an  optimum  null  placement  technique 
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As  an  example  of  this,  we  consider  the  following  Assume  that  we  have  an  array  of 

length  L and  that  there  are  noise  sources  of  level  Sn  (w)/2  located  symmetrically  at  broad- 

uo 

2it 

side,  kn  = ± -~cos(0),  plus  the  background  noise.  Becai  ' the  symmetry  of  the  situation 

A 

we  generate  preformed  beams  of  the  form 


for 


so  that 


G0(oj:C)  = ~ | fi|  < L/2 


G,(co  S)=-£- cos  2 < L/2 


1=  1,M-1 


g0(u:k)  = sinc|k-aa 
9 ](<o:k|  = sine  jk-aa  i-*J 

:hr'] 


(8.39a) 


(8  39b) 


(8  40a) 
(8  40b) 
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We  have  then 


soo(w) = sn0(w)  sinc2(,rLX  cos®n^ + Y" 


Soi(u>)  = si0(«)  = sn  (to)  smc2(rrL^  cos9n)  X [sincfaL^  cosOn-nr) 


- sincfrrL^  cos0n  + nr)) 


fori=  I,M-1, 


(8.41a) 


(8  41b) 
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Sjj(io)  = Sn^(cj)  (sinc(irL^  cos0n-nr)  + sincfaL^  cos0n+nr)) 


X [sinc(irL^  cos0n  + jir)  + sinc(irL^  cos0n-jjr)]  + — Sy, 


forij  = 1,M-1 
and 


[S(co)l 


-1 


(8.41c) 


LS„  (co) 
no 


N0  )l~~  N„ 


M-l 

sinc^(#L^cos9n)  + ^ ^(sind.7rL^cosfln-iir)  + sinc(jrL^cos0n+iir))2 
i=l 


sinc(jrL^cosfin) 

sinc(7rLjlcos9n) 

\ 

sinc(irL^cos8n-ir) 

sinc(irL^cos0n-ir) 

+ sinc(irL^cosdn+ir) 

• 

+ sinc(irL^cos0n+ir) 

1 

} 

• 

sinc(jiL^cosfln-(M-l  )ir) 

sinc(irL^cos0n-(M-l  )n) 

+ sinc(aL^cosan+(M- 1 )ir) 

+ sinc(7rL^cos0n+(M-I  )ir) 

) 

(8.42) 


The  only  beam  which  is  nonzero  at  broadside  is  the  Oth  one;  consequently  we  are  interested  in 
ISH  .or 
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I N„ 


/ M_l  V 

1 + HLS^M/N^  ^ sinc(iL^cos9n-i*)  + sinc(*Lxcosfln+iir)  j 
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We  see  that  the  term 


LSnQ(co)  NM 

— — / sinc(jrL^cos0n-iJr)  + sinc(trL^cos0n+i7r)2 

o Tl 


;-W 


determines  how  fast  we  can  approach  the  white  noise  performance  This  term,  however, 
represents  the  amount  of  noise  energy  that  can  be  cancelled  by  the  off-target  beams  One 
can  also  generate  a general  expression  for  the  beam  pattern. 

It  is  more  illustrative,  however,  to  consider  an  example 


“^n  mJ 


Then  we  have 


sinctri-^cosdj,)  = — 


sincOrL^cosdH-jr)  + sinc(7rL^cos0n-Hr)  = — 


",  4 2 

, + 97  T 

-J7  1+  7 

! 13 

+ 9*  7 


LSn0(w)  /2\2 
r=  N„  U 


The  beam  pattern  follows  as 
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1 i 4 L 

g(u-k  |kT)  = — — |(1  + -7)  sine  k^  - 


2 L LI 

-j7(sinc  k'^Y+^  + sinc  k*o^  — -ir)|  (846) 


We  see  that  when  y is  small  we  have  a conventional  beam;  however,  as  the  directional  noises 
become  strong  we  approach 


(o>;k  |kj)  “sinc^k'Og-jj- 'J  sine  ^k-aa  T*  ^ + slnc(-‘“a  \ ’7r)  ’ (8.47) 


which  places  a null  on  the  direction  of  the  noise  k = — cos 6n  = — • However,  as  we  have  seen 

A A 


in  our  discussion  of  superdirective  array,  it  introduces  severe  side  lobes  at  other  wave  numbers 
with  all  the  attendent  shortcomings  of  this  result.  The  performance  is  given  by 


a1  (cj:k|kf): 


% /i+?A 
L \i+4J 


When  7 is  small,  we  are  limited  by  the  white  noise  level  performance  of — . When  y is  large 

No 


the  performance  approaches 


4 L ’ 


i.e.,  we  are  willing  to  accept  13/4  increase  from  the  white  noise  in  order  to  cancel  the 
directional  noise  sources.  This  represents  the  additional  energy  that  enters  through  the  large 
side  lobes  in  the  beam. 

There  is  another  case  which  is  worthwhile  to  examine.  Assume  that  each  of  the  beams 
are  distributed  about  a wave  number  kj  and  that  the  beam  outputs  are  uncorrelated,  i.e., 


Sy(w)  = Sj(tj)8jj 


m 


to! 
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Such  a situation  could  appear  with  orthogonal  aperture  weightings  and  a frequency  wave 
number  spectrum  which  varies  slowly,  so  that  we  have 

Sy(<y)  = (P(«:ki)  + N0)/S2A«y 


i.e.,  each  beam  has  a power  level  proportional  to  the  frequency  wave  number  spectrum  at 
its  nominal  center  of  its  distribution.  The  optimum  processor  is  now  an  optimal  diversity 
ratio  combiner  with  the  weighting  determined  by  the  signal-to-noise  ratio  in  that  beam. 
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